摘要： GIS矢量数据模型中很多对象都由多部分组成，点、线、面分别对应多点、多线和多面。在构建基础数据模型开发库的时候，除了要合理的屏蔽各种数据驱动可能带来的数据访问接口不一致问题、合理划分操作等级问题（如从点->线段->线串->多线，环->多边形->多多边形；几何对象、要素、要素类、图层等）还需要在不同的操作等级上注意处理一些较为隐蔽的陷阱，多部分图形与单部分图形就是一个。设想一下，图层级别下，一般的点、线、面层均可以同时容纳点和多点、线和多线、面和多面，可是很多情况下我们处理的是图层级对象的一个子集，可能只是一个要素类，因此在设计要素类的类时，必须注意：多面要兼容面， 阅读全文

摘要： 分享一个科技论文写作的资料，虽然主要是讲Computers and GeoSciences这个期刊，但是还是具有较大的指导性，感觉不错。，原文在此：http://about.elsevier.com/pdf/IAMG-author-workshop.pdf 阅读全文

摘要： 标准C++代码，打开OpenMP编译选项得到debug/release版本的可执行程序A.exe，A.exe是一个socket监听服务，监听端口port的请求。A接收到请求会创建一个新的线程t去调用B.dll中的算法完成计算任务，并返回线程t的执行时间。B.dll中会有查询MySQL数据的操作，查询分为2种类型：批量查询和单个查询。线程t会执行部分OpenMP代码，OpenMP for循环代码部分包含了向mysql数据库查询单个数据的语句。批量查询得到的数据放在内存中，因此可以被线程t所启动的多个OpenMP线程(我是4核8线程的i7CPU，一般我设置启动8个OpenMP线程，正常情况下... 阅读全文

摘要： double cls_random::randomBeta( double alpha, double beta){ /*Johnk's beta generator*/ double u, v; double x, y; do { u=cls_random::randomUniform(); v=cls_random::randomUniform(); x=pow(u,1/alpha); y=pow(v,1/beta); } while (x+y>1); return x/(x+y... 阅读全文

摘要： double cls_random::randomRayleigh(double sigma){ double pv = sigma*sqrt(-2 * log(cls_random::randomUniform())); return pv;}期望:E=方差:V=wiki:http://zh.wikipedia.org/wiki/%E7%91%9E%E5%88%A9%E5%88%86%E5%B8%83 阅读全文

摘要： long cls_random::randomPossion(double lambda){ double x = -1, u; double log1, log2; log1 = 0; log2 = -lambda; do { u = cls_random::randomUniform(); log1 += log(u); x++; }while(log1 >= log2); return x; }期望:E=lambda方差:V=lambdawiki:http://zh.wikipe... 阅读全文

摘要： View Code double cls_random::randomGamma( double alpha, double lambda) { double u, v; double x, y; double b, c; double w, z; bool accept = false; double t; if (alpha > 1.0) { /* Best's rejection algorithm XG for gamma random variates (B... 阅读全文

摘要： double cls_random::randomX2(int k){ double dv = 0.0; for (int i=0;i<k;i++) { double p_dv = cls_random::randomGaussian_MarsagliaAndBray(0.0,1.0);//默认生成符合标准正态分布的随机数 dv += p_dv*p_dv; } return dv;}期望：E=k方差：V=2kwiki:http://zh.wikipedia.org/wiki/%E5%8D%A1%E6%96%B9%E5%88%86%E5%... 阅读全文

摘要： View Code long cls_random::randomGeometric( bool shifted, double probability){ long rnd = 0; while(true) { rnd++; double pV = (double)rand()/(double)RAND_MAX; if (pV<probability) { break; } } if (!shifted) { rnd--; }... 阅读全文

摘要： long cls_random::randomNegativeBinomial( double r, double probability){ long rnd = 0; double fr = 0.0; while(true) { double pV = (double)rand()/(double)RAND_MAX; if (pV<probability) { rnd++; } else { fr++; ... 阅读全文