某一地区病例的时空分析心得.doc

利用SatScan软件进行某一地区病例的时空分析
一般步骤：
1、准备数据：进行分析需要准备的数据包括ARCGIS部分和SATSCAN部分。ARCGIS部分主要是分析区域的 .shp文件，而SATSCAN部分主要是分析区域的 .cas、.pop、.geo三种文件。
2、数据的字段属性： .shp文件不涉及字段属性，SATSCAN文件需要按照以下的格式进行字段属性设置才能进行分析：
①Case file: xx.cas
Format: <county> <cases=1> <year> <age group> <sex>
②Population file: nm.pop
Format: <county> <year> <population> <age group> <sex>
③Coordinates file: nm.geo
Format: <county> <x-coordinate> <y-coordinate>
3、向SATSCAN导入准备好的三种数据，分别进行pure spatial, pure temporal, space-time analysis这三种分析，生成分析报告。
4、根据分析结果，在打开的 .shp文件的属性表中添加一个字段，例如病例数目、发病率、发病集中时间段等、然后根据赋给的属性生成一副和疾病有关的专题地图，完成分析。
具体到实例：
1、数据：分析的数据包括长沙市街道的 .shp数据和长沙市H1N1病例的具体数据，但是需要进行一定的修改才可以导入SATSCAN进行分析。
2、制作分析用的SATSCAN数据：原始数据可以通过记事本进行修改然后保存为相应后缀格式的文件。具体来说① .geo文件的制作比较简单，可以从 .dbf文件中获取需要的属性修改成如下的格式：<街道的ID号><Latitude><Longitude>； ② .pop文件的制作可以从长沙市街道数据中提取出必要的信息，修改成如下的格式：<街道的ID><year><population><age group缺省处理><sex缺省处理>; ③ .cas文件相对复杂，需要从长沙市H1N1数据中提取相应的信息，数据设计成如下的格式：<街道的ID><case=1><date><age group缺省处理><sex缺省处理>.
3、导入SATSCAN进行分析，同上在ARCGIS中添加一个字段，生成和H1N1相关的专题地图。
存在的问题：
1、由于软件不认可中文，街道的ID号是采用街道的拼音名还是编号；
2、<date>字段的格式设置为<year day month>;
3、病例所在的位置怎么对应到每一条街道？？很多病例的地址都是写的xx公司、xx街道等等。