bbox = [-170 47;-130 65];           %框定经纬度坐标，只对该范围内的shp断层坐标进行读取，要尽可能的缩小该区域，否则参数太多容易卡掉而且faultname甚至可能都读取不出来
s = shaperead('qfaults.shp', 'BoundingBox',bbox,'Attributes',{'Lon','Lat','faultname'})  
s1=rmfield(s,'Geometry')； %删除掉结构体元素：Geometry
s1=rmfield(s1,'BoundingBox')；%删除掉结构体元素：BoundingBox

 s2 = struct2cell(s1);

%合并经度
Lon = s2{1,1};
for i = 2:1412
    Lon = [Lon,s2{1,i}];
end
Lon = Lon';

%合并纬度
Lat = s2{2,1};
for i = 2:1412
    Lat = [Lat,s2{2,i}];
end
Lat = Lat';

%获得经纬度坐标
coordinate= [Lon,Lat];

gmt psxy AlaskaFault.txt -J$J -R$R -w0.2p,red,solid -K -O >> PS