<![CDATA[z.kklt的记事本 - 常用资料]]> http://www.z.kklt.net/ zh-cn PBlog2 v2.4 z.kklt的记事本 http://www.z.kklt.net/images/logos.gif http://www.z.kklt.net/ z.kklt的记事本 http://www.z.kklt.net/article/常用资料/282.htm <![CDATA[LISP的10进制到16进制转换的函数]]> z.kklt@163.com(z.kklt) Thu,01 Apr 2010 00:36:11 +0800 http://www.z.kklt.net/default.asp?id=282   比如10进制的30转16进制：
    1、30/16，商为1，余为14，16进制的E；
    2、商1/16，商为0，余为1；
    3、将余数逆序，为1E。

    LISP下的10进制转16进制的函数，比较简单，但是有用：

(defun tth (t0 / tl tf i j hd hl l16)
  (setq tf (fix (/ t0 16))
    hl ""
    l16 (list (cons 10 "a") (cons 11 "b") (cons 12 "c") (cons 13 "d") (cons 14 "e") (cons 15 "f")))
  
  (while (/= tf 0)
    (setq tl (append tl (list (rem t0 16)))
      t0 tf
      tf (fix (/ t0 16)))
    )
  (setq tl (append tl (list (rem t0 16))))

  (setq i (length tl)
    j (1- i))
  (repeat i
    (setq hd (nth j tl)
      j (1- j))
    (if (>= hd 10)
      (setq hl (strcat hl (cdr (assoc hd l16))))
      (setq hl (strcat hl (itoa hd))))
    )  
;;;  (princ tl)
;;;  (princ hl)
;;;  (princ "\n")
  )

示例：
(tth 3246)
"cae"
]]> http://www.z.kklt.net/article/常用资料/276.htm <![CDATA[LISP写的高斯投影正反算及七参数投影变换相关程序]]> z.kklt@163.com(z.kklt) Sun,14 Feb 2010 17:08:45 +0800 http://www.z.kklt.net/default.asp?id=276 主要功能：
1、内置54、80、84、2000椭球基本参数
2、上述椭球的标准高斯正反算；
3、上述椭球上具有高程抵偿面的任意带高斯投影正反算；
4、7参数投影变换；
5、上述椭球的大地坐标和空间直角坐标的转换；

辅助功能：
1、D.MMSS格式的度分秒和弧度的互相转换；
2、54、80、84、2000椭球的基本参数及计算常用参数的获取；

利用这些可以做的其他事情：
1、GPS测得的WGS84坐标和地方坐标的转换；
2、换带计算；
3、经过组合，可以实现不同坐标系、不同坐标形式的自动转换。

程序下载z.kklt_cotr.rar:


程序需要在AUTOCAD下运行，运行说明如下：
;**************************************************
;;DUTOHU
;将D.MMSS格式的角度转换为弧度
;用法:
;(dutohu du)

;**************************************************
;;HUTODU
;将弧度计的角度转换为D.MMSS格式的角度
;用法：
;(hutodu hu)

;**************************************************
;;GPP
;获取椭球基本参数,与位置无关
;用法:
;gpp (gn gp)
;gn为椭球名称,可接受54、80、84、2000
;gp为待获取的参数,可接受a、b、e、e1
;如:
;(gpp 54 a)
;意为取54椭球长半轴a的值

;**************************************************
;;GPF
;获取椭球的辅助参数,与位置有关
;gn椭球名
;gp待获取的参数
;gb点的大地纬度B
;gl点的大地经度L
;gu点的单位格式,1为度.分秒;2为十进制度;3为弧度

;gp选择
;w
;v
;dn
;n卯酉圈曲率半径
;r平均曲率半径

;**************************************************
;;BLHTOXYZ
;将大地坐标转为空间直角坐标
;旋转角单位为弧度
;输入大地坐标B、L、H
;用法为：
;(blhtoxyz gn blhb blhl blhh)
;gn可选54、80、84、2000
;输出空间直角坐标，分别为xyzx，xyzy，xyzz，全局变量

;**************************************************
;;XYZTOBLH
;将空间直角坐标转为大地坐标
;输入空间直角坐标
;用法为：
;(xyztoblh gn xyzx xyzy xyzz)
;输出大地坐标B、L、H，分别为BLHB、BLHL、BLHH，全局变量，以弧度计

;**************************************************
;;BURSA
;空间直角坐标七参数转换
;用法为：
;(bursa tx ty tz wx wy wz m gn1x gn1y gn1z)
;tx、ty、tz为三个平移参数，单位为米
;wx、wy、wz为三个旋转参数，单位为弧度
;m为尺度比差异，如0.000004，单位为1，不是ppm
;输出gn2x、gn2y、gn2z为转换后的空间直角坐标，外部变量

;**************************************************
;;BLTOXY
;高斯正算
;子午线弧长lb、xyx、xyy为外部变量，其余为内部变量
;单位均为弧度
;用法为：
;(bltoxy gn gl0 blb bll)
;gn可选54、80、84、2000
;gl0为中央经度
;输入blb、bll为大地经纬度
;输出xyx、xyy为高斯正算后的坐标，全局变量

;**************************************************
;;XYTOBL
;高斯反算
;底点纬度bf、blb、bll为外部变量，其余为内部变量
;单位均为弧度
;用法为：
;(xytobl gn gl0 xyx xyy)
;gn可选54、80、84、2000
;gl0为中央经度
;输入xyx、xyy为平面坐标
;输出blb、bll为高斯反算后的大地经纬度，全局变量，以弧度计

;**************************************************
;;DFBLTOXY
;具有高程抵偿面的任意带高斯投影的高斯正算
;用法为：
;(dfbltoxy gn b0 l0 dx dy m blb bll)
;gn可选54、80、84、2000
;b0为中央纬度
;l0为中央经度
;dx为X加常数
;dy为Y加常数
;m为尺度比
;输入blb、bll为待解算的大地经纬度
;输出xyx、xyy为高斯正算后的坐标，全局变量

;**************************************************
;;DFXYTOBL
;具有高程抵偿面的任意带高斯投影的高斯反算
;用法为：
;(dfxytobl gn b0 l0 dx dy m xyx xyy)
;gn可选54、80、84、2000
;b0为中央纬度
;l0为中央经度
;dx为X加常数
;dy为Y加常数
;m为尺度比
;输入xyx、xyy为平面坐标
;输出xyx、xyy为高斯反算后的大地经纬度，全局变量，以弧度计
]]> http://www.z.kklt.net/article/常用资料/246.htm <![CDATA[二次曲面拟合的副产品：用LISP写的矩阵乘法程序]]> z.kklt@163.com(z.kklt) Mon,14 Sep 2009 22:12:52 +0800 http://www.z.kklt.net/default.asp?id=246
    EXCEL中有现成的矩阵乘法的公式，具体用法为：MMULT(ARRAY1,ARRAY2),选定结果区域，再按Ctrl+Shift+Enter。

    用EXCEL试验下算法可以，要批量转换，还得自己写程序。莫法，只会点LISP，但是LISP里面没有现成的矩阵乘法函数，只有自己写了。

    先复习一下矩阵乘法的算法：
    M行N列矩阵乘N行K列矩阵，得到M行K列矩阵，具体算法为：第一个矩阵的第一行分别乘以第二个矩阵的第一列，再加起来，记在第一行第一列的位子，再把第一个矩阵的第一行分别乘以第二个矩阵的第二列，再加起来，记在第一行的第二列位子上。。。一直乘到K列，得到乘积矩阵的第一行；再换第一个矩阵的第二行、第三行。。。直到第M行，最终得到M行K列矩阵。
    简单说来，就是乘积矩阵中，第I行J列的元素等于第一个矩阵的第I行乘于第二个矩阵的第J列，再求和。

;数组乘法,arraymn与arraynk相乘,得到arraymk
;m行n列数组与n行k列数组相乘，得到m行k列数组
;arraymn(3*2)
;((m1n1 m1n2) (m2n1 m2n2) (m3n1 m3n2))
;arraynk(2*2)
;((n1k1 n1k2) (n2k1 n2k2))
;数组的第m个元素序列为数组的第m行，第m个元素序列的第n个元素为数组的第m行第n列的数据


(defun mmult (arraymn arraynk / am an1 an2 ak ami ani aki)
...
)

用法：
定义ARRAYMN:
(setq arraymn (list (list 1.1 1.2) (list 2.1 2.2) (list 3.1 3.2)))

定义ARRAYNK:
(setq arraynk (list (list 1.1 1.2 1.3) (list 2.1 2.2 2.3) ))

矩阵相乘:
(mmult arraymn arraynk)

乘积矩阵ARRAYMK:
((3.73 3.96 4.19) (6.93 7.36 7.79) (10.13 10.76 11.39))

LISP写的矩阵乘法程序：


]]> http://www.z.kklt.net/article/常用资料/231.htm <![CDATA[推荐一个程序界面设计软件Balsamiq Mockups]]> z.kklt@163.com(z.kklt) Tue,21 Apr 2009 14:35:59 +0800 http://www.z.kklt.net/default.asp?id=231
    这个是官网地址：http://www.balsamiq.com/products/mockups/desktop

      到官网看了下，这东东需要AdobeAIR环境，在官网就有下载，把AIR和Mockups都下下来，试用了一下，画了一个简单的界面设计图，出来的图很简洁，手绘风格很喜欢！

    未注册版本不能保存，隔5分钟弹出提示框。这软件的价格是79美刀，不过有一些license获取方案，比如在BLOG推荐，赞一个！

    注：软件是英文的，刚开始的时候不能输入中文，后来选择使用系统字体就可以了，不过看起来没英文的漂亮，稍稍遗憾。

    这是我用DEMO版画做的一个简单的界面设计，试用版无法保存，是截屏下来的：

]]> http://www.z.kklt.net/article/常用资料/230.htm <![CDATA[GPS功分器，让两台接收机接收一个GPS天线的信号]]> z.kklt@163.com(z.kklt) Wed,15 Apr 2009 19:00:33 +0800 http://www.z.kklt.net/default.asp?id=230     ALDCBS 1X2：


    再一搜，居然发现国内有卖这东西的公司，估计这玩意不便宜，毕竟是进口货，还是给那公司发了邮件索取报价，那边反应倒快，隔天就回复了，只是价格要4400，只得作罢，另外找办法。研究了说明书，发现这个GPS功分器的设定是J1端口给GPS天线和功分器供电，J2端口接了一个200欧的电阻，模拟天线终端，作用的频率范围是1.1GHz-1.7GHz。

    ALDCBS 1X2的关键参数：


    关于给天线的供电问题，电子市场的功分器多半没设计这么复杂，估计是直通的，也就是两台接收机并联给GPS天线供电，用起来应该没问题，这个因素不考虑；只考虑频率范围就好办了，到电子市场搞定，又配了几个转接头，总花费不到200就搞定了。
    代用的GPS功分器：

    回来试验，先只接了一台接收机，没问题，搜星正常，再接一台接收机，几秒钟后，接收机的卫星灯的红灯开始闪烁，再一会，隔一秒闪烁一次，没问题。又试了下存储数据，存了10多分钟数据，解算了一下，一切正常（实际是0基线测试，在另一篇里会讲到）。
]]> http://www.z.kklt.net/article/常用资料/229.htm <![CDATA[利用TTC后处理生成虚拟站点RINEX数据]]> z.kklt@163.com(z.kklt) Tue,14 Apr 2009 23:40:04 +0800 http://www.z.kklt.net/default.asp?id=229     由于平时处理基线和平差都是基于TGO，在TTC中，主要关注TGO中没有的功能，最近在研究后处理VRS模块，利用该模块，可以生成虚拟站点的RINEX观测数据。

    先看看TTC（Trimble Total Control）的帮助文件中对此的简单介绍：
    VRS虚拟参考站
    根据区域网所提供的实时改正所产生的离流动站仅几米的参考站。从而，极大地提高RTK定位的性能。
    使用后处理数据，Trimble Total Control的VRS模块可以为用户网中的每个流动站产生相似的虚拟参考站。
    
      该功能非常有用，只要有参考站的观测数据，可以虚拟任意一点的观测值，但是由于规范中没有对这项技术作出相关规定，没办法正式使用（虽然是2003年的技术，但是对于规范而言，还是太前卫了点，想想看，明明没在这个点摆仪器，居然得到这个点的观测数据了）。不过研究一下、增长增长见识还是很有好处的。

    在GPSNet系列软件中，有一个和GPSWeb集成度很高的服务软件GPServer，利用该软件，可以为用户提供参考站网内任意一点、任意时间的虚拟观测数据，我觉得TTC中的后处理VRS模块就是它的前身。

    过程很简单，导入几个站点的观测数据，输入精确的WGS84坐标，将这些点指定为控制点，然后选择"VRS-产生VRS"，指定移动站位置、观测时间段、采样间隔就可以了。



      为了保证能正确生成RINEX数据，最好选择稍长时间段，比如上面就生成了24小时的观测文件。在使用时，可以用TEQC软件对这些数据进行分割，这里是取了前120分钟的观测值。另外，TTC生成的RINEX只含了观测数据，不包括导航电文，需要另行准备，可以指定已知测站的导航电文。



      导入TGO的时候需要手动指定导航文件，也可以将准备好的导航电文改成和观测数据同名，就不需要手动指定了。



      解算基线


      网图

      
      这个是生成的3个站点的120分钟的观测数据文件（O文件）和用TEQC分割的对应的导航文件（N文件，导入TGO需要手动指定到该导航文件），供有需要的朋友下载，设置个注册限制吧，无它，希望朋友多来看看，也多多交流。
TEST_RINEX.rar：


]]> http://www.z.kklt.net/article/常用资料/192.htm <![CDATA[TGO导出网图的小程序]]> z.kklt@163.com(z.kklt) Mon,17 Nov 2008 12:14:39 +0800 http://www.z.kklt.net/default.asp?id=192     具体的操作方法是：在TGO里导出Trimble交换文件ASC文件；在CAD里加载本程序，运行命令为CDKCATD，选择合适的比例，选择ASC文件即可生成网图！由于是自用，可能有不完善的地方，也未做图框！

    这是TGO里面的网图：


     这是在CAD里面生成的网图：


这是编译过的程序：下载文件 cdkcatd.fas.rar

这是源代码：

]]> http://www.z.kklt.net/article/常用资料/188.htm <![CDATA[TRIMBLE NETRS接收机1.2版固件]]> z.kklt@163.com(z.kklt) Mon,10 Nov 2008 15:48:17 +0800 http://www.z.kklt.net/default.asp?id=188 点击下载此文件
]]> http://www.z.kklt.net/article/常用资料/181.htm <![CDATA[具有高程抵偿面的任意带高斯正形投影正算的补充公式]]> z.kklt@163.com(z.kklt) Tue,04 Nov 2008 11:53:27 +0800 http://www.z.kklt.net/default.asp?id=181

    一般的测量学书上列出了高斯投影正算公式，均是尺度为1，原点纬度为0的情况，而一般地方坐标系的尺度不为1，原点纬度是自定义的，需要补充公式.

    具有高程抵偿面的任意带高斯正形投影正算的补充公式：

    X = mX54 - mX54B0 + X0

         Y  =  mY54 + Y0

式中：  

m为地方坐标系的投影比例;

         X54为待求点按高斯投影正算公式计算的X坐标;

         X54B0为坐标原点按高斯投影正算公式计算的X坐标;

         X0 为坐标原点的X加常数;

    Y54为待求点按高斯投影正算公式计算的Y坐标;

     Y0 为坐标原点的Y加常数.

]]> http://www.z.kklt.net/article/常用资料/176.htm <![CDATA[GPS年积日计算]]> z.kklt@163.com(z.kklt) Fri,10 Oct 2008 10:42:19 +0800 http://www.z.kklt.net/default.asp?id=176
格里高利历转年积日GPS时小软件：下载文件 点击下载此文件

原文地址：
http://aclog.cn/log/?/61.aspx

原文：
GPS数据处理时常需要用到IGS站的观测文件和精密星历。因为两种文件的命名是以不同的时间标示法来取的。格里高利历也就是平时我们所认识的日历，而年积日和GPS时的时间标示法不够直观。在下载IGS站的观测文件和精密星历文件时需要一个小小的计算或者在Linux下用命令进行换算，这些貌似有些繁琐。所以就写了这么一个小的转换软件，它能将输入的格里高利历转换成对应的年积日和GPS时，这样很方便下载对应的IGS文件。

例子：输入2008 6 30 结果为182  1486 1

第一个数182是2008年6月30日对应的年积日，也就是说这一天是2008年的第182天。

第二个数1486是2008年6月30日对应的GPS时，一般也称第1486周。这是从1980年1月6日0时开始起算的周数。

第三个数1表示2008年6月30日这一天是星期一，这个是从星期天起算的，也就是0~6。数字是几就表示星期几。

OK！

]]>