// raster2floatDEM.c
//

#include <stdio.h>
typedef unsigned int guint;
typedef unsigned char guint8;

#include "../Moon_128.c"
// https://astrogeology.usgs.gov/search/map/Moon/LRO/LOLA/Lunar_LRO_LOLA_Global_LDEM_118m_Mar2014
int main(int argc, char* argv)
{
	unsigned int i,ii,j,k,count;
	float h, h0;
	guint8 p;
	//we'll flip the image coord, so we can start at the lower-left and work up
	//  with +xSpacing, +zSpacing
	printf("<GeoElevationGrid DEF='MOON' geoGridOrigin='-90.0 -180 0' ");
	printf(" xDimension='%d' xSpacing='%lf' zDimension='%d' zSpacing='%lf'",
		gimp_image.width+1, 360.0 / (double) gimp_image.width,
		gimp_image.height,  180.0 / (double) (gimp_image.height -1));
	printf(" yScale='1.0' height='");
	count = 0;
	for(i=0;i<gimp_image.height;i++){
		ii = gimp_image.height - i -1; //work from bottom up in the image
		for(j=0;j<gimp_image.width;j++){
			for(k=0;k<gimp_image.bytes_per_pixel;k++){
				p = gimp_image.pixel_data[(ii*gimp_image.width+j)*gimp_image.bytes_per_pixel + k];
				if(k==0) {
					h = (float)p + .5f;
					printf("%.1f ",h);
					if(j == 0) h0 = h;
					count++;
				}
			}
		}
		printf("%.1f ",h0);  //cylindrical projectsions have a seam, and we need to do the seam points 2x
		count ++;
	}
	printf("' />\n");
	if(count != (gimp_image.width+1)*gimp_image.height)
		printf("<!-- warning n(%d) != width(%d)xheight(%d)==(%d) -->",count,(gimp_image.width+1),gimp_image.height,(gimp_image.width+1)*gimp_image.height);
	getchar();
	return 0;
}