Usage examples

let img_file =BufReader::new(File::open("N265E425.tif").unwrap());letmut tiff =GeoTiffReader::open(img_file).unwrap();match tiff.read_pixel(x, y){RasterValue::U16(v) =>println!("Height: {v}"),    _ =>println!("Unexpected pixel type"),};

let img_file =BufReader::new(File::open("N265E425.tif").unwrap());letmut tiff =GeoTiffReader::open(img_file).unwrap();let(x0, y0, w, h) =(2500,3000,100,100);letmut img =ImageBuffer::new(w, h);for(x, y, pixel)in tiff.pixels(x0, y0, w, h){ifletRasterValue::U16(v) = pixel{        img.put_pixel(x - x0, y - y0, image::Luma([v]));}}img.save("dtm.png").unwrap();

Running the examples

cd datamake

cargo run --example info data/tiff/N265E425.tifcargo run --example pixel data/tiff/N265E425.tif 2550 3050cargo run --example crop data/tiff/N265E425.tif 100x100+2500+3000 dtm.pngcargo run --example img2ascii data/tiff/sat.tifcargo run --example http_dtm