关于gshhs安装的问题

Isakting

以前安装过gshhs，也能够正常使用，后来重做系统后再安装的时候就出现了如下问题

>> m_gshhs_l
警告: Coastline file private/gshhs_l.b not found \n(Have you installed it? See the M_Map User's Guide for details)\n ---Using default coastline
instead
> In mu_coast>get_coasts (line 460)
  In mu_coast (line 101)
  In m_gshhs (line 114)
  In m_gshhs_l (line 28)

题主已经按照要求将gshhs相关文件放到了D:\MATLAB2016\toolbox\m_map\private 中，出现了如上错误
题主找到之前有遇到这种问题的，https://www.ilovematlab.cn/thread-277337-1-1.html 他是将每一个m_gshhs_*.m里面的FILENAME都改了
但是我的m_map中的m_gshhs_*.m并没有FILENAME这一行，而是集合到了m_gshhs.m这个脚本里
所以将m_gshhs.m里面的FILENAME也改成了 FILENAME="D:\MATLAB2016\toolbox\m_map\private\"，但是依然无效
再然后题主还在http://compgroups.net/comp.soft-sys.matlab/problem-with-m_map-and-gshhs/393609这篇帖子中看到
I found a bug in function loading the GSHHS. Recently, the number of arguments in the header of polygon stored in GSHHS files have changed.However, I had some trouble loading GSHHS files.Looking inside the code, I found that the number of elements read in the header of each polygon did not match the GSHHS file.It has changed from 8 to 11ftp://ftp.soest.hawaii.edu/pwessel/gshhs/README.TXTIn file 'mu_coast.m' for function '[A,cnt] = get_gheader(fid)'I had to change the line [A,cnt]=fread(fid,8,'int32');to [A,cnt]=fread(fid,11,'int32'); where 11 is the new number of elements to read in the header.Without this change, nothing was displayed.The maintainer of this toolbox is aware of this bug , and will try to solve it.
就是说在上述错误中的mu_coast.m脚本中的读取多边形的元素个数不匹配（从8变成了11）


% Now, in the calling code I have to use A(2),A(3),A(5-7), A(8), A(9) from original.
[A,cnt]=fread(fid,8,'int32');

这行在mu_coast.m中的720行，然而修改为11后没有什么作用。
我使用的是gshhg-bin-2.3.7和最新的m_map，当然gshhs_1.0和2.0版本这些我也都试过，都用不了
我也不知道到底是哪儿出的问题，能试的方法我都试了，请问各位大佬有什么好的解决办法吗？

二爷名声在外

楼主说：

我的m_map中的m_gshhs_*.m并没有FILENAME这一行

建议楼主把程序放出来，我告诉楼主在哪里

Isakting

二爷名声在外 发表于 2019-6-1 10:55
楼主说：
建议楼主把程序放出来，我告诉楼主在哪里

这个是我的m_gshhs.m 脚本文件，其中我认为需要修改的地方已经加粗标红

1. function h=m_gshhs(resolution,varargin)
   
2. % M_GSHHS Add a coastline to a given map using
   
3. %           the Global Self-consistant Hierarchical High-resolution
   
4. %           Shorelines, Rivers, and Borders
   
5. %
   
6. %         M_GSHHS(RES, (standard line option,...,...) ) draws the coastline
   
7. %         river network, or borders as  simple lines.
   
8. %
   
9. %         M_GSHHS(RES,'patch' ( ,standard patch options,...,...) ) draws the
   
10. %         coastline as a number of patches (rivers and borders are not
    
11. %         arranged so patches can be drawn).
    
12. %
    
13. %         M_GSHHS(RES,'save',FILENAME) saves the extracted coastline data
    
14. %         for the current projection in a file FILENAME. This allows
    
15. %         speedier replotting using M_USERCOAST(FILENAME).
    
16. %   
    
17. %         RES: A one-char string (optionally 2 or 3)
    
18. %         
    
19. %         First char: resolution - one of
    
20. %                      'c'  crude
    
21. %                      'l'  low
    
22. %                      'i'  intermediate
    
23. %                      'h'  high
    
24. %                      'f'  full
    
25. %
    
26. %         Second char: type - one of
    
27. %                      'c' GSHHS coastline (default)
    
28. %                      'b' WDB Border
    
29. %                      'r' WDB River
    
30. %  
    
31. %         Third char - if 2nd char is 'b':
    
32. %                      '1' Country borders
    
33. %                      '2' State/Province and Country borders
    
34. %                    - if 2nd char is 'r': '1','2','3','4'
    
35. %                      add successively more tributaries
    
36. %
    
37. %         (also maintained is this optional format:
    
38. %
    
39. %         RES - selections resolution
    
40. %                  1  or 'crude'       
    
41. %                  2  or 'low'         
    
42. %                  3  or 'intermediate'  
    
43. %                  4  or 'high'        
    
44. %                  5  or 'full         
    
45. %
    
46. %         but please don't use this).
    
47. %
    
48. %         See also M_PROJ, M_GRID, M_COAST, M_GSHHS_L, M_GSHHS_H, M_GSHHS_C
    
49. %         M_USERCOAST   
    
50. 
    
51. % Rich Pawlowicz (rich@ocgy.ubc.ca) 15/June/98
    
52. %
    
53. %
    
54. % This software is provided "as is" without warranty of any kind. But
    
55. % it's mine, so you can't sell it.
    
56. %
    
57. %  16/Dec/2005
    
58. %*********************************************************************
    
59. %  Modified after code provided by Bruce Lipphardt (brucel@udel.edu) to
    
60. %  reduce the hierarchy of M_GSHHS_* routines to a single routine with a
    
61. %  variable resolution input:
    
62. % 20/Jan/2008 - added borders and rivers from gshhs v1.10
    
63. % 4/DEc/11 - isstr to ischar
    
64. % Sep/14 - added hierarchy to borders
    
65. % Aug/18 - fixed error that occurred when called m_gshhs_X with 'save'
    
66. %          option (Thanks to H. Grant for pointing this out).
    
67. 
    
68. 
    
69. % Root of directories where all the gshhs_X.b, wdb_borders-X.b and wdb_rivers_X.b
    
70. % files live
    
71. <font color="#ff0000"><b>FILNAME='private/';</b></font>
    
72. 
    
73. 
    
74. %-------------don't change below here----------------------------
    
75. 
    
76. res_list = {'c','l','i','h','f'};
    
77. typ_list={'c','b','r'};
    
78. typ_names={'gshhs_','wdb_borders_','wdb_rivers_'};
    
79. 
    
80. typ=1;
    
81. flaglim='9';
    
82. 
    
83. if ischar(resolution)
    
84. if length(resolution)>=2
    
85.    typ = find(strcmpi(resolution(2),typ_list));
    
86. end
    
87. if length(resolution)>=3
    
88.    flaglim = resolution(3);
    
89. end  
    
90. resolution =  find(strcmpi(resolution(1),res_list));
    
91. end
    
92. 
    
93. 
    
94. if isempty(resolution) || resolution<1 || resolution> length(res_list)
    
95.   error('**Don''t recognize the specified resolution');
    
96. end
    
97. if isempty(typ) || typ<1 || typ> length(res_list)
    
98.   error('**Don''t recognize the specified type');
    
99. end
    
100.   
     
101. res_char = res_list{resolution} ;
     
102. file     = [FILNAME,sprintf('%s%s.b',typ_names{typ},res_char)] ;
     
103. tag_name = sprintf('m_%s%s',typ_names{typ},res_char) ;
     
104. 
     
105. 
     
106. % Set current projection to geographic
     
107. Currentmap=m_coord('set');
     
108. m_coord('geographic');
     
109. 
     
110. 
     
111. if length(varargin)>1 && strcmp(varargin{1},'save')
     
112.   [ncst,Area,k]=mu_coast(res_char,file);
     
113.   save(varargin{2},'ncst','k','Area');
     
114.   h=varargin{2};   % Error if you call m_gshhs_i with 'save' option - thanks HG, Aug/1/2018
     
115. else
     
116.   h=mu_coast([res_char flaglim],file,varargin{:},'tag',tag_name);
     
117. end
     
118. 
     
119. m_coord(Currentmap.name);
     
120. 
     
121. if nargout==0
     
122.     clear h
     
123. end
     

复制代码

这个是我的m_gshhs_i.m脚本文件

1. function h=m_gshhs_i(varargin)
   
2. % M_GSHHS_I Add a coastline to a given map using the 'intermediate' resolution of
   
3. %           the Global Self-consistant Hierarchical High-resolution
   
4. %           Shorelines.
   
5. %
   
6. %         M_GSHHS_I( (standard line option,...,...) ) draws the coastline
   
7. %         as a simple line.
   
8. %         M_GSHHS_I('patch' ( ,standard patch options,...,...) ) draws the
   
9. %         coastline as a number of patches.
   
10. %
    
11. %         M_GSHHS_I('save',FILENAME) saves the extracted coastline data
    
12. %         for the current projection in a file FILENAME. This allows
    
13. %         speedier replotting using M_USERCOAST(FILENAME).
    
14. %   
    
15. %         See also M_PROJ, M_GRID, M_COAST, M_GSHHS_L, M_GSHHS_H, M_GSHHS_C
    
16. %         M_USERCOAST   
    
17. 
    
18. % Rich Pawlowicz (rich@ocgy.ubc.ca) 15/June/98
    
19. %
    
20. %
    
21. % This software is provided "as is" without warranty of any kind. But
    
22. % it's mine, so you can't sell it.
    
23. 
    
24. % Changes
    
25. %  Nov/2017 - changed this into a stub calling gshhs.m (kept
    
26. %             for backwards compatability)
    
27. 
    
28. h=m_gshhs('ic',varargin{:});
    
29. 
    

复制代码

m_gshhs_*.m 其他分辨率的gshhs脚本也和第二个类似，应该是m_gshhs.m来调用对应分辨率的脚本吧？

二爷名声在外

不好意思现在才有时间回复楼主。
我也将我相应的程序的截图放上，供楼主比对

Isakting

二爷名声在外 发表于 2019-7-8 21:59
不好意思现在才有时间回复楼主。
我也将我相应的程序的截图放上，供楼主比对

感谢！已经解决~！

kingjames

Isakting 发表于 2019-9-24 15:28
感谢！已经解决~！

我也遇到了同样的问题，终于解决了  ：FILNAME='D:\matlab2019a\toolbox\m_map\private/';

鼠国三

您好，请问您的m_gshhs_*.b文件是如何打开的？我的用记事本和matlab打开后，都是乱码

Isakting

鼠国三 发表于 2020-7-13 09:12
您好，请问您的m_gshhs_*.b文件是如何打开的？我的用记事本和matlab打开后，都是乱码

不是.b 文件……是.m 脚本

fx805

请问你怎么解决的？我把m_gshhs.m文件里改了FILNAME='C:\Program Files\MATLAB\R2016a\toolbox\m_map\private\gshhs/'，还是不行啊