Add proj for raw geostationary data when writing to geotiff

[tmlunde]

from datetime import datetime, timedelta


geosSettings = {
                "met10": {'sat':"meteosat", 'num':"10", 'instr':"seviri", 'ir':'IR_108', 'freq': 15 }
}
sat = "met10"
fOut = "tst.tif"

cycles = 1
dateNow = datetime.now()
t = dateNow - timedelta(minutes=dateNow.minute % 15, seconds=dateNow.second, microseconds=dateNow.microsecond) - timedelta(15*cycles/(24.*60.))

geosSettings[sat]['meta'] = GeostationaryFactory.create_scene(geosSettings[sat]['sat'], 
                                             geosSettings[sat]['num'], 
                                             geosSettings[sat]['instr'], time_slot=t)

geosSettings[sat]['meta'].load([geosSettings[sat]['ir']])
img = geosSettings[sat]['meta'].image.ir108()
img.save(fOut)

Results in no proj information is written to tst.tif. The proj (and proj4 string) is present in meteosat10.cfg, but unless the scene is projected first (with foo.project("projection") ), this information is omitted. On my wishlist to include the raw geos projection string when writing to geotiff even if project has not been called.

Best wishes
Torleif

[mraspaud]

@tmlunde thanks for reporting this. The problem is probably that geosSettings[sat]['meta'][10.8].area is None...

[mraspaud]

ok, it's actually the opposite, geosSettings[sat]['meta'].area is None, which is quite logical since there could be channels with different areas in the scene object. A workaround is to do

geosSettings[sat]['meta'].area = geosSettings[sat]['meta'][10.8].area

just after the load...

[tmlunde]

Will indeed work. Maybe the missing proj could be set during image creation as the channel problem would self resolve at this stage. But understand if it should ne kept as a feature.

[mraspaud]

I will mark this issue as a feature request for mpop 2.0.

[tmlunde]

Thanks :-)