在Ubuntu上安装GIS软件，可以添加这个PPA源：
deb http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu karmic main
导入GPG KEY
gpg –keyserver subkeys.pgp.net –recv 089EBE08314DF160
gpg –export –armor 089EBE08314DF160 | sudo apt-key add -
这个源中目前包含了最新版本的GIS软件：
针对地理数据的ORM，有一个Hibernate的扩展HibernateSpatial项目可以将JTS对象映射到MySQL/PostGIS/Oracle中。这个扩展同样可以用在Grails里，这里有一篇简单的介绍，关于在Grails和MySQL中管理地理数据：
http://www.grails.org/MySQL+GIS-Geometry+with+Grails
不过按照这个文章里介绍的方法用，很可能会遭遇这样的报错：
org.hibernate.MappingException: No Dialect mapping for JDBC type: 2003
这个问题最终在这里得到了解答：
http://n2.nabble.com/No-Dialect-mapping-for-JDBC-type-2003-td1141106.html
按照邮件列表里的反映，上面的配置在Postgis里是可以work的，但是如果用Mysql还需要指定JPA的columnDefinition，对应的Hibernate属性是sql-type。虽然作者承诺会在今后的版本里修改这个问题，不过眼下的M2版本还没有修正这个问题。为此，Grails的用户特地提出在Grails中加入sql-type的支持：
http://jira.codehaus.org/browse/GRAILS-3201
现在按照下面文档的说明，可以在mapping里指定sqlType了：
http://grails.org/doc/latest/ref/Database%20Mapping/column.html
实例代码里的domain定义应该改成：
import com.vividsolutions.jts.geom.Polygon
import org.hibernatespatial.GeometryUserType
public class MyPoly {
String name
Polygon poly
static mapping = {
poly type: GeometryUserType， sqlType:"GEOMETRY"
}
}
于是，再也没有莫名其妙的No Dialect报错了。
废话不说了，直接上code吧：
import com.drew.imaging.jpeg.*;
import com.drew.metadata.*;
import com.drew.metadata.exif.*;
file = new File("sample.jpg");
meta = JpegMetadataReader.readMetadata(file);
gpsdir = meta.getDirectory(GpsDirectory.class);
lat = gpsdir.getRationalArray(GpsDirectory.TAG_GPS_LATITUDE);
lon = gpsdir.getRationalArray(GpsDirectory.TAG_GPS_LONGITUDE);
lats = lat[0].doubleValue() + lat[1].doubleValue()/60 + lat[2].doubleValue()/3600;
lons = lon[0].doubleValue() + lon[1].doubleValue()/60 + lon[2].doubleValue()/3600;
println(lats);
println(lons);
metadata-extractor似乎是现在惟一的读取exif信息的java库。在ivy中可以添加这样一个依赖：
<dependency org="com.drewnoakes" name="metadata-extractor" rev="2.4.0-beta-1" conf="runtime"/>
续昨，虽然昨天升级到了gnome2.28，也安装了libchamplain，但是却没有发现这个库应用。Ubuntu上的Empathy似乎是没有把location support编译进去。看不到效果怎么办，好在libchamplain的作者还提供了另一个应用。这是一个eog的插件，用于读取照片的exif信息中的gps latitude和longitude，进而将他显示在libchamplain的地图里。
安装这个插件可以follow官方网站上的做法，当然开发包、头文件之类的东西要齐备。
libchamplain使用了clutter，在载入地图的时候还有fadein的效果。clutter是将要应用在下一代桌面gnome shell中的图形特效库。
插件是用C写的，可以看成是libchamplain的一个例子，作者说已经有了其他语言的binding，如此未来在gnome桌面上开发tile地图程序会变得简单。
今天偶然看到一个Ruby的GeoHash库，功能非常简单，就是将经纬度坐标哈希成字符串，并且，利用这个字符串可以粗略地比较两点的距离。这个库的代码host在GitHub上，要安装，需要将github加入gem：
sudo gem sources -a http://gems.github.com
sudo gem install davotroy-geohash
用法，irb：
irb(main):001:0> require ‘rubygems’
=> true
irb(main):002:0> require ‘geohash’
=> true
irb(main):003:0> GeoHash.encode(32.168, 118.54)
=> “wtsr12n0nj”
irb(main):004:0> GeoHash.decode(‘wtsr12n0nj’)
=> [32.168, 118.54]
代码可以从github上获得
git clone git://github.com/davetroy/geohash.git
GeoHash的核心部分用C实现，仔细看一下可以发现是一个四叉树的结构：
[codesyntax lang="c"]
#define BASE32 "0123456789bcdefghjkmnpqrstuvwxyz"
static void encode_geohash(double latitude, double longitude, int precision, char *geohash) {
int is_even=1, i=0;
double lat[2], lon[2], mid;
char bits[] = {16,8,4,2,1};
int bit=0, ch=0;
lat[0] = -90.0; lat[1] = 90.0;
lon[0] = -180.0; lon[1] = 180.0;
while (i < precision) {
if (is_even) {
mid = (lon[0] + lon[1]) / 2;
if (longitude > mid) {
ch |= bits[bit];
lon[0] = mid;
} else
lon[1] = mid;
} else {
mid = (lat[0] + lat[1]) / 2;
if (latitude > mid) {
ch |= bits[bit];
lat[0] = mid;
} else
lat[1] = mid;
}
is_even = !is_even;
if (bit < 4)
bit++;
else {
geohash[i++] = BASE32[ch];
bit = 0;
ch = 0;
}
}
geohash[i] = 0;
}
[/codesyntax]
可惜在学校的时候没有看到这么精致的代码。