需求大抵是这样的，客户用igraph进行社会网络分析，然后生成了xml格式的数据文件，但并没有很好的网络展示平台，manager就来找我，希望通过我们研发的WebGIS系统展示。在我们的WebGIS系统中，矢量数据都是保存在PostgreSQL数据库里，一般来说，都是通过shapefile 生成sql files，然后导入数据库。由于没有shapefile，但是分析graphicML文件发现结构还是很简单的，而且对于社会网络的两个node之间的关系（relationship），从几何的角度讲，两个nodes构成的关系就是一条两点直线（edge）。基于此，我希望通过python对graphML进行分析，然后直接生成sql file，导入到postgreSQL里面，此后就能用uDig就能直接转换成为shapefile，我不清楚新版的ArcGIS 9.3有没有相似的功能，好久不用Arctoolbox了。

原始数据片段：

< ?xml version="1.0" encoding="UTF-8"?>
<graphml xmlns="http://graphml.graphdrawing.org/xmlns"
         xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
         xsi:schemaLocation="http://graphml.graphdrawing.org/xmlns
         http://graphml.graphdrawing.org/xmlns/1.0/graphml.xsd">
<!-- Created by igraph -->
  <key id="filestub" for="graph" attr.name="filestub" attr.type="string"/>
  <key id="id" for="node" attr.name="id" attr.type="string"/>
  <key id="xcoord" for="node" attr.name="xcoord" attr.type="double"/>
  <key id="ycoord" for="node" attr.name="ycoord" attr.type="double"/>
  <key id="h_loc" for="edge" attr.name="h_loc" attr.type="string"/>
  <key id="h_zip" for="edge" attr.name="h_zip" attr.type="string"/>
  <key id="h_city" for="edge" attr.name="h_city" attr.type="string"/>
  <key id="t_loc" for="edge" attr.name="t_loc" attr.type="string"/>
  <key id="t_zip" for="edge" attr.name="t_zip" attr.type="string"/>
  <key id="t_city" for="edge" attr.name="t_city" attr.type="string"/>
  <key id="appyear" for="edge" attr.name="appyear" attr.type="string"/>
  <key id="patent" for="edge" attr.name="patent" attr.type="string"/>
  <key id="pat_type" for="edge" attr.name="pat_type" attr.type="string"/>
  <graph id="G" edgedefault="undirected" parse.nodes="1022" parse.edges="9399">
    <data key="filestub">/tmp/VDC/DSB/iGraph.322801</data>
    <node id="n0">
      <data key="id">20190035</data>
      <data key="xcoord">394.439</data>
      <data key="ycoord">245.759</data>
    </node>
    <node id="n1">
      <data key="id">21810055</data>
      <data key="xcoord">452.308</data>
      <data key="ycoord">-49.6771</data>
    </node>
    <edge source="n1020" target="n1021">
      <data key="h_loc">MA</data>
      <data key="h_zip">01532</data>
      <data key="h_city">NORTHBORO</data>
      <data key="t_loc">MA</data>
      <data key="t_zip">01601</data>
      <data key="t_city">WORCESTER</data>
      <data key="appyear">2004</data>
      <data key="patent">7459547</data>
      <data key="pat_type">U</data>
    </edge>
  </graph>
</graphml>

通过python对XML进行解析，我还是向大家推荐pyQuery，大概是喜欢Jquery的缘故吧，我觉得pyQuery用起来很方便，很顺手。大家可以到pyquery.org查找帮助和教程。需要注意的是，pyQuery对于命名空间的解析似乎不是很好，我记得以前我写过一篇帖子叫做两个GIS编程技巧，里面提及了用@nodeName=namespace:elementName的方式来对异域命名空间的元素进行解析。但是pyQuery似乎不能进行该操作，所以我在进行解释之前，将根节点的所有命名空间都删除了。
解析本省是一件很简单的事情，随后写入到另一个文件中，其次对每一个node和edge进行分析，按照SQL语法，生成INSERT INTO的语句。

python 代码片段

from pyquery import PyQuery as pq
from lxml import etree
d = pq(filename="data/rna.xml")
file = open('data/output.sql', 'w+')
content = ''
nodes =  d('node')
for item in nodes:
  content += "INSERT INTO sn_node ( the_geom, nid) VALUES ( GeomFromText('POINT(" + item[1].text.__str__() + " " + item[2].text.__str__() +")', 900913), '" + item.attrib['id'] + "');\n"
file.write(content)
file.close()

从而生成output.sql,如下片段，

INSERT INTO sn_node ( the_geom, nid) VALUES ( GeomFromText('POINT(394.439 245.759)', 900913), 'n0');
INSERT INTO sn_node ( the_geom, nid) VALUES ( GeomFromText('POINT(452.308 -49.6771)', 900913), 'n1');
INSERT INTO sn_node ( the_geom, nid) VALUES ( GeomFromText('POINT(172.295 72.8809)', 900913), 'n2');
INSERT INTO sn_node ( the_geom, nid) VALUES ( GeomFromText('POINT(-217.813 -581.18)', 900913), 'n3');
INSERT INTO sn_edge ( the_geom, snid, tnid, s_zip, t_zip, appyear, patent, pat_type) VALUES ( GeomFromText('LINESTRING(-206.829 790.104,-200.639 791.055)', 900913), 'n121', 'n365', '05403', '05452', '1997', '5902044', 'U');
INSERT INTO sn_edge ( the_geom, snid, tnid, s_zip, t_zip, appyear, patent, pat_type) VALUES ( GeomFromText('LINESTRING(-206.829 790.104,-231.181 757.468)', 900913), 'n121', 'n214', '05403', '05403', '1998', '6134704', 'U');

生成了SQL的INSERT INTO 语句后，然后加上创建数据库和构建索引的语句

CREATE TABLE sn_node (
id serial,
nid VARCHAR(6)
);
SELECT AddGeometryColumn('', 'sn_node','the_geom',900913,'POINT',2);
CREATE INDEX sn_node_index ON sn_node USING GIST ( the_geom );

所有的SQL语句组成了一个SQL文件，如果你有PgAdminIII，导入不是繁琐的事情。不再熬述了。uDig可以直接读取PostgreSQL里面的数据，然后我在uDig里面可以绘制出社会网络的模型，也可uDig输出shapefile，再用ArcGIS读取并展示。如图所示：

This application was demonstrated as part of the GIS and Web 2.0 Application Development Demo Theater for the ESRI UC 2009.

The sample illustrates how to:

1. Integrate the ArcGIS Silverlight API with Twitter
2. Geocode your friends with Bing asynchronously
3. Use control templates and animation in XAML
4. Construct point and line graphics 
5. Consume WCF services

中国古代就有吃不到葡萄就说葡萄酸的传统，但是我并不是想说葡萄酸，只是想公正的思考一下为什么这两个人能够获得胜利。Alper做的是将ArcGIS JS API 和EXT结合的构架。说句实话，我真的不觉得这个应用有什么创新，这种概念在2006年时就有了geoEXT，更不说现在的Mapfish。同时，在我所做的应用中，是很少用地图JS API 和EXT进行结合使用的，把这个两个都疯狂耗费客户端内存的东西放在一起使用，是一般用户所难以接受的，虽然界面美观，但是从另一个角度说，你的界面也很单一。在编程方面，我不是个爱慕虚荣的人，所以我不会为了UI而消耗大量的内存。

其次，第二个来自Matthew Petre的，足足四个人的团队，，为什么做出的东西就像是在用simple viewer team的东西呀。。而且，功能很简单，我真的后悔没有把我去年当internet GIS助教时给学生们做的实例拿出来当参赛作品。

也许这就是ESRI的风格，昨天我忍不住骂了ESRI的ArcGIS，今天又不得不颇感无奈，也许自己钟情的开源和ESRI本来就是两条路，想要ESRI支持开源，给你点financial support，估计是相当impossible的stuff。

anyway,最后David Martinez告诉我，说他们明年也要开发heatmap，简直是心头再砍一刀，是自己的思路太超前了？？还是我抢走了ESRI 工程师的饭碗，做了他们应该做的事情？ 呵呵， 就当这次是和ESRI交了一次手，大哥，俺领教您了。下次，俺会长个心眼。

最后还是要感谢这次帮助我的兄弟姐妹们，每每我在QQ或者MSN上倡议帮俺投票，兄弟们总是义无反顾，以后在道上混，还需要大家多多照顾：）

今天很想用geoprocessing做点事情，毕竟illir老师希望我们基于ArcGIS做点相应的开发，加之自己以前做spatial analyst还算颇有经验。谁知arcgisscripting的 version hell真的让人崩溃，俺总不能为了满足arcgis的需要删除shapely，numpy，scipy以及python for GDAL吧。。

一句长叹，我本将心向明月，奈何明月照沟渠。

Hi Guys！ I participate in the 2009 ArcGIS Server Mashup Code Challenge, I need your help to vote for my mashup work at http://surveys2.esri.com/mashup, my work’s title is HeatMap Plugin for ArcGIS Flex API, Flex, Bo Zhao!!! Thank you very much!

1. /lab/heatmap/wifi.html
2. /lab/heatmap/live_earthquakes.html.

And you can download this plugin at http://resources.esri.com/arcgisserver/apis/flex/index.cfm?fa=codeGalleryDetails&scriptID=16048 .

PS: Hope you could support my heatmap plugin. And just vote for me if you think this plugin is useful. And Now I am trying to convert it to meet Umap.

datasource: feed://earthquake.usgs.gov/eqcenter/catalogs/7day-M2.5.xml

I decided to make a plugin for ArcGIS Flex API first to make heatmaps. It might cost me couple of days, however, it will faciliate to make heatmaps then.

Paul said that If you want people to switch you need to make the transition as painless as possible. Firefox got people to switch to IE by

• *Making better software
• *Not making user learn a new UI for interacting with the web
• *importing all their IE favorites
• *THEN building in cool new features that keep people around

So, uDig Team have to do is make something better than ArcGIS, but not so much better that it is not familiar to the existing user base, that works transparently with all their existing data and presumably their .mxd files too. And meanwhile give it away for free.

Paul said, Open source is not about users, it is about developers. It is only about users in so far as users become sufficiently engaged in the project that they either become developers themselves, or support developers through careful bug finding or documentation.

At last, Paul describe the correct models by his point of view. He said that the FireFox or Open Office are incorrect modles, because they are something a litte aberrate. Paul said the correct models are Linux, Apache, Perl, PostgreSQL，not user-friendly, but still very useful.

如果影像文件没有金字塔或者金字塔文件无法访问，ArcReader并不能自己创建金字塔。在ArcGIS中，金字塔文件对于影像的绘制性能非常重要。通常，在ArcMap的作者使用栅格时会创建金字塔文件，该文件被放到与栅格文件同一个文件夹。如果ArcMap的作者没有对栅格文件所在文件夹的写权限，金字塔文件就会创建到该作者的临时文件夹下。那么如果在其它电脑上的ArcRead的用户打开一个带栅格的PMF文件，就无法访问到金字塔文件了。没有金字塔文件，绘画的性能就很差，质量也不理想。出现这种情况时，可以将栅格的金字塔复制到本地的temprasterproxies文件夹，但是最好还是将金字塔和栅格文件一起放在一个单独的源文件夹下。

 这时，我们用ArcGIS的hillshade命令，将DEM数据处理成hillshade图。选择一定的太阳高度角和方向角。如下面的左图。这样，把hillshade层放到用过渡色表示的DEM数据的下面，然后再给他赋一个透明度。赋透明度的目的是为了将将DEM和hillshade进行叠加。经过试验，透明度为30%，效果比较好。如下右图所示。

呵呵，看到了吧，效果还是不错滴。其实大家会说Musicman不关心技术，只关心新闻，其实musicman每天都要看很多资料，每当夜深人静的时候，一人在床前，盯着自己的X40，就有说不尽的话，哎，真的很想出国，有时候发现自己留在国内是个错误，但是留下了就努力取干吧，与其在蛮怨中生活，还不如自己坚持取努力，让自己变得更加强悍。