颜色模型转换公式
2010-07-01 20:08
颜色模型转换公式
为了用计算机来表示和处理颜色，必须采用定量的方法来描述颜色，即建立颜色模型。目前广泛使用的颜色模型觉颜色模型。计算颜色模型又称为色度学颜色模型，主要应用于纯理论研究和计算推导；工业颜色模型侧重于实直接接口的颜色模型和控制。
1、计算颜色模型有CIE的RGB、XYZ、Luv、LCH、LAB、UCS、UVW。
2、工业颜色模型NTSC的RGB、YUV、YIQ、CMYK、YCbCr。
3、视觉颜色模型有HS*系列，包含HSL、HSV(B)。
常见Color Formula/Matrix如下：
1、RGB<->CIE XYZ
a、RGB<->CIE XYZ REC601
|X| | 0.607 0.174 0.201| |R|
|Y| = | 0.299 0.587 0.114| * |G|
|Z| | 0.000 0.066 1.117| |B|
|R| | 1.910 -0.532 -0.288| |X|
|G| = |-0.985 1.999 -0.028| * |Y|
|B| | 0.058 -0.118 0.898| |Z|
b、RGB<->CIE XYZ REC709
|X| | 0.412 0.358 0.180| |R|
|Y| = | 0.213 0.715 0.072| * |G|
|Z| | 0.019 0.119 0.950| |B|
|R| | 3.241 -1.537 -0.499| |X|
|G| = |-0.969 1.876 -0.042| * |Y|
|B| | 0.056 -0.204 1.057| |Z|
c、RGB<->CIE XYZ ITU
|X| = | 0.431 0.342 0.178| |R|
|Y| = | 0.222 0.707 0.071| * |G|
|Z| = | 0.020 0.130 0.939| |B|
|R| | 3.063 -1.393 -0.476| |X|
|G| = |-0.969 1.876 0.042| * |Y|
|B| | 0.068 -0.229 1.069| |Z|
2、RGB<->CMYK
a、RGB -> CMYK
K = min(1-R, 1-G, 1-B)
C = (1-R-K) / (1-K)
M = (1-G-K) / (1-K)
Y = (1-B-K) / (1-K)
b、CMYK->RGB
R = 1 - min(1, C *(1-K) + K)
G = 1 - min(1, M *(1-K) + K)
B = 1 - min(1, Y *(1-K) + K)
3、RGB->CIE RGB
|Rcie| | 1.167 -0.146 -0.151| |R|
|Gcie| = | 0.144 0.753 0.159| * |G|
|Bcie| |-0.001 0.059 1.128| |B|
4、CIE XYZ->CIE Lab
D65 xn = 0.312713 , yn = 0.329016
L = 116 * ( ( Y / Yn ) ^ ( 1 / 3 ) ) Y/Yn > 0.008856
L = 903.3 * Y / Yn Y/Yn <= 0.008856
a = 500*(f(X/Xn)-f(Y/Yn))
b = 200*(f(Y/Yn)-f(Z/Zn))
5、RGB->HSV
V=max(R,G,B)
S=(V-min(R,G,B))*255/V if V!=0, 0 otherwise
(G - B)*60/S, if V=R
H= 180+(B - R)*60/S, if V=G
240+(R - G)*60/S, if V=B
若H<0，则H=H+360
使用上面从0°到360°变化的公式计算色调（hue）值，确保它们被2 除后能试用于8位。
6、RGB<->YIQ
YIQ是北美电视系统（NTSC）所采用的。Y指亮度，I和Q指色调，描述图像的色彩和饱和度。
|Y| | 0.299 0.587 0.114| |R|
|I| ＝| 0.596 -0.274 -0.322| * |G|
|Q| | 0.211 -0.523 0.312| |B|
|R| | 1 0.956 0.621| |Y|
|G| = | 1 -0.272 -0.647| * |I|
|B| |-1 -1.106 -1.703| |Q|
7、RGB<->YUV
YUV是欧洲的电视系统所采用（PAL）。Y指亮度，U和V指色调。
|Y| | 0.299 0.587 0.114| |R|
|U| = |-0.148 -0.289 0.437| * |G|
|V| | 0.615 -0.515 -0.100| |B|
|R| | 1 0 1.140| |Y|
|G| = | 1 -0.395 -0.581| * |I|
|B| |-1 2.032 0 | |Q|
8、RGB<->YCbCr
JPEG采用的颜色模型是YCbCr。它是从YUV颜色模型衍生来的。其中Y指亮度，而Cb和Cr是将U和V
|Y | | 0.2990 0.5870 0.1140 0 | |R|
|Cb| |-0.1687 -0.3313 0.5000 128| * |G|
|Cr | = | 0.5000 -0.4187 -0.0813 128| |B|
|1 | | 0 0 0 1 | |1|
|R| | 1 1.40200 0 | |Y |
|G| = | 1 -0.34414 -0.71414| * |Cb-128|
|B| | 1 1.77200 0 | |Cr-128 |
标准参考代码：
https://www.doczj.com/doc/cd6159666.html,/math.php?MATH=M22#text22
出处：https://www.doczj.com/doc/cd6159666.html,/jonfei/archive/2007/04/25/1583885.aspx
CMYK_NORMAL_COLORSPACE:
r = ((maxval-k)*(maxval-c))/maxval;
g = ((maxval-k)*(maxval-m))/maxval;
b = ((maxval-k)*(maxval-y))/maxval;
CMYK_ADOBE_COLORSPACE:
r = (k*c)/maxval;
g = (k*m)/maxval;
b = (k*y)/maxval;
This differs from the formula used in groff (as submitted by Gaius; to
be available soon from the CVS):
r = maxval - min(maxval, c*(maxval-k)/maxval + k);
g = maxval - min(maxval, m*(maxval-k)/maxval + k);
b = maxval - min(maxval, y*(maxval-k)/maxval + k);
Lab 颜色空间
Lab颜色空间是由CIE(国际照明委员会)制定的一种色彩模式。自然界中任何一点色都可以在Lab空间中表达出这种模式是以数字化方式来描述人的视觉感应，与设备无关，所以它弥补了RGB和CMYK模式必须依赖于设备RGB模式和CMYK模式的色彩空间大。这就意味着RGB以及CMYK所能描述的色彩信息在Lab空间中都能亮度；a的正数代表红色，负端代表绿色；b的正数代表黄色，负端代表兰色(a,b)有L=116f(y)-16, a=500[f(x/0 f(x)=7.787x+0.138, x<0.008856; f(x)=(x)1/3,x>0.008856
第六讲图像类型与 彩色模型的转换 【目录】 一、图像类型的转换 (1) 1、真彩图像→索引图像 (3) 2、索引图像→真彩图像 (3) 3、真彩图像→灰度图像 (4) 4、真彩图像→二值图像 (4) 5、索引图像→灰度图像 (5) 6、灰度图像→索引图像 (6) 7、灰度图像→二值图像 (7) 8、索引图像→二值图像 (8) 9、数据矩阵→灰度图像 (9) 二、彩色模型的转换 (9) 1、图像的彩色模型 (10) 2、彩色转换函数 (10) 三、纹理映射 (13) 【正文】 一、图像类型的转换
1、真彩图像→索引图像 【格式】X =d i t h e r (R G B ,m a p ) 【说明】按指定的颜色表m a p 通过颜色抖动实现转换 【输入】R G B 可以是d o u b l e 或u i n t 8类型 【输出】X 超过256色则为d o u b l e 类型，否则输出为u i n t 8型 【例】 C L F ,R G B =i m r e a d ('f l o w e r s .t i f '); 100 200 300 400 500 50100150200250300350 100 200 300 400 500 50100150200250300350 【输出】R G B 为d o u b l e 类型 【例】 C L F ,l o a d t r e e s ; R G B =i n d 2r g b (X ,m a p ); s u b p l o t (1,2,1);s u b i m a g e (X ,m a p );t i t l e ('索引图') s u b p l o t (1,2,2);s u b i m a g e (R G B );t i t l e ('真彩图')
§2.3.1 坐标系的分类 正如前面所提及的,所谓坐标系指的是描述空间位置的表达形式,即采用什么方法来表示空间位置。人们为了描述空间位置，采用了多种方法，从而也产生了不同的坐标系，如直角坐标系、极坐标系等。 在测量中常用的坐标系有以下几种： 一、空间直角坐标系 空间直角坐标系的坐标系原点位于参考椭球的中心，Z 轴指向参考椭球的北极，X 轴指向起始子午面与赤道的交点，Y 轴位于赤道面上且按右手系与X 轴呈90°夹角。某点在空间中的坐标可用该点在此坐标系的各个坐标轴上的投影来表示。空间直角坐标系可用图2-3来表示： 图2-3 空间直角坐标系 二、空间大地坐标系 空间大地坐标系是采用大地经、纬度和大地高来描述空间位置的。纬度是空间的点与参考椭球面的法线与赤道面的夹角；经度是空间中的点与参考椭球的自转轴所在的面与参考椭球的起始子午面的夹角；大地高是空间点沿参考椭球的法线方向到参考椭球面的距离。空间大地坐标系可用图2-4来表示：
图2-4空间大地坐标系 三、平面直角坐标系 平面直角坐标系是利用投影变换，将空间坐标空间直角坐标或空间大地坐标通过某种数学变换映射到平面上，这种变换又称为投影变换。投影变换的方法有很多，如横轴墨卡托投影、UTM 投影、兰勃特投影等。在我国采用的是高斯-克吕格投影也称为高斯投影。UTM 投影和高斯投影都是横轴墨卡托投影的特例，只是投影的个别参数不同而已。 高斯投影是一种横轴、椭圆柱面、等角投影。从几何意义上讲，是一种横轴椭圆柱正切投影。如图左侧所示，设想有一个椭圆柱面横套在椭球外面，并与某一子午线相切（此子午线称为中央子午线或轴子午线），椭球轴的中心轴CC ’通过椭球中心而与地轴垂直。 高斯投影满足以下两个条件： 1、 它是正形投影； 2、 中央子午线投影后应为x 轴，且长度保持不变。 将中央子午线东西各一定经差（一般为6度或3度）范围内的地区投影到椭圆柱面上，再将此柱面沿某一棱线展开，便构成了高斯平面直角坐标系，如下图２－５右侧所示。 图2-5 高斯投影 x 方向指北，y 方向指东。 可见，高斯投影存在长度变形，为使其在测图和用图时影响很小，应相隔一定的地区，另立中央子午线，采取分带投影的办法。我国国家测量规定采用六度带和三度带两种分带方法。六度带和三度带与中央子午线存在如下关系： 366 N L ＝中； n L 33＝中 其中，N 、n 分别为6度带和3度带的带号。
空间直角坐标系与大地坐标系转换程序 #include
Excel中日期格式转换问题 1.数字格式转年月日格式 工作表中出生年月日格式格式是如：198707141987090620090902如果想改成1987-7-14或改成1987年7月14日。 1、就用此函数插入一列如B列，B1=(LEFT(A1,4)&"-"&MID(A1,5,2)&"-"&RIGHT(A1,2));其它单元格可以用此公式自动填充，然后把此单元格设置成日期1987-7-14； 2、插入一列如B列，B1=(LEFT(A1,4)&"年"&MID(A1,5,2)&"月"&RIGHT(A1,2)&"日");其它单元格可以用此公式自动填充，然后把此单元格设置成日期1987年7月14日。 把这种格式2009.9.2转换成2009-09-02这样的格式：假如原始数据在A列从A1开始，则在B1输入公式=TEXT(SUBSTITUTE(A1,".","-"),"yyyy-mm-dd")，然后下拉复制公式即可。 2.在Excel中将非日期格式数据转换为日期格式 如果Excel表格中有大量的类似“19841006”、“19841029”非日期格式的数据要转换为日期格式，如“1984-10-6”，可以用分列的方法： 1）.选择需要转换单元格或区域，单击菜单“数据→分列”，弹出“文本分列向导-3之步骤1”对话框。如果是Excel2007，则在功能区中选择“数据”选项卡，在“数据工具”组中，单击“分列”按钮。
2）.单击两次“下一步”，在“文本分列向导-3之步骤3”中，在“列数据格式”中选择“日期”。还可以根据需要，在其后的下拉列表中选择一种日期格式。本例为默认的“YMD”格式。 3）.单击“完成”按钮，Excel会在原单元格或区域中将数值格式的“19841006”和文本格式的“1984.10.6”转换为日期格式“1984-10-6”。 另外，对于“1984.10.6”文本格式的数据，还可以将“.”替换为“-”，也可以转换为日期格式。 3.将日期格式转换为数字格式 在日期后插入一列取函数=text(a1,"yyyymmddhhmm")yyyy表示年mm表示月dd表示日hh表示时mm表示分a1表示日期位置（第一行第一列），然后向下拖曳填充即可改变其他各行日期 4.Excel文本日期转换为日期格式 什么是序列号? Excel可将日期存储为可用于计算的序列号。默认情况下，1900年1月1日的序列号为1，2008年1月1日的序列号为39,448，这是因为它距1900年1月1日有39,448天。 2、若要将转换公式复制到相邻单元格的区域中，请选择在其中键入公式的单元格，然后拖动填充柄，使之覆盖大小与包含文本日期的单元格区域匹配的空单元格区域。 拖动填充柄之后，应有序列号与包含文本日期的单元格区域相对应的单元格区域。
坐标转换之计算公式 一、参心大地坐标与参心空间直角坐标转换 1名词解释： A ：参心空间直角坐标系： a) 以参心0为坐标原点； b) Z 轴与参考椭球的短轴（旋转轴）相重合； c) X 轴与起始子午面和赤道的交线重合； d) Y 轴在赤道面上与X 轴垂直，构成右手直角坐标系0-XYZ ； e) 地面点P 的点位用（X ，Y ，Z ）表示； B ：参心大地坐标系： a) 以参考椭球的中心为坐标原点，椭球的短轴与参考椭球旋转轴重合； b) 大地纬度B ：以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B ； c) 大地经度L ：以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L ； d) 大地高H ：地面点沿椭球法线至椭球面的距离为大地高H ； e) 地面点的点位用（B ，L ，H ）表示。 2 参心大地坐标转换为参心空间直角坐标： ?? ???+-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中，N 为椭球面卯酉圈的曲率半径，e 为椭球的第一偏心率，a 、b 椭球的长短半 径，f 椭球扁率，W 为第一辅助系数 a b a e 2 2-= 或 f f e 1*2-= W a N B W e =-=22sin *1( 3 参心空间直角坐标转换参心大地坐标
[]N B Y X H H e N Y X H N Z B X Y L -+=+-++==cos ))1(**)()(*arctan( )arctan(2 2222 二 高斯投影及高斯直角坐标系 1、高斯投影概述 高斯-克吕格投影的条件：1. 是正形投影；2. 中央子午线不变形 高斯投影的性质：1. 投影后角度不变；2. 长度比与点位有关，与方向无关； 3. 离中央子午线越远变形越大 为控制投影后的长度变形，采用分带投影的方法。常用3度带或6度带分带，城市或工 程控制网坐标可采用不按3度带中央子午线的任意带。 2、高斯投影正算公式： 5 2224253 2236 4254 42232)5814185(cos 120 )1(cos 6 cos )5861(cos sin 720 495(cos sin 24 cos sin 2l t t t B N l t B N Bl N y l t t B B N l t B B N Bl B N X x ηηηηη-++-++-+=+-+++-++=） 3、高斯投影反算公式：
空间坐标转换说明 集团文件版本号：（M928-T898-M248-WU2669-I2896-DQ586-M1988）
坐标转换说明 GPS 接收机接收到GPS （大地坐标：经度、纬度和高度值）信号后，并不利于显示，需要将大地坐标进行转换，现选用东北天坐标系（也叫站心坐标系）作为显示的依据。 GPS 接收机接收到的第一个信号L （经度）、B （纬度）和H （高度），作为东北天坐标系的原点。当接收到第二个信号时L 1、B 1和H 1，应用坐标转换公式，转换到东北天坐标系下进行显示。依次类推，凡是接收到的GPS 信号都转换到东北天坐标系下进行显示，在东北天坐标系下预测出来的坐标值通过坐标转换公式在显示屏上显示大地坐标（经度、纬度和高度）。 1.大地坐标与直角坐标的相互转化 对空间某一点，大地坐标系(L ,B ,H )到直角坐标系(X ,Y ,Z )的转换关系如下： ?? ? ?? +-=+=+=B H e N Z L B H N Y L B H N X sin ])1([sin cos )(cos cos )(2（1） 由直角坐标系(X ,Y ,Z )转化到大地坐标系(L ,B ,H )的公式如下： ??? ? ??? --=+-++==)1(sin /]})1((/[)(arctan{) /arctan(2222e N B Z H H e N Y X H N Z B X Y L （2） 式中：B e a N 22sin 1/-=，N 为该点的卯酉圈曲率半径；2222/)(a b a e -=，a 、 b 、e 分别为该大地坐标系对应参考椭球的长半轴、短半轴和第一偏心率。长半轴 a =6378137±2m ，短半轴 b =6356.7523142km ，90130066943799.02=e 。 从公式（2）看出，经度比较容易求得，纬度和高度必须通过迭代计算获直接计算得到。迭代计算的次序为：N H B →→，通常迭代四次可以达到H 优于0.001m ，B 优于0.00001''的计算精度；教科书中给出的直接法计算公式比较繁琐，有的计算公式的应用条件受到一定限制，例如要求大地高度小于10000m 时，才能使B 、H 达到上述计算精度，有的直接计算公式精度较低。 根据[张华海]提供的方法，本文建议采用该方法将直角坐标（X ,Y ,Z ）转变成大地坐标（L ,B ,H ）。该方法的公式形式比较简便，B 、H 的计算精度高；用计算出
2011-04-22 15:10 Excel中日期格式转换问题 Excel中日期格式转换问题 1.数字格式转年月日格式 工作表中出生年月日格式格式是如：19870714 19870906 20090902如果想改成1987-7-14或改成1987年7月14日。 1、就用此函数插入一列如B列， B1=(LEFT(A1,4)&"-"&MID(A1,5,2)&"-"&RIGHT(A1,2)); 其它单元格可以用此公式自动填充，然后把此单元格设置成日期1987-7-14 ； 2、插入一列如B列，B1=(LEFT(A1,4)&"年"&MID(A1,5,2)&"月"&RIGHT(A1,2)&"日"); 其它单元格可以用此公式自动填充，然后把此单元格设置成日期1987年7月14 日。 把这种格式2009.9.2转换成2009-09-02这样的格式：假如原始数据在A 列从A1开始，则在B1输入公式=TEXT(SUBSTITUTE(A1,".","-"),"yyyy-mm-dd")，然后下拉复制公式即可。 2.在Excel中将非日期格式数据转换为日期格式 如果Excel表格中有大量的类似“19841006”、“19841029”非日期格式的数据要转换为日期格式，如“1984-10-6”，可以用分列的方法： 1）.选择需要转换单元格或区域，单击菜单“数据→分列”，弹出“文本分列向导- 3之步骤1”对话框。如果是Excel 2007，则在功能区中选择“数据”选项卡，在“数据工具”组中，单击“分列”按钮。
大地(高斯平面)坐标系工程坐标系转换大地坐标系--->工程坐标系 ======================== 待转换点为P，大地坐标为：Xp、Yp 工程坐标系原点o: 大地坐标：Xo、Yo 工程坐标：xo、yo 工程坐标系x轴之大地方位角：a dX=Xp-Xo dY=Yp-Yo P点转换后之工程坐标为xp、yp: xp=dX*COS(a)+dY*SIN(a)+xo yp=-dX*SIN(a)+dY*COS(a)+yo 工程坐标系--->大地坐标系 ======================== 待转换点为P，工程坐标为：xp、yp 工程坐标系原点o: 大地坐标：Xo、Yo 工程坐标：xo、yo 工程坐标系x轴之大地方位角：a dx=xp-xo dy=yp-yo P点转换后之工程坐标为xp、yp: xp=Xo+dx*COS(a)-dy*SIN(a)
yp=Yo+dx*SIN(a)+dy*COS(a) 坐标方位角计算程序 置镜点坐标：ZX ZY 后视点坐标：HX HY 方位角：W 两点间距离: S Lb1 0← ｛A, B, C, D｝← A〝ZX=〞:B〝ZY=〞:C〝HX=〞:D 〝HY=〞:W=tg1((D-B)÷(C-A)):(D-B)>0=>(C-A)>0=>W=W:∟∟(D-B)>0=>(C-A)<0=>W=W+180:∟∟(D-B)<0=>(C-A)<0=>W=W+180:∟∟(D-B)<0=>(C-A)>0=>W=360+W∟∟W=W◢ S=√((D-B)2+(C-A)2) ◢ Goto 0← CASIO fx－4500p坐标计算程序 根据坐标计算方位角 W＝W＋360△W：“ALF(1~2)＝”L1 A“X1＝”：B“Y1＝”：Pol(C“X2”－A，D“Y2”－B：“S＝”▲W＜0 直线段坐标计算 L1 X“X(0)”：Y“Y(0)”：S“S(0)”：A“ALF” L2 Lb1 2 L3 {L}：L“LX”
坐标转换方法 空间直角坐标系如果其原点不动，绕着某一个轴旋转而构成的新的坐标系，这个过程就叫做坐标旋转。在旧坐标系中的坐标与在旋转后新坐标系中的坐标有一定的转换关系，这种转换关系可以用转换矩阵来表示。 如图5.7，直角坐标系XYZ，P点的坐标为（x, y, z），其相应的在XY 平面，XZ平面，YZ平面分别为M（x, y,0），Q（x,0, z）和N（0, y, z）。 图5.7直角坐标系XYZ 设?表示第j 轴的旋转角度，R j (?) 表示绕第j 轴的旋转，其正方向是沿坐标轴向原点看去的逆时针方向。很明显当j 轴为旋转轴时，它对应的坐标中的j 分量是不变的。由于直角坐标系是对称的，下面我们以绕Z轴旋转为例推导其旋转变换矩阵，其它两个轴推导和它是一样的。 设图5.7的坐标绕Z轴逆时针旋转θ角度，新坐标为X 'Y'Z'，如图5.8所示： 图5.8 坐标绕Z 轴逆时针旋转θ角度 由于坐标中的z 分量不变，我们可以简化地在XY 平面进行分分析，如图
5.9所示： 图5.9坐标绕Z 轴逆时针旋转θ 角度的XY 平面示意图 点 M X 和点M X ' 分别是M 点在X 轴和X '轴的投影。如图5.9 cos cos() sin sin() X X X X x OM OM MOM OM y MM OM MOM OM ?θ?θ==∠=-??==∠=-? （5-1） cos cos sin sin X X X X x OM OM MOM OM y MM OM MOM OM ? ?'''''==∠=??'==∠=? （5-2） 把（5-1）式按照三角函数展开得： cos cos sin sin sin cos cos sin x OM OM y OM OM ?θ?θ ?θ?θ=+??=+? （5-3） 把（5-2）式代入（5-3）式得： cos sin sin cos x x y y x y θθ θθ''=+??''=-+? （5-4） 坐标中的z 分量不变，即z = z'这样整个三维坐标变换就可以写成（用新坐标表 示旧坐标） cos sin sin cos x x y y x y z z θθ θθ''=+? ?''=-+??' =? （5-5） 把式（5-5）用一个坐标旋转变换矩阵R Z (θ) 表示可以写成：
由传递函数转换成状态空间模型——方法多!!! SISO 线性定常系统 高阶微分方程化为状态空间表达式 SISO ()()()()()()m n u b u b u b y a y a y a y m m m n n n n ≥+++=++++--- 1102211 )(2 211110n n n n m m m a s a s a s b s b s b s G +++++++=--- 假设1+=m n 外部描述 ←—实现问题：有了部结构—→模拟系统 部描述 SISO ? ??+=+=du cx y bu Ax x 实现问题解决有多种方法，方法不同时结果不同。 一、 直接分解法 因为 1 0111 11()()()()()()()() 1m m m m n n n n Y s Z s Z s Y s U s Z s U s Z s b s b s b s b s a s a s a ----?=? =?++++++++ ???++++=++++=----) ()()() ()()(11 11110s Z a s a s a s s U s Z b s b s b s b s Y n n n n m m m m 对上式取拉氏反变换，则 ? ??++++=++++=----z a z a z a z u z b z b z b z b y n n n n m m m m 1) 1(1)(1)1(1)(0 按下列规律选择状态变量，即设)1(21,,,-===n n z x z x z x ，于是有
?????? ?+----===-u x a x a x a x x x x x n n n n 12113 221 写成矩阵形式 式中，1-n I 为1-n 阶单位矩阵，把这种标准型中的A 系数阵称之为友阵。只要系统状态方程的系数阵A 和输入阵b 具有上式的形式，c 阵的形式可以任意，则称之为能控标准型。 则输出方程 121110x b x b x b x b y m m n n ++++=-- 写成矩阵形式 ??????? ? ????????=--n n m m x x x x b b b b y 12101 1][ 分析c b A ,,阵的构成与传递函数系数的关系。 在需要对实际系统进行数学模型转换时，不必进行计算就可以方便地写出状态空间模型的A 、b 、c 矩阵的所有元素。 例：已知SISO 系统的传递函数如下，试求系统的能控标准型状态空间模型。 4 2383)()(2 3++++=s s s s s U s Y 解：直接得到系统进行能控标准型的转换，即
电子表格中有关日期的几个变换与计算 从身份证中提取出生日期=date(mid(A1,7,4),mid(A1,11,2),mid(A1,13,2))直接提取出日期。 身份证算年龄 如要算指定日期的年龄,则设A1文本格式的18位身份证号,B1为指定日期,在C1输入公式如下: =YEAR(TODAY())-V ALUE(MID(A1,7,4)) =YEAR(TODAY())-IF(LEN(TRIM(A1))>15,V ALUE(MID(A1,7,4)),V ALUE("1 9"&MID(A6,7,2))) 此公式15位、18位通用 到指定日期年龄 假设A1为身份号，15、18位都可，在B1显示此人年龄（至2009年底计算），公式如下： =2009-LEFT(IF(LEN(A1)=15,19&TEXT(MID(A1,7,6),"00-00-00"),TEXT(MID(A1, 7,8),"0000")),4) 到指定日期年龄 =DA TEDIF(RIGHT(TEXT(MID(A1,7,11)-500,"#-00-00,"),10),指定日期,"Y") =MOD(YEAR(指定日期)-MID(A1,7,LEN(A1)/2-5),190) 身份证号快速输入 方法一、单元格设置---------数字--------自定义------------类型
“输入41062100000000000#”（6位数加11个0加#，这样设置的单元格，只需要输入后12位就行了，前面的会自动添加） 方法二、用公式“=A1&B1做” 从身份证中提取数字串 例如，某人身份证为310123************（18位），输入到"A1"单元格，建议用把单元格改成"文本"，不然会显示为科学计数法。 一、提取数字 方法一、在"B1"单元格输入：=left(A1,14)，出现结果是“31012319010203”，然后在"C1"单元格里输入=right(B1,8)，出现结果是“19010203”，基本工程完成了。 方法二、两个命令合并使用，在"D1"单元格里输入： =right(left(A1,14),8)。 方法三、函数【mid】在"E1"单元格输入公式=mid(A1,7,8)试试，就算出了正确的结果“19010203” 二、转化成日期格式函数【date】 例如，=date(1901,02,03) 输出结果为日期格式的“1901-2-3” 三、也可以这样：按照如下操作提取出年月日：提取年，在"F1"单元格输入：=mid(A1,7,4)，输出结果为“1901”。提取月，在"G1"单元格输入：=mid(A1,11,2)，输出结果为“02”。提取日，在"H1"单元格
灰度图像处理程序代码代码 1.二值图像 function erzhi_Callback(hObject, eventdata, handles) % hObject handle to erzhi (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) axes(handles.axes2); x=(handles.img); if isrgb(x) msgbox('这是彩色图像,不能转换为二值图像','转换失败'); else j=im2bw(x); imshow(j); end 2.图像腐蚀 function fushi_Callback(hObject, eventdata, handles) % hObject handle to fushi (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) axes(handles.axes2); x=(handles.img); if isrgb(x) msgbox('这是彩色图像,不能进行图像腐蚀','失败'); else j=im2bw(x); se=eye(5); bw=bwmorph(j,'erode'); imshow(bw); 3.创建索引图像 function chuanjian_Callback(hObject, eventdata, handles) % hObject handle to chuanjian (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) axes(handles.axes2); x=(handles.img); if isrgb(x) msgbox('这是彩色图像,不能创建索引图像','创建失败'); else y=grayslice(x,16); axes(handles.axes2); imshow(y,jet(16)); end 4.轮廓图
大地坐标与直角空间坐标转换计算公式 一、参心大地坐标与参心空间直角坐标转换 1名词解释： A ：参心空间直角坐标系： a) 以参心0为坐标原点； b) Z 轴与参考椭球的短轴（旋转轴）相重合； c) X 轴与起始子午面和赤道的交线重合； d) Y 轴在赤道面上与X 轴垂直，构成右手直角坐标系0-XYZ ； e) 地面点P 的点位用（X ，Y ，Z ）表示； B ：参心大地坐标系： a) 以参考椭球的中心为坐标原点，椭球的短轴与参考椭球旋转轴重合； b) 大地纬度B ：以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B ； c) 大地经度L ：以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L ； d) 大地高H ：地面点沿椭球法线至椭球面的距离为大地高H ； e) 地面点的点位用（B ，L ，H ）表示。 2 参心大地坐标转换为参心空间直角坐标： ?? ? ?? +-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中，N 为椭球面卯酉圈的曲率半径，e 为椭球的第一偏心率，a 、b 椭球的长短半径，f 椭球扁率，W 为第一辅助系数 a b a e 2 2-= 或 f f e 1 *2-= W a N B W e = -=22 sin *1( 西安80椭球参数： 长半轴a=6378140±5（m ）
短半轴b=6356755.2882m 扁 率α=1/298.257 3 参心空间直角坐标转换参心大地坐标 [ ] N B Y X H H e N Y X H N Z B X Y L -+= +-++==cos ))1(**)() (*arctan() arctan(2 22 2 2 二 高斯投影及高斯直角坐标系 1、高斯投影概述 高斯-克吕格投影的条件：1. 是正形投影；2. 中央子午线不变形 高斯投影的性质：1. 投影后角度不变；2. 长度比与点位有关，与方向无关； 3. 离中央子午线越远变形越大 为控制投影后的长度变形，采用分带投影的方法。常用3度带或6度带分带，城市或工程控制网坐标可采用不按3度带中央子午线的任意带。 2、高斯投影正算公式： 52224253 2236 425442232)5814185(cos 120 )1(cos 6 cos )5861(cos sin 720 495(cos sin 24cos sin 2l t t t B N l t B N Bl N y l t t B B N l t B B N Bl B N X x ηηηηη-++-++-+=+-+++-++ =） 3、高斯投影反算公式：
PIM->PSM 模型转换的途径 mdaSky UML软件工程组织 由MDA 的PIM（平台独立模型）向PSM（平台特定模型）转换的方法目前尚未实现标准化。因此目前市售的工具不得不利用自主方法进行这部分的处理。由PIM 向PSM 的转换方法由于将在2004 年实现标准化，只有这个重要的步骤标准化了，才更加有利于MDA 这项技术的推广。 2004 年将是MDA 大发展的一年，为什么这样说，我们来看看业界一些重要的公司是如何应对MDA 这项技术的。最近，美国Compuware 的OptimalJ 等基于对象技术标准化团体美国OMG （Object Management Group ）倡导的模型驱动架构（MDA）的Java 开发工具业已亮相。那么Java 工具阵营的老大哥Borland 公司的JBuilder 是否会支持MDA 那？看看他们是怎么说：“我们也在关注MDA， 但是目前仍在观察其动向。比如说第一点，OptimalJ 等产品与JBuilder，包括价格在内，不属于同一类产品。要是支持MDA 的话，Together 更好一些。JBuilder X 在能够轻松构筑Web 应用的角度上，以比这些工具更低的成本实现了相同的功能。同样，即便1 行代码都不写，也能够自动生成可访问数据库的Web 应用架构，在开发过程中及开发完成后均可轻松变更Web 应用服务器等平台。由PIM 向PSM 的转换方法由于将在2004 年实现标准化，因此到时准备在Together 中配备基于MDA 的模型自动生成功能。”看来Borland 公司也不会轻视MDA 这项技术，准备在Together 产品中支持MDA。 MDA 技术是否会取得较大的成功，让我们拭目以待。 下面简单讲述一下从PIM 到PSM 转化的5 种途径： 1. Marking
Excel提取出生及日期格式转换 年龄正确公式 =IF(LEN(B2)=15,YEAR(NOW())-1900-V ALUE(MID(B2,7,2)), IF(LEN(B2)=18,YEAR(NOW())-V ALUE(MID(B2,7,4)),"身份证错")) 根据身份证号码（15位和18位通用）自动提取性别和出生年月的自编公式，供需要的网友参考：说明：公式中的B2是身份证号 1、根据身份证号码求性别： =IF(LEN(B2)=15,IF(MOD(V ALUE(RIGHT(B2,3)),2)=0,"女","男"),IF(LEN(B2)=18,IF(MOD(V ALUE(MID(B2,15,3)),2)=0,"女","男"),"身份证错")) 2、根据身份证号码求出生年月： =IF(LEN(B2)=15,CONCATENATE("19",MID(B2,7,2),".",MID (B2,9,2)),IF(LEN(B2)=18,CONCATENATE(MID(B2,7,4),".", MID(B2,11,2)),"身份证错")) =MID(AF5,7,8) 3、根据身份证号码求年龄： =IF(LEN(B2)=15,YEAR(NOW())-1900-V ALUE(MID(B2,7,2)), IF(LEN(B2)=18,YEAR(NOW())-V ALUE(MID(B2,7,4)),"身份
证错")) 如何将yyyy-mm-dd转为yyyymm:函数法： =TEXT(A1,"emm") 电子表格YYMM转化YYYY-MM-DD ="20"&LEFT(C14,2)&"-"&RIGHT(C14,2)&"-"&"01" 如何将yyyy-mm-dd转为yyyymmdd =text(A1,"emmdd") 在excel中设A1=196303或者19630302，则 B1=TEXT(LEFT(A1,4)&"-"&MID(A1,5,2)&"-"&RIGHT(A1,2 ),"YYYY-MM-DD"),回车即可转换成1963-03或者 1963-03-02，，下拉可以批量转换。 excel日期yyyy-mm-dd转换成yyyy.mm.dd 查找-替换为.单击全部替换按钮。或者用公式：=text(a1,"yyyy.mm.dd")
2013-2014学年第二学期图像通信课程设计报告设计题目：图像的各种颜色空间转换
摘要 所谓三基色原理，是指自然界常见的各种颜色光都可由红、绿、蓝三种色光按照不同比例相配而成。同样，绝大多数颜色也可以分解成红、绿、蓝三种色光。这就是色度学中的最基本的原理。 彩色模型的用途是在某些标准下用通常课接受的方式简化彩色规范。常常涉及到用几种不同的彩色空间表示图形和图像的颜色，以应对不同的场合和应用。因此，在数字图像的生成、存储、处理及显示时，对应不同的彩色空间，需要作不同的处理和转换。现在主要的彩色模型有RGB模型、CMY模型、YUV模型、YIQ 模型、YcbCr模型、HSI模型等。本设计主要使用MATLAB编程的方法，实现RGB与其余四种模型之间的互化。即使用不同的色彩模型表示同一图形或图像。通过转换实现色彩模型的变换之后，可以让同一幅图像以各种模式在全球范围内流通，所以本设计具有一定的实际意义。一般的图像原始都为RGB—加色混合色彩模型，它与剩下的几个色彩模型之间存在着函数对应关系，通过矩阵运算改变模型的参数就可以实现不同色彩模型之间的相互转换。例如CMY—减色混合色彩模型，就是利用青色、深红色、黄色这三种彩色按照一定比例来产生想要的 彩色，CMY是RGB三基色的补色，它与RGB存在如下关系：C M Y = 1 1 1 - R G B ， 使用MATLAB编程时，读入三个通道的数值，按照对应关系进行矩阵变换就可以转换成CMY色彩模型。其他色彩模型转换原理与此相似。 关键词：MATLAB，RGB、YUV、YIQ、YCbCr、HSI、色彩模型
一、设计任务、目的和要求 任务：实现RGB模型、CMY模型、YUV模型、YIQ模型、YcbCr模型、HSI 模型这几种不同色彩模型之间的相互转换 要求：最终结果用图像显示 二、总体方案设计 系统运行环境：WINDOWS 7操作系统 编程软件平台：MATLAB2012b 编码算法原理：将原图的三基色数值读入，根据不同色彩模型之间的相互关系，通过矩阵运算改变不同的亮度和色度等信息来实现色彩模型的转换，然后将变换后的图像导出 流程图： 三、设计实现
坐标转换模型 1.空间直角坐标系间的转换模型（七参数模型） ①公式（布尔莎模型）： ②分析： (1)将O-XYZ中的长度单位缩放l+m倍，使其与O'-X'Y'Z'的长度单位一致； (2)从X反向看向原点O，以O为旋转点，让O-XYZ绕X轴顺时针旋转Wx角，使经过旋转后的Y轴与O'-X'Y'Z’平面平行； (3)从Y反向看向原点O，以O为旋转点，让O-XYZ绕Y轴顺时针旋转Wy角，使经过旋转后的X轴与O'-X'Y'Z'平面平行。显然，此时Z轴也与Z'轴平行； (4)从Z反向看向原点O，以O点为旋转点，O-XYZ绕Z轴顺时针旋转Wz角，使经过旋转后的X轴与X’轴平行。显然，此时O-XYZ的三个坐标轴己与O'-X'Y'Z’中相应的坐标轴平行； 原坐标为O-XYZ，转换到新坐标O-X’Y’Z’.(两坐标系都为空间直角坐标系）其中（dX dY dZ）为坐标原点的平移参数，即将坐标O-XYZ的原点分别沿三个坐标轴平移-dX,-dY,-dZ，使原坐标轴与O-X’Y’Z’的点重合。m为尺度参数，（w1 w2 w3）分别为坐标轴的旋转参量（角度），构成的旋转矩阵分别为： 分别将R1 R2 R3代入上式，可得：
当旋转角度w1 w2 w3很小时（<=10),cos(w)=1,sin(w)=0;在误差允许范围内可以将模型简化为：(同样七参数模型） 四参数模型是在七参数模型的特例，没有考虑坐标轴的旋转量，只考虑坐标轴的平移。 总结： 类似布尔莎模型（以坐标原点为参考点），还有莫洛金斯基坐标模型（以目标点为变换中心）、武测转换模型和范士转换模型（以控制网参考点的站心地平坐标系的三个坐标轴为旋转轴），这些坐标转换模型很容易实现相关坐标在不同坐标系的转换，但是参考位置的偏移向量的相关参数，在实际运用中这些参量是很难测定的，并且受地球重力等物理因素的影响，两个坐标系统即使经过相似变换，仍可能存在较大的残差，所以这些模型适用于简单且规则模型中。 ④程序： clc clear all dX=input('please input value of dX=');
在EXCEL中如何将日期格式转换 如何将格式转换成1900-01-01的格式？ 如果日期所在单元格在A列，并从A1开始。 1、插入一辅助列B列， 2、在B1输入：=MID(A1,1,4)&-MID(A1,5,2)&-MID(A1,7,2)回车， 3、向下填充该公式。 4、选中B列→复制→粘贴→选择性粘贴→数值→确定 5、删除A列 在EXCEL中如何利用身份证号码提取出生时间 假设身份证号在A1，时间格式为1900-1-1 =IF(LEN(A1)=15,"19"&MID(A1,7,2)&"-"&MID(A1,9,2)&"-"&MID(A1,11,2),MID(A1,7,4)&"-"&MID(A1 ,11,2)&"-"&MID(A1,13,2)) 在EXCEL中如何利用身份证号码提取性别 =IF(LEN(A1)=15,IF(MOD(MID(A1,15,1),2)=1,"男","女") , IF(MOD(MID(A1,17,1),2)=1,"男","女") ) 在EXCEL中如何利用身份证号码计算周岁 假设身份证号在C1，时间格式为1900-1-1，出生时间为X 答案1： =DATEDIF(IF(LEN(A1)=15,"19"&MID(A1,7,2)&"-"&MID(A1,9,2)&"-"&MID(A1,11,2),MID(A1,7,4)&"-" &MID(A1,11,2)&"-"&MID(A1,13,2)),TODAY(),”y”) 简化之：=DATEDIF(X,TODAY(),”y”) 答案2： =TEXT(TODAY()-IF(LEN(C1)=15,"19"&MID(C1,7,2)&"-"&MID(C1,9,2)&"-"&MID(C1,11,2),MID(C1,7,4) &"-"&MID(C1,11,2)&"-"&MID(C1,13,2 )),"yy") 简化之：=TEXT(TODAY()-X,”yy”) 在EXCEL中如何利用身份证号码计算周岁 假设身份证号在C1: =IF(C1="","",IF(IF(LEN(C1)=15,"19"&MID(C1,7,2)&"-"&MID(C1,9,2)&"-"&MID(C1,11,2),MID(C1,7,4) &"-"&MID(C1,11,2)&"-"&MID(C1,13,2))>"2009-12-31","",TEXT(TODAY()-IF(LEN(C1)=15,"19"&MID( C1,7,2)&"-"&MID(C1,9,2)&"-"&MID(C1,11,2),MID(C1,7,4)&"-"&MID(C1,11,2)&"-"&MID(C1,13,2))," yy")))
色彩空间介绍 颜色模型是指某个三维颜色空间中的一个可见光子集,它包含某个颜色域的所有颜色。如我们所熟知的三原色光模式.三原色光模式（RGB color model），又称RGB颜色模型或红绿蓝颜色模型，是一种加色模型，将红（Red）、绿（Green）、蓝（Blue）三原色的色光以不同的比例相加，以产生多种多样的色光(如图1所示)。 图1 在大多数的彩色图形显示设备一般都是使用红、绿、蓝三原色，我们的真实感图形学中的主要的颜色模型也是RGB模型，但是红、绿、蓝颜色模型用起来不太方便，它与直观的颜色概念如色调、饱和度和亮度等没有直接的联系。为了更便于颜色的直观表示,一些学者提出了其它的颜色模型,如HSV、HSI、CHL、LAB、CMY等。 RGB颜色模型 RGB（Red,Green,Blue）颜色模型通常使用于彩色阴极射线管等彩色光栅图形显示设备中，彩色光栅图形的显示器都使用R、G、B数值来驱动R、G、B电子枪发射电子，并分别激发荧光屏上的R、G、B三种颜色的荧光粉发出不同亮度的光线，并通过相加混合产生各种颜色。RGB颜色模型称为与设备相关的颜色模型，RGB颜色模型所覆盖的颜色域取决于显示设备荧光点的颜色特性，是与硬件相关的。它是我们使用最多，最熟悉的颜色模型。它采用三维直角坐标系。红、绿、蓝原色是加性原色，各个原色混合在一起可以产生复合色。RGB颜色模型通常采用如图2所示的单位立方体来表示。在正方体的主对角线上，各原色的强度相等，产生由暗到明的白色，也就是不同的灰度值。目前在计算机硬件中采取每一象素用24比特表示的方法,（0，0，0）为黑色，（255，255，255）为白色。正方体的其他六个角点分别为红、黄、绿、青、蓝和品红。