[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN106357766A - Medical display device remote correction method based on lookup table rapid generation algorithm - Google Patents

Medical display device remote correction method based on lookup table rapid generation algorithm Download PDF

Info

Publication number
CN106357766A
CN106357766A CN201610828258.XA CN201610828258A CN106357766A CN 106357766 A CN106357766 A CN 106357766A CN 201610828258 A CN201610828258 A CN 201610828258A CN 106357766 A CN106357766 A CN 106357766A
Authority
CN
China
Prior art keywords
jnd
value
brightness
ddl
display
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610828258.XA
Other languages
Chinese (zh)
Other versions
CN106357766B (en
Inventor
李岳衡
薛团结
徐荣蓉
潘进勇
孙得娣
郭臣
孙蔓
任松
居美艳
黄平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201610828258.XA priority Critical patent/CN106357766B/en
Publication of CN106357766A publication Critical patent/CN106357766A/en
Application granted granted Critical
Publication of CN106357766B publication Critical patent/CN106357766B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L67/00Network arrangements or protocols for supporting network services or applications
    • H04L67/01Protocols
    • H04L67/02Protocols based on web technology, e.g. hypertext transfer protocol [HTTP]
    • H04L67/025Protocols based on web technology, e.g. hypertext transfer protocol [HTTP] for remote control or remote monitoring of applications
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L67/00Network arrangements or protocols for supporting network services or applications
    • H04L67/01Protocols
    • H04L67/02Protocols based on web technology, e.g. hypertext transfer protocol [HTTP]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L67/00Network arrangements or protocols for supporting network services or applications
    • H04L67/01Protocols
    • H04L67/12Protocols specially adapted for proprietary or special-purpose networking environments, e.g. medical networks, sensor networks, networks in vehicles or remote metering networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Digital Computer Display Output (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

本发明公开了一种基于查找表快速生成算法的医用显示设备远程校正方法,包括以下步骤:(1)利用互联网络通过http协议经加密处理后,实现远程指令发布,根据远程指令内容实现设备商服务器与医院服务器的连接;(2)根据Barten视觉模型和标准灰阶显示函数GSDF之间的关系,获得显示器亮度曲线;对标准化亮度显示的数字驱动等级DDL进行记录,通过查找表快速生成算法生成查找表格;(3)根据所得的查找表,通过标准灰阶显示矫正方法对问题显示器进行远程校正。本发明提高了医用显示设备校正效率,降低了医用显示设备的人工维护成本,有效降低了由于显示器显示不准确所带来的医生误诊,提高了医学数字影像系统的有效性和准确性。

The invention discloses a remote calibration method for medical display equipment based on a look-up table rapid generation algorithm, which comprises the following steps: (1) using the Internet to implement remote command release after encrypted processing through the http protocol, and realizing the device manufacturer according to the content of the remote command The connection between the server and the hospital server; (2) Obtain the brightness curve of the display according to the relationship between the Barten visual model and the standard gray scale display function GSDF; record the digital drive level DDL of the standardized brightness display, and generate it through a lookup table rapid generation algorithm Lookup table; (3) According to the obtained lookup table, remote correction is performed on the problematic display through a standard grayscale display correction method. The invention improves the correction efficiency of the medical display device, reduces the manual maintenance cost of the medical display device, effectively reduces the doctor's misdiagnosis caused by inaccurate display display, and improves the effectiveness and accuracy of the medical digital image system.

Description

一种基于查找表快速生成算法的医用显示设备远程校正方法A remote calibration method for medical display equipment based on a fast generation algorithm of look-up table

技术领域technical field

本发明涉及一种适用于复杂网络环境下的医用显示设备远程校正方法,具体涉及一种基于查找表快速生成算法的医用显示设备远程校正方法,属于医疗设备维护技术领域。The invention relates to a remote calibration method for medical display equipment applicable in a complex network environment, in particular to a remote calibration method for medical display equipment based on a lookup table rapid generation algorithm, and belongs to the technical field of medical equipment maintenance.

背景技术Background technique

医用显示设备管理与显示的标准化校正在新型医疗机构网络管理领域和医疗设备维护领域有着很强的现实意义。医用显示器的显示校正是指对医用显示设备的显示特征抽样出曲线特点(也即将显示设备的显示特征以曲线的形式表现出来),从而可以依据相关规范对显示设备显示特性进行检验,并进而对显示特性不符合标准特性的显示设备进行标准化显示校正。传统医用显示设备的校准都是在医用显示器出现故障报修之后,通过派遣技术人员到现场进行故障排除和显示校准。这种处理方式的缺陷很明显,即无法在发生故障时及时处理和预警,容易产生误诊并影响诊断准确率。The standardized correction of medical display equipment management and display has a strong practical significance in the field of new medical institution network management and medical equipment maintenance. The display calibration of the medical display refers to sampling the curve characteristics of the display characteristics of the medical display equipment (that is, the display characteristics of the display equipment are expressed in the form of a curve), so that the display characteristics of the display equipment can be tested according to the relevant specifications, and then the A display device whose display characteristics do not conform to the standard characteristics is subjected to standardized display correction. The calibration of traditional medical display equipment is to send technicians to the site for troubleshooting and display calibration after the medical display fails and is reported for repair. The defect of this processing method is obvious, that is, it cannot be processed and warned in time when a fault occurs, and it is easy to cause misdiagnosis and affect the accuracy of diagnosis.

发明内容Contents of the invention

针对现有技术存在的不足,本发明目的是提供一种基于查找表快速生成算法的医用显示设备远程校正方法,通过互联网远程监控医用显示器,并在出现故障时及时通过网络进行远程校正,不但大大提高了工作效率,同时有效降低了设备维护的人工成本。Aiming at the deficiencies in the prior art, the purpose of the present invention is to provide a remote calibration method for medical display equipment based on a look-up table rapid generation algorithm, remotely monitor the medical display through the Internet, and perform remote calibration through the network in time when a fault occurs, which not only greatly Improve work efficiency and effectively reduce the labor cost of equipment maintenance.

为了实现上述目的,本发明是通过如下的技术方案来实现:In order to achieve the above object, the present invention is achieved through the following technical solutions:

本发明的一种基于查找表快速生成算法的医用显示设备远程校正方法,包括以下几个步骤:A remote calibration method for medical display equipment based on a look-up table rapid generation algorithm of the present invention comprises the following steps:

(1)利用互联网络,通过http协议将设备商服务器发送的指令信号经加密处理(现有加密技术)后,发送到医院本地服务器以实现远程指令发布,医院本地服务器根据远程指令内容实现设备商服务器与医院服务器的连接,从而使得处在远端的管理员能够对医院的医用显示设备实现远程控制;(1) Using the Internet, the command signal sent by the equipment vendor server is encrypted through the http protocol (existing encryption technology), and then sent to the hospital’s local server to realize remote command release. The hospital’s local server realizes the equipment vendor’s The connection between the server and the hospital server, so that the remote administrator can remotely control the hospital's medical display equipment;

(2)根据Barten视觉模型和标准灰阶显示函数GSDF之间的关系,获得显示器亮度曲线;并对标准化亮度显示的数字驱动等级DDL进行记录,通过查找表快速生成算法生成查找表格;(2) Obtain the brightness curve of the display according to the relationship between the Barten visual model and the standard gray scale display function GSDF; and record the digital drive level DDL of the standardized brightness display, and generate a lookup table through a lookup table rapid generation algorithm;

(3)根据步骤(2)所得的查找表,通过标准灰阶显示矫正方法对问题显示器进行远程校正。(3) According to the look-up table obtained in step (2), remote correction is performed on the problematic display through a standard gray scale display correction method.

步骤(1)中,所述设备商服务器与医院服务器的连接方法如下:In step (1), the connection method between the equipment vendor server and the hospital server is as follows:

首先,设置本地协议地址和主机端口,通过TCP协议设定远程设备IP地址和端口号,创建套接字连接;一旦创建套接字,即可以通过Connect函数连接到指定的服务器并通过Send/SendTo函数向远端服务器发送数据,同时以Receive/ReceiveFrom函数从服务端接收数据,服务器端则需要使用Bind函数绑定所指定的接口使套接字与一个本地终结点相联,并通过Listen函数侦听该接口上的请求,当侦听到用户端的连接时,则调用Accept函数完成连接的操作,并创建新的套接字处理传入的连接请求;First, set the local protocol address and host port, set the IP address and port number of the remote device through the TCP protocol, and create a socket connection; once the socket is created, you can connect to the specified server through the Connect function and pass Send/SendTo The function sends data to the remote server, and at the same time receives data from the server with the Receive/ReceiveFrom function, and the server needs to use the Bind function to bind the specified interface to connect the socket to a local endpoint, and use the Listen function to detect Listen to the request on this interface, when listening to the connection of the client, call the Accept function to complete the connection operation, and create a new socket to handle the incoming connection request;

随后,调用Begin Receive函数用于异步接收数据,遍历Socket Flag接收到的所有数据和控制信息参数,并存储到指定偏移量位置的缓冲区,以便对接收到的指令进行解析;以解释为命令字符标记区分每条指令;根据不同的命令码,按照Telnet协议数据格式转换为相应的字符串;Subsequently, call the Begin Receive function to receive data asynchronously, traverse all data and control information parameters received by Socket Flag, and store them in the buffer at the specified offset position, so as to analyze the received instructions; interpret them as commands The character mark distinguishes each instruction; according to different command codes, it is converted into a corresponding string according to the data format of the Telnet protocol;

接收完数据后,根据取得的不同控制信息参数发送回返数据;根据实现的Telnet协议类得到回返数据和接收数据,然后由网络设备的命令参数实现对应的功能模块。After receiving the data, send the return data according to the obtained different control information parameters; obtain the return data and received data according to the implemented Telnet protocol class, and then realize the corresponding function modules by the command parameters of the network equipment.

步骤(2)中,所述显示器亮度曲线通过以下方法获得:In step (2), the display brightness curve is obtained by the following method:

符号j代表人眼可识别的最小可觉差JND级别,采用公式(1)来计算调制闭值S(j);The symbol j represents the minimum perceivable difference JND level recognizable by the human eye, and formula (1) is used to calculate the modulation closing value S(j);

SS (( jj )) == qq 11 ·· Mm oo pp tt (( LL uu mm )) qq 22 dd 22 LL uu mm ++ qq 33 -- -- -- (( 11 ))

式中,q1、q2、q3为调制系数,d为瞳孔直径,Lum为目标光照亮度,即光源在同一方向的光强与发光面在该方向上投影表面积之比;是高斯点扩展函数,其中,μ为空间频率,Csph为球面差,σ0为σ的参照值,σ是高斯点扩展函数的参数;In the formula, q 1 , q 2 , and q 3 are the modulation coefficients, d is the pupil diameter, and Lum is the target illumination brightness, that is, the ratio of the light intensity of the light source in the same direction to the projected surface area of the light-emitting surface in this direction; is the Gaussian point spread function, Among them, μ is the spatial frequency, C sph is the spherical difference, σ 0 is the reference value of σ, and σ is the parameter of the Gaussian point spread function;

人眼可识别的每一级亮度均值都是通过上一级的亮度均值计算而来,计算关系如公式(2)所示:The average brightness of each level recognizable by the human eye is calculated from the average brightness of the previous level, and the calculation relationship is shown in formula (2):

LL uu mm (( jj ++ 11 )) == LL uu mm (( jj )) 11 ++ SS (( jj )) 11 -- SS (( jj )) -- -- -- (( 22 ))

瞳孔直径d的计算公式为:The formula for calculating the pupil diameter d is:

d=4.6-2.8·tanh[0.4·log10(0.625·Lum)]。d = 4.6-2.8·tanh [0.4·log 10 (0.625·Lum)].

通过所述查找表快速生成算法生成查找表格的方法如下:The method for generating the lookup table by the fast generation algorithm of the lookup table is as follows:

(2-1)通过获得的显示器亮度曲线,导出数字驱动等级DDL与亮度关系表;(2-1) Deriving the relationship table between digital drive level DDL and brightness through the obtained display brightness curve;

(2-2)获取1023个JND值所一一对应的亮度值,导出最小可觉差JND与亮度关系表;(2-2) Obtain the luminance values corresponding to 1023 JND values one by one, and derive the relationship table between the minimum perceivable difference JND and luminance;

(2-3)依据最小可觉差JND与亮度关系表,求出每个最小可觉差JND值一一对应三次样条插值基值,同时导出最小可觉差JND对样条基值表JND-Spline,当JNDindex=n时,其中,JNDindex表示JND的序列值变量,n是JND的具体数值;设立中间变量{U(n),1≤n≤1022},其中,U(n)是构造的中间变量且U(0)=0,U(1023)=0,S(1023)=0,Y(0)=0,其中,Y(0)是中间变量Y(n)在时的数值,建立调制闭值S(n)数列递推公式,此处公式中的JND(n)代表最小可觉差JND的序列值n:(2-3) According to the relationship table between the minimum perceivable difference JND and brightness, each minimum perceivable difference JND value corresponds to the cubic spline interpolation base value one by one, and at the same time derive the minimum perceivable difference JND pair spline base value table JND -Spline, when JND index = n, wherein, JND index represents the serial value variable of JND, and n is the specific value of JND; set up an intermediate variable {U(n), 1≤n≤1022}, wherein, U(n) is the constructed intermediate variable and U(0)=0, U(1023)=0, S(1023)=0, Y(0)=0, where Y(0) is the intermediate variable Y(n) in When the numerical value is established, the modulation closed value S(n) sequence recursion formula is established, where JND(n) in the formula represents the sequence value n of the minimum perceivable difference JND:

QQ (( nno )) == JJ NN DD. (( nno )) -- JJ NN DD. (( nno -- 11 )) JJ NN DD. (( nno ++ 11 )) -- JJ NN DD. (( nno -- 11 )) -- -- -- (( 44 ))

P(n)=Q(n)*Y(n-1)+2 (5)P(n)=Q(n)*Y(n-1)+2 (5)

YY (( nno )) == QQ (( nno )) -- 11 QQ (( nno )) -- -- -- (( 66 ))

KK (( nno )) == LL uu mm (( JJ NN DD. (( nno ++ 11 )) )) -- LL uu mm (( JJ NN DD. (( nno )) )) JJ NN DD. (( nno ++ 11 )) -- JJ NN DD. (( nno )) -- LL uu mm (( JJ NN DD. (( nno )) )) -- LL uu mm (( JJ NN DD. (( nno -- 11 )) )) JJ NN DD. (( nno )) -- JJ NN DD. (( nno -- 11 )) -- -- -- (( 77 ))

Uu (( nno )) == 66 ** KK (( nno )) JJ NN DD. (( nno ++ 11 )) -- JJ NN DD. (( nno -- 11 )) -- QQ (( nno )) ** Uu (( nno -- 11 )) PP (( nno )) -- -- -- (( 88 ))

S(n)=K(n)*S(n-1)+U(n) (9)S(n)=K(n)*S(n-1)+U(n) (9)

其中Q(n),P(n),Y(n),K(n)是为了得到调制闭值S(n)数列所构造的中间变量;Among them, Q(n), P(n), Y(n), K(n) are intermediate variables constructed in order to obtain the sequence of modulated closed values S(n);

(2-4)根据数字驱动等级DDL与亮度关系表查找DDLmin、DDLmax对应的亮度值LumDDLmin、LumDDLmax,其中,DDLmin为设定定义范围内DDL的最小值,DDLmax为设定定义范围内DDL的最大值,进而由最小可觉差JND与亮度关系表获得LumAmb+LumDDLmin对应的最小可用最小可觉差JND值JNDmin,以及LumAmb+LumDDLmax对应的最大可用JND值JNDmax,LumAmb为测量得到的环境光亮度值,在此基础上获得显示对比度系数,计算公式如下:(2-4) Look up the brightness values Lum DDLmin and Lum DDLmax corresponding to DDL min and DDL max according to the digital drive level DDL and brightness relationship table, where DDL min is the minimum value of DDL within the defined range of the setting, and DDL max is the setting The maximum value of DDL within the defined range, and then obtain the minimum available minimum perceivable difference JND value JND min corresponding to Lum Amb +Lum DDLmin and the maximum available JND value corresponding to Lum Amb +Lum DDLmax from the table of the relationship between the minimum perceivable difference JND and brightness JND max , Lum Amb is the measured ambient light brightness value, on this basis, the display contrast coefficient is obtained, and the calculation formula is as follows:

DD. ii sthe s pp aa ythe y __ cc oo nno tt rr aa sthe s tt == JNDJND mm aa xx -- JNDJND mm ii nno DDLDDL mm aa xx -- DDLDDL mm ii nno -- -- -- (( 1010 ))

(2-5)将JNDmin到JNDmax的Barten视觉模型曲线分成JNDmax至JNDmin区间,再利用对样条基值表JND-Spline得到每个区间内三次插值多项式,然后在[JNDmax,JNDmin]内做2Pbits次等间隔取样,其中,Pbits为查找表的位宽,运用假设检验中常用的P值方法生成假定值P值对JND值的查找表P-JND;对于第i个P值所对应的最小可觉差JND序列值(2-5) Divide the Barten visual model curve from JND min to JND max into intervals from JND max to JND min , and then use the paired spline base value table JND-Spline to obtain a cubic interpolation polynomial in each interval, and then in [JND max , JND min ] to do 2 Pbits equal interval sampling, wherein, Pbits is the bit width of the lookup table, use the P value method commonly used in hypothesis testing to generate the lookup table P-JND of the assumed value P value to JND value; for the i-th The least noticeable difference JND sequence value corresponding to the P value

JJ (( ii )) == JNDJND mm aa xx -- JNDJND mm ii nno 22 PP bb ii tt sthe s ** (( ii -- 11 )) ++ JNDJND mm ii nno ,, 11 ≤≤ ii ≤≤ 22 PP bb ii tt sthe s -- -- -- (( 1111 ))

确定了区间n,亮度值就可以通过J(i)所在最小可觉差JND值的区间n来求得,也即它的对应亮度值可以通过公式(12)计算出:Once the interval n is determined, the luminance value can be obtained through the interval n of the minimum perceivable difference JND value where J(i) is located, that is, its corresponding luminance value can be calculated by formula (12):

LL uu mm (( JJ (( ii )) )) == (( JJ NN DD. (( nno )) -- JJ (( ii )) )) ·&Center Dot; LL uu mm (( JJ NN DD. (( nno )) )) ++

(( JJ NN DD. (( ii )) -- JJ NN DD. (( nno -- 11 )) )) ·&Center Dot; LL uu mm (( JJ NN DD. (( nno -- 11 )) )) ++ SS (( nno )) 66 ·&Center Dot; [[ (( JJ NN DD. (( nno )) -- JJ (( ii )) )) 33 -- JJ NN DD. (( nno )) ++ JJ (( ii )) ]] ++ SS (( nno -- 11 )) 66 ·&Center Dot; [[ (( JJ (( ii )) -- JJ NN DD. (( nno -- 11 )) )) 33 ++ JJ NN DD. (( nno -- 11 )) -- JJ (( ii )) ]] -- -- -- (( 1212 ))

(2-6)根据最小可觉差JND与亮度关系表和DDL与亮度关系表找到每一个P值对应的DDL值,从而获得据最小可觉差JND值和亮度查找DDL值的查找表。(2-6) Find the DDL value corresponding to each P value according to the relationship table between JND and luminance and the relationship table between DDL and luminance, so as to obtain a lookup table for finding the DDL value according to the JND value and luminance.

步骤(3)中,通过标准灰阶显示矫正方法对问题显示器进行远程校正具体包括以下几个步骤:In step (3), the remote correction of the problematic display through the standard grayscale display correction method specifically includes the following steps:

(3-1)使用亮度计量设备采集得到不同DDL值所对应的亮度值,然后获得需要矫正的医院本地显示器设备的特征曲线,也即是设备亮度随DDL变化而变化的曲线(类似于图2);(3-1) Use brightness metering equipment to collect brightness values corresponding to different DDL values, and then obtain the characteristic curve of the local display device in the hospital that needs to be corrected, that is, the curve of the device brightness changing with the change of DDL (similar to Figure 2 );

(3-2)标准灰阶显示函数GSDF曲线被定义为12bit亮度级别数学插值,每个级别的亮度值都能够根据Barten视觉模型得到;通过Barten视觉模型中的关系计算公式(Barten视觉模型中的关系计算公式已被公认,此处不再赘述),可以计算出不同的JND级别值所对应的不同亮度值,其所对应亮度如下:(3-2) The standard grayscale display function GSDF curve is defined as a mathematical interpolation of 12bit brightness levels, and the brightness value of each level can be obtained according to the Barten visual model; through the relationship calculation formula in the Barten visual model (Barten visual model in The relationship calculation formula has been recognized and will not be repeated here), and different brightness values corresponding to different JND level values can be calculated, and the corresponding brightness values are as follows:

loglog 1010 LL uu mm (( jj )) == aa ++ cc ·&Center Dot; LL nno (( jj )) ++ ee ·&Center Dot; [[ LL nno (( jj )) ]] 22 ++ gg ·&Center Dot; [[ LL nno (( jj )) ]] 33 ++ mm ·&Center Dot; [[ LL nno (( jj )) ]] 44 11 ++ bb ·&Center Dot; LL nno (( jj )) ++ dd ·· [[ LL nno (( jj )) ]] 22 ++ ff ·· [[ LL nno (( jj )) ]] 33 ++ hh ·· [[ LL nno (( jj )) ]] 44 ++ kk ·· [[ LL nno (( jj )) ]] 55 -- -- -- (( 33 ))

上式中,Ln表示取自然对数;j表示人眼可识别的最小可觉差JND级别,取值范围为从1至1023的自然数;a、b、c、d、e、f、g、h、k、m均为公式计算所需的常数;In the above formula, Ln represents the natural logarithm; j represents the JND level of the minimum perceivable difference that can be recognized by the human eye, and the value range is a natural number from 1 to 1023; a, b, c, d, e, f, g, h, k, and m are constants required for formula calculation;

(3-3)由以上各个对应关系,能够导出JND值和DDL之间的关系,然后以JND值与亮度值呈线性关系为目标进行相互补偿(现有方法,此处不再赘述);(3-3) From the above correspondences, the relationship between the JND value and the DDL can be derived, and then the JND value and the luminance value are linearly related to each other for mutual compensation (existing methods, not repeated here);

(3-4)通过对比调整前后的本地显示器设备特征曲线的相似度来评估调整符合度,其相似度越高,调整符合度就越好;反之则不好,需要继续调整(评估方法已经存在,此处不再赘述)。(3-4) Evaluate the adjustment compliance by comparing the similarity of the local display device characteristic curve before and after adjustment. The higher the similarity, the better the adjustment compliance; otherwise, it is not good and needs to continue to adjust (the evaluation method already exists , will not be repeated here).

本发明提高了医用显示设备校正效率,有效降低了医用显示设备的人工维护成本,最重要的是有效降低了由于显示器显示不准确所带来的医生误诊,提高了整个医学数字影像系统的有效性和准确性。The present invention improves the calibration efficiency of medical display equipment, effectively reduces the manual maintenance cost of medical display equipment, and most importantly effectively reduces the doctor's misdiagnosis caused by inaccurate display display, and improves the effectiveness of the entire medical digital imaging system and accuracy.

附图说明Description of drawings

图1为基于互联网技术的医用显示器远程校准系统原理架构图;Figure 1 is a schematic diagram of the remote calibration system for medical displays based on Internet technology;

图2为医用显示器亮度与最小可觉差之间关系曲线图;Fig. 2 is a graph showing the relationship between the brightness of a medical display and the minimum perceivable difference;

图3为标准灰阶显示函数GSDF算法校准过程图;Figure 3 is a calibration process diagram of the standard grayscale display function GSDF algorithm;

图4为本发明医用显示设备标准校正流程图。Fig. 4 is a flow chart of the standard calibration of the medical display device of the present invention.

具体实施方式detailed description

为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。In order to make the technical means, creative features, goals and effects achieved by the present invention easy to understand, the present invention will be further described below in conjunction with specific embodiments.

医用显示器的显示校正是指对医用显示设备的显示特征抽样出曲线特点,从而可以依据明确规范对显示设备显示特性进行标准验证,从而对显示特性不符合标准特性的显示设备进行标准化显示校正。The display calibration of medical displays refers to sampling the curve characteristics of the display characteristics of medical display equipment, so that the standard verification of the display characteristics of the display equipment can be carried out according to the clear specifications, so as to standardize the display calibration of the display equipment whose display characteristics do not meet the standard characteristics.

参见图4,本发明提出了一种医用显示设备远程校正方法。首先,利用网络通过http协议经加密处理后实现远程指令发布,根据指令内容对医用显示器实现远程控制;其次,本地校准部分根据Barten视觉模型对人的视觉特点的阐述和标准灰阶显示函数理论,引入人眼最小可觉差的概念,对标准显示和眼球感知特点之间的函数关系进行推导,使亮度的均匀变化在人眼感知中呈现线性变化,并将亮度显示的DDL进行记录,绘制出一个最小可觉差与显示器二极管亮度驱动电压的表格;最后,根据该查找表格完成显示器的远程标准显示校正。Referring to Fig. 4, the present invention proposes a remote calibration method for medical display equipment. First, use the network to release remote commands after encrypted processing through the http protocol, and realize remote control of the medical display according to the content of the commands; secondly, the local calibration part is based on the elaboration of the human visual characteristics of the Barten visual model and the standard grayscale display function theory, Introduce the concept of the minimum perceivable difference of the human eye, deduce the functional relationship between the standard display and the characteristics of eyeball perception, so that the uniform change of brightness presents a linear change in human perception, record the DDL of the brightness display, and draw the A table of minimum perceivable difference and display diode luminance driving voltage; finally, the remote standard display calibration of the display is completed according to the look-up table.

本发明的实施例可广泛应用于医疗机构显示设备的远程校准系统中,即通过设计一种远程互联网络连接和管理员本地校正的方法,对远端的医用显示设备进行指令操作,并运用快速查找表算法实现远程显示设备的本地快速校正。本发明的发明思想重点在本地的校准算法部分,也即提出了一种基于查找表快速生成算法的医用显示设备远程校正方法,以方便设备商管理人员在本地对医疗机构网络内的问题显示设备进行远程诊断与校正。The embodiment of the present invention can be widely used in the remote calibration system of display equipment in medical institutions, that is, by designing a method of remote Internet connection and local calibration by administrators, the remote medical display equipment can be commanded and operated, and the fast A look-up table algorithm enables local fast calibration of remote display devices. The inventive idea of the present invention focuses on the local calibration algorithm, that is, a remote calibration method for medical display equipment based on a look-up table rapid generation algorithm is proposed, so as to facilitate the management personnel of the equipment manufacturer to locally display the equipment for problems in the medical institution network Carry out remote diagnosis and correction.

一般医用显示设备在医疗机构内的网络结构及其远程互联校正控制示意图如图1所示:医用显示器位于医院局域网内,可以通过互联网与处于远端的设备商管理员进行连接。具体连接可通过Telnet协议实现,Telnet协议是TCP/IP协议族中的一员,是Internet远程登陆服务的标准协议和主要方式。The schematic diagram of the network structure and remote interconnection correction control of general medical display equipment in medical institutions is shown in Figure 1: the medical display is located in the hospital LAN and can be connected to the remote equipment vendor administrator through the Internet. The specific connection can be realized through the Telnet protocol, which is a member of the TCP/IP protocol family and is the standard protocol and main method of Internet remote login service.

首先,设置本地协议地址和主机端口,通过TCP协议设定远程设备IP地址和端口号,创建套接字连接。一旦创建套接字,即可以通过Connect函数连接到指定的服务器并通过Send/SendTo函数向远端服务器发送数据,同时以Receive/ReceiveFrom函数从服务端接收数据,服务器端则需要使用Bind函数绑定所指定的接口使套接字与一个本地终结点相联,并通过Listen函数侦听该接口上的请求,当侦听到用户端的连接时,则调用Accept函数完成连接的操作,并创建新的套接字处理传入的连接请求(以上所提及的函数均是Linux C下的功能函数)。随后调用Begin Receive函数用于异步接收数据,遍历Socket Flag接收到的所有数据和控制信息参数,并存储到指定偏移量位置的缓冲区,以便对接收到的指令进行解析。以解释为命令字符(IAC,Interpreted As Command)标记区分每条指令;根据不同的命令码,按照Telnet协议数据格式转换为相应的字符串。接收完数据后,根据取得的不同控制信息参数发送回返数据。根据实现的Telnet协议类得到回返数据和接收数据,然后由网络设备的命令参数实现对应的功能模块。First, set the local protocol address and host port, set the IP address and port number of the remote device through the TCP protocol, and create a socket connection. Once the socket is created, you can connect to the specified server through the Connect function and send data to the remote server through the Send/SendTo function, and at the same time receive data from the server through the Receive/ReceiveFrom function, and the server needs to use the Bind function to bind The specified interface connects the socket with a local endpoint, and listens to the request on the interface through the Listen function. When the connection of the client is detected, the Accept function is called to complete the connection operation and create a new one. The socket handles incoming connection requests (the functions mentioned above are all functional functions under Linux C). Then call the Begin Receive function to receive data asynchronously, traverse all the data and control information parameters received by Socket Flag, and store them in the buffer at the specified offset position, so as to analyze the received instructions. Each instruction is distinguished by interpreting as a command character (IAC, Interpreted As Command) mark; according to different command codes, it is converted into a corresponding string according to the data format of the Telnet protocol. After receiving the data, send the return data according to the obtained different control information parameters. According to the implemented Telnet protocol class, the return data and received data are obtained, and then the corresponding function modules are realized by the command parameters of the network equipment.

GSDF查找表是本算法中最为核心的一步,它的作用就是使亮度值与JND值之间表现出近似的线性关系。The GSDF lookup table is the most core step in this algorithm, and its function is to make an approximate linear relationship between the brightness value and the JND value.

因此,需要通过实验的方法来获得显示器的亮度曲线,还有与之对应的JND值的范围。然后设定定义范围内DDL的最小值为DDLmin,最大值为DDLmax,查找表的位宽为Pbits,测量得到的环境光亮度值为LumAmbTherefore, it is necessary to obtain the luminance curve of the display and the range of the corresponding JND value through an experimental method. Then set the minimum value of DDL within the defined range as DDL min , the maximum value as DDL max , the bit width of the lookup table as Pbits, and the measured value of ambient light brightness as Lum Amb .

快速查找表算法流程:Fast lookup table algorithm flow:

(2-1)通过实验的方法获得显示器的亮度曲线,导出DDL与亮度关系表,如图3中右下图曲线所示意;(2-1) Obtain the luminance curve of the display by an experimental method, and derive the relationship table between DDL and luminance, as shown in the lower right curve in Figure 3;

(2-2)利用医学图像和相关信息的国际标准,即医学数字成像和通信标准中第十四部分内的公式进行反向推导计算,可以得出1023个JND值所一一对应的亮度值,导出JND与亮度关系表(如图3中左下角曲线所示意);(2-2) Use the international standard for medical images and related information, that is, the formula in the fourteenth part of the medical digital imaging and communication standard to perform reverse derivation and calculation, and the brightness value corresponding to 1023 JND values can be obtained one-to-one , export the JND and brightness relationship table (as shown in the lower left corner curve in Figure 3);

(2-3)依据JND与亮度关系表,求出每个JND值一一对应三次样条插值基值,同时导出JND对样条基值表JND-Spline,当JNDindex=n时,其中,JNDindex表示JND的序列值变量,n是JND的具体数值;设立中间变量{U(n),1≤n≤1022},其中,U(n)是构造的中间变量且U(0)=0,U(1023)=0,S(1023)=0,Y(0)=0,其中,Y(0)是中间变量Y(n)在时的数值,建立调制闭值S(n)数列递推公式,此处公式中的JND(n)代表最小可觉差JND的序列值n:(2-3) According to the JND and luminance relationship table, find out that each JND value corresponds to the cubic spline interpolation base value one by one, and simultaneously derive the JND pair spline base value table JND-Spline, when JND index =n, wherein, JND index represents the serial value variable of JND, n is the specific value of JND; set up an intermediate variable {U(n), 1≤n≤1022}, where U(n) is a constructed intermediate variable and U(0)=0 , U(1023)=0, S(1023)=0, Y(0)=0, where Y(0) is the intermediate variable Y(n) in When the numerical value is established, the modulation closed value S(n) sequence recursion formula is established, where JND(n) in the formula represents the sequence value n of the minimum perceivable difference JND:

QQ (( nno )) == JJ NN DD. (( nno )) -- JJ NN DD. (( nno -- 11 )) JJ NN DD. (( nno ++ 11 )) -- JJ NN DD. (( nno -- 11 )) -- -- -- (( 44 ))

P(n)=Q(n)*Y(n-1)+2 (5)P(n)=Q(n)*Y(n-1)+2 (5)

YY (( nno )) == QQ (( nno )) -- 11 QQ (( nno )) -- -- -- (( 66 ))

KK (( nno )) == LL uu mm (( JJ NN DD. (( nno ++ 11 )) )) -- LL uu mm (( JJ NN DD. (( nno )) )) JJ NN DD. (( nno ++ 11 )) -- JJ NN DD. (( nno )) -- LL uu mm (( JJ NN DD. (( nno )) )) -- LL uu mm (( JJ NN DD. (( nno -- 11 )) )) JJ NN DD. (( nno )) -- JJ NN DD. (( nno -- 11 )) -- -- -- (( 77 ))

Uu (( nno )) == 66 ** KK (( nno )) JJ NN DD. (( nno ++ 11 )) -- JJ NN DD. (( nno -- 11 )) -- QQ (( nno )) ** Uu (( nno -- 11 )) PP (( nno )) -- -- -- (( 88 ))

S(n)=K(n)*S(n-1)+U(n) (9)S(n)=K(n)*S(n-1)+U(n) (9)

其中Q(n),P(n),Y(n),K(n)是为了得到调制闭值S(n)数列所构造的中间变量。Among them, Q(n), P(n), Y(n), and K(n) are intermediate variables constructed in order to obtain the modulated closed-value S(n) sequence.

(2-4)根据DDL与亮度关系表查找DDLmin、DDLmax对应的亮度值LumDDLmin、LumDDLmax,进而由JND与亮度关系表获得LumAmb+LumDDLmin对应的最小可用JND值JNDmin,以及LumAmb+LumDDLmax对应的最大可用JND值JNDmax,在此基础上获得显示对比度系数(DisplayContrast),计算公式如下:(2-4) Look up the luminance values Lum DDLmin and Lum DDLmax corresponding to DDL min and DDL max according to the DDL and luminance relationship table, and then obtain the minimum available JND value JND min corresponding to Lum Amb + Lum DDLmin from the JND and luminance relationship table, and The maximum available JND value JND max corresponding to Lum Amb +Lum DDLmax , on this basis, the display contrast coefficient (DisplayContrast) is obtained, and the calculation formula is as follows:

DD. ii sthe s pp aa ythe y __ cc oo nno tt rr aa sthe s tt == JNDJND mm aa xx -- JNDJND mm ii nno DDLDDL mm aa xx -- DDLDDL mm ii nno -- -- -- (( 1010 ))

(2-5)将JNDmin到JNDmax的Bartern模型曲线分成JNDmax至JNDmin区间,再利用JND-Spline得到每个区间内三次插值多项式,然后在[JNDmax,JNDmin]内做2Pbits次等间隔取样,运用假设检验中常用的P值方法生成假定值P值对JND值的查找表P-JND;对于第i个P值所对应的JND序列值(2-5) Divide the Bartern model curve from JND min to JND max into intervals from JND max to JND min , then use JND-Spline to obtain cubic interpolation polynomials in each interval, and then do 2 Pbits within [JND max , JND min ] Sampling at sub-equal intervals, using the P value method commonly used in hypothesis testing to generate the lookup table P-JND of the assumed value P value to JND value; for the JND sequence value corresponding to the i-th P value

JJ (( ii )) == JNDJND mm aa xx -- JNDJND mm ii nno 22 PP bb ii tt sthe s ** (( ii -- 11 )) ++ JNDJND mm ii nno ,, 11 ≤≤ ii ≤≤ 22 PP bb ii tt sthe s -- -- -- (( 1111 ))

只要确定了区间n,亮度值就可以通过J(i)所在JND值的区间n来求得,也即它的对应亮度值可以通过公式(12)计算出:As long as the interval n is determined, the luminance value can be obtained through the interval n of the JND value where J(i) is located, that is, its corresponding luminance value can be calculated by formula (12):

LL uu mm (( JJ (( ii )) )) == (( JJ NN DD. (( nno )) -- JJ (( ii )) )) ·&Center Dot; LL uu mm (( JJ NN DD. (( nno )) )) ++ (( JJ NN DD. (( ii )) -- JJ NN DD. (( nno -- 11 )) )) ·· LL uu mm (( JJ NN DD. (( nno -- 11 )) )) ++ SS (( nno )) 66 ·&Center Dot; [[ (( JJ NN DD. (( nno )) -- JJ (( ii )) )) 33 -- JJ NN DD. (( nno )) ++ JJ (( ii )) ]] ++ SS (( nno -- 11 )) 66 ·&Center Dot; [[ (( JJ (( ii )) -- JJ NN DD. (( nno -- 11 )) )) 33 ++ JJ NN DD. (( nno -- 11 )) -- JJ (( ii )) ]] -- -- -- (( 1212 ))

(2-6)根据JND与亮度关系表和DDL与亮度关系表找到每一个变量对应的DDL值,这个DDL值激发产生的亮度最近似于按照GSDF规定的亮度,从而获得据JND值和亮度查找DDL值的查找表。(2-6) Find the DDL value corresponding to each variable according to the relationship table between JND and brightness and the relationship table between DDL and brightness. The brightness generated by the excitation of this DDL value is the closest to the brightness specified by GSDF, so as to obtain the JND value and brightness search A lookup table of DDL values.

在远程连接完成之后,通过具体分析标准灰阶显示函数(Grayscale StandardDisplay Function,GSDF)与Barten视觉模型之间的关系,完成医用显示器的标准化显示。Barten视觉模型是一种研究人眼视觉特性的模型,根据Barten视觉模型,在人眼可识别的亮度范围内,人眼的对比敏感度是非线性的,人眼在明处的对比敏感度比在暗处高,即人眼在亮处能分辨出相对更小的亮度变化。已有的研究和分析表明,在可识别的亮度范围内人眼可分辨出1023个亮度级别JND(最小可觉差);我们用符号j代表人眼可识别的JND级别,因此可采用公式(1)来计算调制闭值S(j)。After the remote connection is completed, the standardized display of the medical display is completed by analyzing the relationship between the standard gray scale display function (Grayscale Standard Display Function, GSDF) and the Barten visual model. The Barten visual model is a model for studying the visual characteristics of the human eye. According to the Barten visual model, within the brightness range that the human eye can recognize, the contrast sensitivity of the human eye is nonlinear, and the contrast sensitivity of the human eye in bright places is higher than that in dark places. The height is high, that is, the human eye can distinguish relatively smaller brightness changes in bright places. Existing research and analysis have shown that the human eye can distinguish 1023 brightness levels JND (minimum noticeable difference) within the recognizable brightness range; we use the symbol j to represent the JND level recognizable by the human eye, so the formula ( 1) to calculate the modulation closing value S(j).

SS (( jj )) == qq 11 ·&Center Dot; Mm oo pp tt (( LL uu mm )) qq 22 dd 22 LL uu mm ++ qq 33 -- -- -- (( 11 ))

式中q1=0.1183034375,q2=3.962774805×10-5,q3=l.356243499×10-7是调制系数;d为瞳孔直径,单位:毫米(mm);Lum为目标光照亮度,单位:坎德拉/平方米(cd/m2),为光源在同一方向的光强与发光面在该方向上投影表面积之比;d的计算公式为d=4.6-2.8·tanh[0.4·log10(0.625·Lum)]。是高斯点扩展函数,是高斯点扩展函数的参数,在45°时,σ0=0.0133deg,Csph=0.0001deg/mm3,deg为角度单位。其中,μ为空间频率,单位:周/度(c/deg),Csph为球面差,σ0为σ的参照值。In the formula, q 1 =0.1183034375, q 2 =3.962774805×10 -5 , q 3 =l.356243499×10 -7 is the modulation coefficient; d is the pupil diameter, unit: millimeter (mm); Lum is the target illumination brightness, unit: Candela/square meter (cd/m 2 ), is the ratio of the light intensity of the light source in the same direction to the projected surface area of the light-emitting surface in this direction; the calculation formula of d is d=4.6-2.8·tanh[0.4·log 10 (0.625 Lum)]. is the Gaussian point spread function, is the parameter of the Gaussian point spread function, at 45°, σ 0 =0.0133deg, C sph =0.0001deg/mm 3 , and deg is the angle unit. Among them, μ is the spatial frequency, unit: cycle/degree (c/deg), C sph is the spherical aberration, and σ 0 is the reference value of σ.

人眼可识别的每一级亮度均值都是通过上一级的亮度均值计算而来,计算关系如公式(2)所示:The average brightness of each level recognizable by the human eye is calculated from the average brightness of the previous level, and the calculation relationship is shown in formula (2):

LL uu mm (( jj ++ 11 )) == LL uu mm (( jj )) 11 ++ SS (( jj )) 11 -- SS (( jj )) -- -- -- (( 22 ))

相邻两级调制阈值也就是由人眼视觉关系引入的最小可觉差概念(JND:Just-Noticeable-Difference),最小可觉差与亮度间的关系图如图2所示。The adjacent two-level modulation threshold is the concept of just-noticeable-difference (JND: Just-Noticeable-Difference) introduced by the human visual relationship. The relationship between the minimum noticeable difference and brightness is shown in Figure 2.

如图3所示,医用显示器标准灰阶校正流程主要包括以下过程:As shown in Figure 3, the standard grayscale calibration process for medical displays mainly includes the following processes:

(3-1)使用亮度计量设备采集得到不同DDL值所对应的亮度值,然后获得改变设备的特征曲线;(3-1) Use brightness metering equipment to collect brightness values corresponding to different DDL values, and then obtain the characteristic curve of changing equipment;

(3-2)GSDF曲线被定义为12bit亮度级别数学插值,每个级别的亮度值都能够根据Barten视觉模型得到。通过Barten视觉模型中的关系计算公式,可以计算出不同的JND级别值所对应的不同亮度值,其所对应亮度如下:(3-2) The GSDF curve is defined as 12bit brightness level mathematical interpolation, and the brightness value of each level can be obtained according to the Barten visual model. Through the relationship calculation formula in the Barten visual model, different brightness values corresponding to different JND level values can be calculated, and the corresponding brightness values are as follows:

loglog 1010 LL uu mm (( jj )) == aa ++ cc ·· LL nno (( jj )) ++ ee ·· [[ LL nno (( jj )) ]] 22 ++ gg ·&Center Dot; [[ LL nno (( jj )) ]] 33 ++ mm ·· [[ LL nno (( jj )) ]] 44 11 ++ bb ·&Center Dot; LL nno (( jj )) ++ dd ·&Center Dot; [[ LL nno (( jj )) ]] 22 ++ ff ·· [[ LL nno (( jj )) ]] 33 ++ hh ·· [[ LL nno (( jj )) ]] 44 ++ kk ·&Center Dot; [[ LL nno (( jj )) ]] 55 -- -- -- (( 33 ))

上式中Ln表示取自然对数;j表示人眼可识别的最小可觉差JND级别,取值范围为从1至1023的自然数;a=-1.3011877,b=-2.584019lE-2,c=8.0242636E-2,d=-l.0320229E-l,e=1.3646699E-l,f=2.8745620E-2,g=-2.5468404E-2,h=-3.1978977E-3,k=1.2992634E-4,m=1.3635334E-3,均为公式计算所需的常数系数。In the above formula, Ln represents the natural logarithm; j represents the JND level of the minimum perceivable difference that can be recognized by the human eye, and the value range is a natural number from 1 to 1023; a=-1.3011877, b=-2.584019lE-2, c= 8.0242636E-2, d=-l.0320229E-l, e=1.3646699E-l, f=2.8745620E-2, g=-2.5468404E-2, h=-3.1978977E-3, k=1.2992634E-4, m =1.3635334E-3, both are constant coefficients required for formula calculation.

(3-3)由以上各个对应关系,能够导出JND值和DDL之间的关系,然后以JND值与亮度值呈线性关系为目标进行简单调整;(3-3) From the above correspondences, the relationship between the JND value and DDL can be derived, and then a simple adjustment is made with the goal of a linear relationship between the JND value and the brightness value;

(3-4)通过对比调整前后的本地显示器设备的特征曲线相似度来评估调整符合度,其相似度越高,调整符合度就越好;反之则不好,需要继续调整。(3-4) Evaluate the adjustment compliance by comparing the similarity of the characteristic curve of the local display device before and after adjustment. The higher the similarity, the better the adjustment compliance; otherwise, it is not good and needs to continue to adjust.

由以上本发明给出的具体实施过程可以看出,对医用显示器进行远程快速专业校正的方法在不同的显示亮度条件下皆具有非常强的适用性,因而能够实现通过远程网络控制的方式对医用显示器进行高效、快速的本地校正,有效降低医用显示设备的维护成本。It can be seen from the specific implementation process given by the present invention above that the method of remote and fast professional calibration of medical displays has very strong applicability under different display brightness conditions, so it can be realized through remote network control. The monitor performs efficient and fast local calibration, effectively reducing the maintenance cost of medical display equipment.

以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。The basic principles and main features of the present invention and the advantages of the present invention have been shown and described above. Those skilled in the industry should understand that the present invention is not limited by the above-mentioned embodiments. What are described in the above-mentioned embodiments and the description only illustrate the principle of the present invention. Without departing from the spirit and scope of the present invention, the present invention will also have Variations and improvements are possible, which fall within the scope of the claimed invention. The protection scope of the present invention is defined by the appended claims and their equivalents.

Claims (6)

1. a kind of Medical display equipment remotely correcting method based on look-up table Fast Generation it is characterised in that include with Under several steps:
(1) utilize internet, after the encrypted process of command signal equipment vendor's server being sent by http agreement, send To hospital's home server to realize teleinstruction issue, hospital's home server realizes equipment vendor's clothes according to teleinstruction content Business device and the connection of hospital server, so that the manager being in far-end can realize far to the Medical display equipment of hospital Process control;
(2) according to the relation between barten vision mode and standard gray scale explicit function gsdf, obtain display brightness curve; And digital drive grade ddl that normalized brightness shown records, look-up table is generated by look-up table Fast Generation Lattice;
(3) by standard gray scale, the look-up table according to step (2) gained, shows that antidote carries out long-range school to problem display Just.
2. the Medical display equipment remotely correcting method based on look-up table Fast Generation according to claim 1, its It is characterised by, in step (1), described equipment vendor server is as follows with the method for attachment of hospital server:
First, setting native protocol address and host port, by tcp protocol remote equipment ip address and port numbers, creates Build socket to connect;Once establishment socket, you can and to be passed through to the server specified by connect functional link Send/sendto function sends data to far-end server, is received from service end with receive/receivefrom function simultaneously Data, server end then needs to make socket be connected with a local destination node using the specified interface of bind function binding, And the request on this interface is intercepted by listen function, when listening to the connection of user side, then call accept function complete Become the operation connecting, and create the incoming connection request of new socket word processing;
Subsequently, call begin receive function for Asynchronous Reception data, travel through all numbers that socket flag receives According to control information parameter, and store the relief area of the amount of specifying Offsets position, to parse to the instruction receiving;With It is construed to command character labelling and distinguish every instruction;According to different command codes, be converted to according to telnet protocol data form Corresponding character string;
After having received data, the different control information parameters according to obtaining send return data;According to the telnet agreement realized Class obtains return data and receiving data, then realizes corresponding functional module by the command parameter of the network equipment.
3. the Medical display equipment remotely correcting method based on look-up table Fast Generation according to claim 1, its It is characterised by, in step (2), described display brightness curve negotiating following methods obtain:
Symbol j represents human eye discernible just noticeable difference jnd rank, to calculate modulation using formula (1) and to close value s (j);
s ( j ) = q 1 · m o p t ( l u m ) q 2 d 2 l u m + q 3 - - - ( 1 )
In formula, q1、q2、q3For the index of modulation, d is pupil diameter, and lum is target light degree of illuminating, and that is, light source is in unidirectional light intensity Ratio with light-emitting area proj ected surface areas in this direction;It is Gaussian point spread function, Wherein, μ is spatial frequency, csphFor spherical aberration, σ0Reference value for σ, σ is the parameter of Gaussian point spread function;
The discernible every one-level luminance mean value of human eye is all to be calculated by the luminance mean value of upper level, calculated relationship such as formula (2) shown in:
l u m ( j + 1 ) = l u m ( j ) 1 + s ( j ) 1 - s ( j ) - - - ( 2 ) .
4. the Medical display equipment remotely correcting method based on look-up table Fast Generation according to claim 3, its It is characterised by, the computing formula of pupil diameter d is:
D=4.6-2.8 tanh [0.4 log10(0.625·lum)].
5. the Medical display equipment remotely correcting method based on look-up table Fast Generation according to claim 3, its It is characterised by, the method generating lookup table by described look-up table Fast Generation is as follows:
(2-1) pass through the display brightness curve obtaining, derive digital drive grade ddl and brightness relationship table;
(2-2) obtain 1023 jnd value institute brightness values correspondingly, derive just noticeable difference jnd and brightness relationship table;
(2-3) according to just noticeable difference jnd and brightness relationship table, obtain each just noticeable difference jnd value and correspond three samples Bar interpolation base value, derives just noticeable difference jnd to batten base value table jnd-spline simultaneously, works as jndindexDuring=n, wherein, jndindexRepresent the sequential value variable of jnd, n is the concrete numerical value of jnd;Set up intermediate variable { u (n), 1≤n≤1022 }, its In, u (n) is intermediate variable and u (0)=0, u (1023)=0, s (1023)=0, the y (0)=0 of construction, wherein, during y (0) is Between variable y (n) existWhen numerical value, set up modulation and close value s (n) ordered series of numbers recurrence formula, jnd (n) in formula represents herein The sequential value n of just noticeable difference jnd:
q ( n ) = j n d ( n ) - j n d ( n - 1 ) j n d ( n + 1 ) - j n d ( n - 1 ) - - - ( 4 )
P (n)=q (n) * y (n-1)+2 (5)
y ( n ) = q ( n ) - 1 q ( n ) - - - ( 6 )
k ( n ) = l u m ( j n d ( n + 1 ) ) - l u m ( j n d ( n ) ) j n d ( n + 1 ) - j n d ( n ) - l u m ( j n d ( n ) ) - l u m ( j n d ( n - 1 ) ) j n d ( n ) - j n d ( n - 1 ) - - - ( 7 )
u ( n ) = 6 * k ( n ) j n d ( n + 1 ) - j n d ( n - 1 ) - q ( n ) * u ( n - 1 ) p ( n ) - - - ( 8 )
S (n)=k (n) * s (n-1)+u (n) (9)
Wherein q (n), p (n), y (n), k (n) are to close, in order to obtain modulation, the intermediate variable that value s (n) ordered series of numbers is constructed;
(2-4) according to digital drive grade ddl and brightness relationship table search ddlmin、ddlmaxCorresponding brightness valueWherein, ddlminFor the minima of ddl as defined in the range of setting, ddlmaxFor setting the range of definition The maximum of interior ddl, and then obtained with brightness relationship table by just noticeable difference jndCorresponding I With just noticeable difference jnd value jndmin, andCorresponding maximum available jnd value jndmax, lumambFor surveying The environment bright angle value measuring, obtains on this basis and displays contrast coefficient, computing formula is as follows:
d i s p a y _ c o n t r a s t = jnd m a x - jnd m i n ddl m a x - ddl m i n - - - ( 10 )
(2-5) by jndminTo jndmaxBarten vision mode curve be divided into jndmaxTo jndminInterval, recycles to batten Base value table jnd-spline obtains cubic interpolation polynomial in each interval, then in [jndmax,jndmin] in do 2pbitsInferior Interval sampling, wherein, pbits is the bit wide of look-up table, generates assumed value p value pair with conventional p value method in hypothesis testing The look-up table p-jnd of jnd value;For the just noticeable difference jnd sequential value corresponding to i-th p value
j ( i ) = jnd m a x - jnd m i n 2 p b i t s * ( i - 1 ) + jnd m i n , 1 ≤ i ≤ 2 p b i t s - - - ( 11 )
Determine interval n, brightness value just can be tried to achieve by the interval n of j (i) place just noticeable difference jnd value, namely it Corresponding brightness value can be calculated by formula (12):
l u m ( j ( i ) ) = ( j n d ( n ) - j ( i ) ) · l u m ( j n d ( n ) ) + ( j n d ( i ) - j n d ( n - 1 ) ) · l u m ( j n d ( n - 1 ) ) + s ( n ) 6 · [ ( j n d ( n ) - j ( i ) ) 3 - j n d ( n ) + j ( i ) ] + s ( n - 1 ) 6 · [ ( j ( i ) - j n d ( n - 1 ) ) 3 + j n d ( n - 1 ) - j ( i ) ] - - - ( 12 )
(2-6) the corresponding ddl of each p value is found with brightness relationship table and ddl with brightness relationship table according to just noticeable difference jnd Value, thus obtain the look-up table searching ddl value according to just noticeable difference jnd value and brightness.
6. the Medical display equipment remotely correcting method based on look-up table Fast Generation according to claim 1, its Be characterised by, in step (3), by standard gray scale show antidote problem display is remotely corrected specifically include with Under several steps:
(3-1) collect the corresponding brightness value of different ddl values using brightness measuring equipment, then obtain the doctor needing to correct The characteristic curve of institute's local display equipment, that is to say the curve that device brightness changes with ddl and changes;
(3-2) standard gray scale explicit function gsdf curve is defined as 12bit gray scale mathematical interpolation, the brightness of each rank Value can obtain according to barten vision mode;By the relation computing formula in barten vision mode, can calculate Different brightness values corresponding to different jnd class value, brightness corresponding to it is as follows:
log 10 l u m ( j ) = a + c · l n ( j ) + e · [ l n ( j ) ] 2 + g · [ l n ( j ) ] 3 + m · [ l n ( j ) ] 4 1 + b · l n ( j ) + d · [ l n ( j ) ] 2 + f · [ l n ( j ) ] 3 + h · [ l n ( j ) ] 4 + k · [ l n ( j ) ] 5 - - - ( 3 )
In above formula, ln represents and takes natural logrithm;J represents human eye discernible just noticeable difference jnd rank, and span is from 1 To 1023 natural number;A, b, c, d, e, f, g, h, k, m are formula and calculate required constant;
(3-3) by each corresponding relation above, the relation between jnd value and ddl can be derived, then with jnd value and brightness value Linear mutually compensated for for target;
(3-4) similarity of the local display equipment characteristic curve before and after being adjusted by contrast assesses adjustment goodness of fit, its Similarity is higher, and adjustment goodness of fit is better;Otherwise then bad, need to continue adjustment.
CN201610828258.XA 2016-09-18 2016-09-18 A remote calibration method for medical display equipment based on a fast generation algorithm of look-up table Expired - Fee Related CN106357766B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610828258.XA CN106357766B (en) 2016-09-18 2016-09-18 A remote calibration method for medical display equipment based on a fast generation algorithm of look-up table

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610828258.XA CN106357766B (en) 2016-09-18 2016-09-18 A remote calibration method for medical display equipment based on a fast generation algorithm of look-up table

Publications (2)

Publication Number Publication Date
CN106357766A true CN106357766A (en) 2017-01-25
CN106357766B CN106357766B (en) 2019-04-30

Family

ID=57858043

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610828258.XA Expired - Fee Related CN106357766B (en) 2016-09-18 2016-09-18 A remote calibration method for medical display equipment based on a fast generation algorithm of look-up table

Country Status (1)

Country Link
CN (1) CN106357766B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108717847A (en) * 2018-05-18 2018-10-30 青岛海信医疗设备股份有限公司 The method of DICOM calibrations, medical display device and computer storage media
CN114417246A (en) * 2022-02-28 2022-04-29 昇显微电子(苏州)有限公司 Method and device for automatically evaluating Demura effect
WO2024044957A1 (en) * 2022-08-30 2024-03-07 京东方科技集团股份有限公司 Display apparatus and brightness adjustment method therefor

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1845232A (en) * 2005-04-08 2006-10-11 浙江大学 Liquid crystal display GAMMA curve and color temperature automatic measurement and correction method and system
CN101620835A (en) * 2009-03-27 2010-01-06 青岛海信电器股份有限公司 GAMMA correction method of display equipment and device therefor
CN101697230A (en) * 2009-11-03 2010-04-21 蒋慧琴 Device and method for adaptive enhancement treatment of medical image
CN103680438A (en) * 2013-11-22 2014-03-26 武汉精立电子技术有限公司 Synchronous correction method for gamma curve and flicker of LCD
CN104183229A (en) * 2013-07-09 2014-12-03 上海联影医疗科技有限公司 Display device correcting method and correcting method
CN105489193A (en) * 2016-01-29 2016-04-13 明基电通有限公司 Image adjusting method and display device

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1845232A (en) * 2005-04-08 2006-10-11 浙江大学 Liquid crystal display GAMMA curve and color temperature automatic measurement and correction method and system
CN101620835A (en) * 2009-03-27 2010-01-06 青岛海信电器股份有限公司 GAMMA correction method of display equipment and device therefor
CN101697230A (en) * 2009-11-03 2010-04-21 蒋慧琴 Device and method for adaptive enhancement treatment of medical image
CN104183229A (en) * 2013-07-09 2014-12-03 上海联影医疗科技有限公司 Display device correcting method and correcting method
CN103680438A (en) * 2013-11-22 2014-03-26 武汉精立电子技术有限公司 Synchronous correction method for gamma curve and flicker of LCD
CN105489193A (en) * 2016-01-29 2016-04-13 明基电通有限公司 Image adjusting method and display device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孙安玉,江贵平,吕庆文: "DICOM灰阶图像软拷贝快速标准化显示方法", 《医学物理与工程学》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108717847A (en) * 2018-05-18 2018-10-30 青岛海信医疗设备股份有限公司 The method of DICOM calibrations, medical display device and computer storage media
CN108717847B (en) * 2018-05-18 2020-11-20 青岛海信医疗设备股份有限公司 DICOM calibration method, medical display device and computer storage medium
CN114417246A (en) * 2022-02-28 2022-04-29 昇显微电子(苏州)有限公司 Method and device for automatically evaluating Demura effect
CN114417246B (en) * 2022-02-28 2024-12-06 昇显微电子(苏州)股份有限公司 A method and device for automatically evaluating Demura effect
WO2024044957A1 (en) * 2022-08-30 2024-03-07 京东方科技集团股份有限公司 Display apparatus and brightness adjustment method therefor

Also Published As

Publication number Publication date
CN106357766B (en) 2019-04-30

Similar Documents

Publication Publication Date Title
CN106357766A (en) Medical display device remote correction method based on lookup table rapid generation algorithm
CN105740856A (en) Method for reading readings of pointer instrument based on machine vision
CN106053024B (en) A kind of LED light source preference degree prediction technique towards monochromatic system object
CN111896549A (en) A system and method for monitoring building cracks based on machine learning
CN115035364A (en) A pointer meter reading method based on deep neural network
CN209132626U (en) A kind of wisdom classroom management system based on Internet of Things
CN111783664A (en) Folding line compensation method and device for folding screen, electronic equipment and storage medium
CN113834585B (en) Method, system, device, medium and electronic equipment for detecting silicon rod in reduction furnace
CN108535684A (en) The method and system that the wireless transmitter spatially separated is positioned
CN102547372B (en) Method for determining color lookup table parameters of liquid crystal display television and system device therefore
CN111198225A (en) Intelligent pesticide residue rapid detection instrument based on block chain traceability system
CN104008527B (en) Method for defogging single image
Ullah et al. Gray‐Level Image Transformation of Paved Road Cracks with Metaphorical and Computational Analysis
CN103793868A (en) Data transmission method, device and system
CN112102252B (en) Method and device for detecting appearance defects of welding spots of microstrip antenna
Liu et al. Local dimming algorithm of automotive LCD instrument based on otsu and maximum entropy
CN114332106A (en) Pointer meter dial plate unfolding method based on approximate pole center
KR20180082172A (en) Method and system for skin analysis consifering illumination characteristic
WO2020034249A1 (en) Cloud data transmission apparatus and method based on crude oil nuclear magnetic resonance analyzer
CN113095293A (en) Factory calibration method for Internet of things water meter based on image recognition
CN112001738A (en) Method for constructing multi-factor logistics product price model and application method thereof
CN119005117A (en) Design method of visible light emission equalizer based on LED equivalent circuit
CN115361290B (en) Configuration comparison method, device, electronic equipment and storage medium
CN119233036B (en) Remote Ultrasound Image Transmission System
CN104793494B (en) A kind of automatic installation and locating method of fatigue tester test specimen based on machine vision

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190430

CF01 Termination of patent right due to non-payment of annual fee