摘要： 算法一 #include #include #include struct POINT { int x; int y; }; bool IsLineSegmentCross(POINT pFirst1, POINT pFirst2, POINT pSecond1, POINT pSecond2) { //每个线段的两点都在另一个线段的左右不同侧，则能断定线段相交 //公式... 阅读全文

摘要： 1. 矩阵Y对标量x求导： 相当于每个元素求导数后转置一下，注意M×N矩阵求导后变成N×M了 Y = [y(ij)] --> dY/dx = [dy(ji)/dx] 2. 标量y对列向量X求导： 注意与上面不同，这次括号内是求偏导，不转置，对N×1向量求导后还是N×1向量 y = f(x1,x2,..,xn) --> dy/dX = (Dy/Dx1,Dy/Dx2,..,Dy/Dxn)' ... 阅读全文

摘要： 前几天做了有关gis方面的东西，用到了geos库，发现很多以前不了解的，特此总结下。1 gdal和geos、proj三个联合编译的时候，难免会出现各种问题，源代码有问题就在source中改正，环境问题就换编译环境等等等，release的版本有时候也会出问题，尤其是联合编译…有的版本本身就有bug，只能找到fix文件或者换稳定的版本2 gdal编译的时候eg：unint_64!=double就是vc6.0的编译环境问题，换成vs环境下就搞定了，貌似还有capi/geos_c.h里面size_t的未定义问题，直接注释掉就搞定了3 linux下gdal编译是出现 libgdal.la.lnkscri 阅读全文