网站开发 一个页面多少钱,建设部网站查资质中裕隆,怎么做网站做站点,百度网站认证v1clc;close all;clear;
%地理坐标和像素坐标的相互转换
[pic,R]geotiffread(boston.tif); %读取带地理坐标信息的tif影像
[m,n,~]size(pic);%像素坐标转换为地理坐标
figure(1),imshow(pic),title(不带地理坐标的静态图片);
hold on; scatter(n/4,m/4,500,r.); %选择1/4处像素坐…clc;close all;clear;
%地理坐标和像素坐标的相互转换
[pic,R]geotiffread(boston.tif); %读取带地理坐标信息的tif影像
[m,n,~]size(pic);%像素坐标转换为地理坐标
figure(1),imshow(pic),title(不带地理坐标的静态图片);
hold on; scatter(n/4,m/4,500,r.); %选择1/4处像素坐标并在图上标示
[lon,lat]pix2map(R,m/4,n/4); %转换为地理坐标
figure(2),mapshow(pic,R);
mapshow(lon,lat,Marker,.,MarkerEdgeColor,r);
title(带地理坐标的tif影像);
axis off;
%打印结果
disp([(,num2str(m/4),,,num2str(n/4),) - (,num2str(lon),,,num2str(lat),)]); %地理坐标转换为像素坐标。注意地理坐标从左下角开始而像素坐标从左上角开始
%左下角开始点地理坐标为[R.XLimWorld(1),R.YLimWorld(1)]
%R.RasterWidthInWorld表示图像的地理宽度
xR.XLimWorld(1)(3/4)*R.RasterWidthInWorld; %获取图像3/4处的地理坐标
yR.YLimWorld(1)(1/4)*R.RasterHeightInWorld;
figure(3),mapshow(pic,R),axis off;
mapshow(x,y,Marker,*,MarkerEdgeColor,r);
[row,col]map2pix(R,x,y);
figure(4),imshow(pic);
hold on;
scatter(col,row,100,r*);
%命令行中打印结果
disp([(,num2str(x),,,num2str(y),) - (,num2str(row),,,num2str(col),)]); 像素坐标转地理坐标pix2map
地理坐标转像素坐标map2pix
大地坐标与三维直角坐标系的互相转换两个程序