carray 1.3.6 → 1.5.3

Sign up to get free protection for your applications and to get access to all the features.
Files changed (315) hide show
  1. checksums.yaml +4 -4
  2. data/.yardopts +32 -0
  3. data/LICENSE +22 -0
  4. data/NEWS.md +65 -0
  5. data/README.md +23 -26
  6. data/Rakefile +31 -0
  7. data/TODO.md +17 -0
  8. data/carray.gemspec +24 -31
  9. data/{ca_iter_block.c → ext/ca_iter_block.c} +16 -18
  10. data/{ca_iter_dimension.c → ext/ca_iter_dimension.c} +20 -21
  11. data/{ca_iter_window.c → ext/ca_iter_window.c} +12 -14
  12. data/{ca_obj_array.c → ext/ca_obj_array.c} +451 -176
  13. data/{ca_obj_bitarray.c → ext/ca_obj_bitarray.c} +18 -23
  14. data/{ca_obj_bitfield.c → ext/ca_obj_bitfield.c} +12 -16
  15. data/{ca_obj_block.c → ext/ca_obj_block.c} +47 -54
  16. data/{ca_obj_fake.c → ext/ca_obj_fake.c} +10 -12
  17. data/{ca_obj_farray.c → ext/ca_obj_farray.c} +21 -23
  18. data/{ca_obj_field.c → ext/ca_obj_field.c} +30 -32
  19. data/{ca_obj_grid.c → ext/ca_obj_grid.c} +32 -33
  20. data/{ca_obj_mapping.c → ext/ca_obj_mapping.c} +11 -13
  21. data/{ca_obj_object.c → ext/ca_obj_object.c} +40 -42
  22. data/{ca_obj_reduce.c → ext/ca_obj_reduce.c} +3 -5
  23. data/{ca_obj_refer.c → ext/ca_obj_refer.c} +38 -40
  24. data/{ca_obj_repeat.c → ext/ca_obj_repeat.c} +45 -47
  25. data/{ca_obj_select.c → ext/ca_obj_select.c} +4 -6
  26. data/{ca_obj_shift.c → ext/ca_obj_shift.c} +26 -28
  27. data/{ca_obj_transpose.c → ext/ca_obj_transpose.c} +26 -28
  28. data/{ca_obj_unbound_repeat.c → ext/ca_obj_unbound_repeat.c} +106 -160
  29. data/{ca_obj_window.c → ext/ca_obj_window.c} +33 -35
  30. data/{carray.h → ext/carray.h} +90 -59
  31. data/{carray_access.c → ext/carray_access.c} +244 -109
  32. data/{carray_attribute.c → ext/carray_attribute.c} +161 -207
  33. data/{carray_call_cfunc.c → ext/carray_call_cfunc.c} +1 -3
  34. data/{carray_cast.c → ext/carray_cast.c} +350 -149
  35. data/{carray_cast_func.rb → ext/carray_cast_func.rb} +1 -2
  36. data/{carray_class.c → ext/carray_class.c} +28 -36
  37. data/{carray_conversion.c → ext/carray_conversion.c} +58 -70
  38. data/{carray_copy.c → ext/carray_copy.c} +34 -50
  39. data/{carray_core.c → ext/carray_core.c} +75 -62
  40. data/ext/carray_data_type.c +66 -0
  41. data/{carray_element.c → ext/carray_element.c} +34 -53
  42. data/{carray_generate.c → ext/carray_generate.c} +137 -50
  43. data/{carray_iterator.c → ext/carray_iterator.c} +53 -53
  44. data/{carray_loop.c → ext/carray_loop.c} +77 -106
  45. data/{carray_mask.c → ext/carray_mask.c} +95 -114
  46. data/{carray_math.rb → ext/carray_math.rb} +47 -31
  47. data/ext/{mathfunc/carray_mathfunc.c → carray_mathfunc.c} +1 -3
  48. data/{carray_numeric.c → ext/carray_numeric.c} +43 -46
  49. data/{carray_operator.c → ext/carray_operator.c} +72 -36
  50. data/{carray_order.c → ext/carray_order.c} +73 -217
  51. data/{carray_sort_addr.c → ext/carray_sort_addr.c} +14 -21
  52. data/{carray_stat.c → ext/carray_stat.c} +6 -8
  53. data/{carray_stat_proc.rb → ext/carray_stat_proc.rb} +25 -27
  54. data/{carray_test.c → ext/carray_test.c} +59 -51
  55. data/{carray_undef.c → ext/carray_undef.c} +1 -11
  56. data/{carray_utils.c → ext/carray_utils.c} +12 -4
  57. data/{extconf.rb → ext/extconf.rb} +1 -6
  58. data/{mkmath.rb → ext/mkmath.rb} +13 -3
  59. data/{ruby_carray.c → ext/ruby_carray.c} +20 -8
  60. data/{ruby_ccomplex.c → ext/ruby_ccomplex.c} +1 -3
  61. data/{ruby_float_func.c → ext/ruby_float_func.c} +1 -3
  62. data/ext/version.h +16 -0
  63. data/{version.rb → ext/version.rb} +0 -0
  64. data/lib/carray.rb +52 -10
  65. data/lib/carray/{base/autoload.rb → autoload.rb} +48 -8
  66. data/lib/carray/autoload/autoload_base.rb +1 -5
  67. data/lib/carray/autoload/autoload_gem_cairo.rb +9 -0
  68. data/lib/carray/autoload/autoload_gem_ffi.rb +9 -0
  69. data/lib/carray/autoload/autoload_gem_gnuplot.rb +2 -0
  70. data/lib/carray/autoload/autoload_gem_io_csv.rb +14 -0
  71. data/lib/carray/autoload/autoload_gem_io_pg.rb +6 -0
  72. data/lib/carray/autoload/autoload_gem_io_sqlite3.rb +12 -0
  73. data/lib/carray/autoload/autoload_gem_narray.rb +10 -0
  74. data/lib/carray/autoload/autoload_gem_numo_narray.rb +15 -0
  75. data/lib/carray/autoload/autoload_gem_opencv.rb +16 -0
  76. data/lib/carray/autoload/autoload_gem_random.rb +8 -0
  77. data/lib/carray/autoload/autoload_gem_rmagick.rb +23 -0
  78. data/lib/carray/autoload/{autoload_graphics_zimg.rb → autoload_gem_zimg.rb} +0 -0
  79. data/lib/carray/basic.rb +191 -0
  80. data/lib/carray/broadcast.rb +82 -0
  81. data/lib/carray/compose.rb +315 -0
  82. data/lib/carray/construct.rb +460 -0
  83. data/lib/carray/convert.rb +115 -0
  84. data/lib/carray/info.rb +1 -3
  85. data/lib/carray/{base/inspect.rb → inspect.rb} +9 -11
  86. data/lib/carray/io/imagemagick.rb +2 -4
  87. data/lib/carray/{base/iterator.rb → iterator.rb} +6 -6
  88. data/lib/carray/mask.rb +102 -0
  89. data/lib/carray/{base/math.rb → math.rb} +21 -53
  90. data/lib/carray/math/histogram.rb +8 -10
  91. data/lib/carray/math/recurrence.rb +1 -3
  92. data/lib/carray/mkmf.rb +9 -3
  93. data/lib/carray/object/ca_obj_iterator.rb +1 -3
  94. data/lib/carray/object/ca_obj_link.rb +1 -3
  95. data/lib/carray/object/ca_obj_pack.rb +9 -11
  96. data/lib/carray/obsolete.rb +270 -0
  97. data/lib/carray/ordering.rb +181 -0
  98. data/lib/carray/{base/serialize.rb → serialize.rb} +50 -66
  99. data/lib/carray/string.rb +188 -0
  100. data/lib/carray/{base/struct.rb → struct.rb} +19 -21
  101. data/lib/carray/{io/table.rb → table.rb} +1 -9
  102. data/lib/carray/testing.rb +51 -0
  103. data/lib/carray/time.rb +76 -0
  104. data/lib/carray/transform.rb +98 -0
  105. data/misc/Methods.ja.md +182 -0
  106. data/{NOTE → misc/NOTE} +0 -0
  107. data/spec/Classes/CABitfield_spec.rb +58 -0
  108. data/spec/Classes/CABlockIterator_spec.rb +114 -0
  109. data/spec/Classes/CABlock_spec.rb +205 -0
  110. data/spec/Classes/CAField_spec.rb +39 -0
  111. data/spec/Classes/CAGrid_spec.rb +75 -0
  112. data/spec/Classes/CAMap_spec.rb +0 -0
  113. data/{test/test_CAMapping.rb → spec/Classes/CAMapping_spec.rb} +35 -36
  114. data/spec/Classes/CAObject_attribute_spec.rb +33 -0
  115. data/spec/Classes/CAObject_spec.rb +33 -0
  116. data/spec/Classes/CARefer_spec.rb +93 -0
  117. data/spec/Classes/CARepeat_spec.rb +65 -0
  118. data/spec/Classes/CASelect_spec.rb +22 -0
  119. data/spec/Classes/CAShift_spec.rb +16 -0
  120. data/spec/Classes/CAStruct_spec.rb +71 -0
  121. data/{test/test_CATranspose.rb → spec/Classes/CATranspose_spec.rb} +20 -21
  122. data/spec/Classes/CAUnboudRepeat_spec.rb +102 -0
  123. data/spec/Classes/CAWindow_spec.rb +54 -0
  124. data/spec/Classes/CAWrap_spec.rb +8 -0
  125. data/{test/test_CArray.rb → spec/Classes/CArray_spec.rb} +48 -92
  126. data/spec/Classes/CScalar_spec.rb +55 -0
  127. data/spec/Features/feature_130_spec.rb +19 -0
  128. data/spec/Features/feature_attributes_spec.rb +280 -0
  129. data/spec/Features/feature_boolean_spec.rb +98 -0
  130. data/spec/Features/feature_broadcast.rb +116 -0
  131. data/spec/Features/feature_cast_function.rb +19 -0
  132. data/spec/Features/feature_cast_spec.rb +33 -0
  133. data/spec/Features/feature_class_spec.rb +84 -0
  134. data/spec/Features/feature_complex_spec.rb +42 -0
  135. data/{test/test_composite.rb → spec/Features/feature_composite_spec.rb} +17 -18
  136. data/spec/Features/feature_convert_spec.rb +46 -0
  137. data/spec/Features/feature_copy_spec.rb +123 -0
  138. data/spec/Features/feature_creation_spec.rb +84 -0
  139. data/spec/Features/feature_element_spec.rb +144 -0
  140. data/spec/Features/feature_extream_spec.rb +54 -0
  141. data/spec/Features/feature_generate_spec.rb +74 -0
  142. data/spec/Features/feature_index_spec.rb +69 -0
  143. data/spec/Features/feature_mask_spec.rb +580 -0
  144. data/spec/Features/feature_math_spec.rb +97 -0
  145. data/spec/Features/feature_order_spec.rb +146 -0
  146. data/spec/Features/feature_ref_store_spec.rb +209 -0
  147. data/spec/Features/feature_serialization_spec.rb +125 -0
  148. data/spec/Features/feature_stat_spec.rb +397 -0
  149. data/spec/Features/feature_virtual_spec.rb +48 -0
  150. data/spec/Features/method_eq_spec.rb +81 -0
  151. data/spec/Features/method_is_nan_spec.rb +12 -0
  152. data/spec/Features/method_map_spec.rb +54 -0
  153. data/spec/Features/method_max_with.rb +20 -0
  154. data/spec/Features/method_min_with.rb +19 -0
  155. data/spec/Features/method_ne_spec.rb +18 -0
  156. data/spec/Features/method_project_spec.rb +188 -0
  157. data/spec/Features/method_ref_spec.rb +27 -0
  158. data/spec/Features/method_round_spec.rb +11 -0
  159. data/spec/Features/method_s_linspace_spec.rb +48 -0
  160. data/spec/Features/method_s_span_spec.rb +14 -0
  161. data/spec/Features/method_seq_spec.rb +47 -0
  162. data/spec/Features/method_sort_with.rb +43 -0
  163. data/spec/Features/method_sorted_with.rb +29 -0
  164. data/spec/Features/method_span_spec.rb +42 -0
  165. data/spec/Features/method_wrap_readonly_spec.rb +43 -0
  166. data/{test → spec/UnitTest}/test_CAVirtual.rb +0 -0
  167. data/spec/spec_all.rb +0 -1
  168. data/utils/convert_test.rb +73 -0
  169. data/utils/{extract_rdoc.rb → extract_yard.rb} +7 -12
  170. data/{devel → utils}/guess_shape.rb +0 -0
  171. data/utils/{diff_method.rb → monkey_patch_methods.rb} +17 -7
  172. metadata +173 -316
  173. data/COPYING +0 -56
  174. data/GPL +0 -340
  175. data/Gemfile +0 -8
  176. data/Gemfile.lock +0 -33
  177. data/LEGAL +0 -50
  178. data/Makefile +0 -451
  179. data/TODO +0 -5
  180. data/carray_cast_func.c +0 -8670
  181. data/carray_config.h +0 -26
  182. data/carray_math.c +0 -24805
  183. data/carray_random.c +0 -531
  184. data/carray_stat_proc.c +0 -12170
  185. data/ext/calculus/carray_calculus.c +0 -931
  186. data/ext/calculus/carray_interp.c +0 -358
  187. data/ext/calculus/extconf.rb +0 -12
  188. data/ext/calculus/lib/autoload/autoload_math_calculus.rb +0 -2
  189. data/ext/calculus/lib/math/calculus.rb +0 -119
  190. data/ext/calculus/lib/math/interp/adapter_interp1d.rb +0 -31
  191. data/ext/fortio/extconf.rb +0 -3
  192. data/ext/fortio/lib/carray/autoload/autoload_fortran_format.rb +0 -5
  193. data/ext/fortio/lib/carray/io/fortran_format.rb +0 -43
  194. data/ext/fortio/lib/fortio.rb +0 -3
  195. data/ext/fortio/lib/fortio/fortran_format.rb +0 -605
  196. data/ext/fortio/lib/fortio/fortran_format.tab.rb +0 -536
  197. data/ext/fortio/lib/fortio/fortran_format.y +0 -215
  198. data/ext/fortio/lib/fortio/fortran_namelist.rb +0 -151
  199. data/ext/fortio/lib/fortio/fortran_namelist.tab.rb +0 -470
  200. data/ext/fortio/lib/fortio/fortran_namelist.y +0 -213
  201. data/ext/fortio/lib/fortio/fortran_sequential.rb +0 -345
  202. data/ext/fortio/ruby_fortio.c +0 -182
  203. data/ext/fortio/test/test_H.rb +0 -5
  204. data/ext/fortio/test/test_T.rb +0 -7
  205. data/ext/fortio/test/test_fortran_format.rb +0 -86
  206. data/ext/fortio/test/test_namelist.rb +0 -25
  207. data/ext/fortio/test/test_namelist_write.rb +0 -10
  208. data/ext/fortio/test/test_sequential.rb +0 -13
  209. data/ext/fortio/test/test_sequential2.rb +0 -13
  210. data/ext/fortio/work/test.rb +0 -10
  211. data/ext/fortio/work/test_e.rb +0 -19
  212. data/ext/fortio/work/test_ep.rb +0 -10
  213. data/ext/fortio/work/test_parse.rb +0 -12
  214. data/ext/imagemap/carray_imagemap.c +0 -495
  215. data/ext/imagemap/doc/call_graph.dot +0 -28
  216. data/ext/imagemap/draw.c +0 -567
  217. data/ext/imagemap/extconf.rb +0 -13
  218. data/ext/imagemap/lib/autoload/autoload_graphics_imagemap.rb +0 -1
  219. data/ext/imagemap/lib/graphics/imagemap.rb +0 -273
  220. data/ext/imagemap/lib/image_map.rb +0 -4
  221. data/ext/imagemap/test/swath_index.rb +0 -83
  222. data/ext/imagemap/test/swath_warp.rb +0 -99
  223. data/ext/imagemap/test/test.rb +0 -23
  224. data/ext/imagemap/test/test_image.rb +0 -42
  225. data/ext/imagemap/test/test_line.rb +0 -14
  226. data/ext/imagemap/test/test_rotate.rb +0 -17
  227. data/ext/imagemap/test/test_triangle.rb +0 -20
  228. data/ext/imagemap/test/test_warp.rb +0 -26
  229. data/ext/mathfunc/extconf.rb +0 -18
  230. data/ext/mathfunc/lib/autoload/autoload_math_mathfunc.rb +0 -1
  231. data/ext/mathfunc/lib/math/mathfunc.rb +0 -15
  232. data/ext/mathfunc/test/test_hypot.rb +0 -5
  233. data/ext/mathfunc/test/test_j0.rb +0 -22
  234. data/ext/mathfunc/test/test_jn.rb +0 -8
  235. data/ext/mathfunc/test/test_sph.rb +0 -9
  236. data/ext/narray/README +0 -22
  237. data/ext/narray/ca_wrap_narray.c +0 -500
  238. data/ext/narray/carray_narray.c +0 -21
  239. data/ext/narray/extconf.rb +0 -57
  240. data/ext/narray/lib/autoload/autoload_math_narray.rb +0 -1
  241. data/ext/narray/lib/autoload/autoload_math_narray_miss.rb +0 -11
  242. data/ext/narray/lib/math/narray.rb +0 -17
  243. data/ext/narray/lib/math/narray_miss.rb +0 -45
  244. data/lib/carray/autoload/autoload_graphics_gnuplot.rb +0 -2
  245. data/lib/carray/autoload/autoload_io_csv.rb +0 -14
  246. data/lib/carray/autoload/autoload_io_pg.rb +0 -6
  247. data/lib/carray/autoload/autoload_io_sqlite3.rb +0 -12
  248. data/lib/carray/autoload/autoload_io_table.rb +0 -1
  249. data/lib/carray/autoload/autoload_math_interp.rb +0 -4
  250. data/lib/carray/base/basic.rb +0 -1090
  251. data/lib/carray/base/obsolete.rb +0 -131
  252. data/lib/carray/graphics/gnuplot.rb +0 -2131
  253. data/lib/carray/graphics/zimg.rb +0 -296
  254. data/lib/carray/io/csv.rb +0 -572
  255. data/lib/carray/io/pg.rb +0 -101
  256. data/lib/carray/io/sqlite3.rb +0 -215
  257. data/lib/carray/math/interp.rb +0 -57
  258. data/lib/carray/math/interp/adapter_gsl_spline.rb +0 -47
  259. data/mkmf.log +0 -375
  260. data/mt19937ar.c +0 -182
  261. data/mt19937ar.h +0 -86
  262. data/rdoc_ext.rb +0 -1115
  263. data/rdoc_main.rb +0 -27
  264. data/rdoc_math.rb +0 -5
  265. data/rdoc_stat.rb +0 -31
  266. data/spec/CABlockIterator/CABlockIterator_spec.rb +0 -113
  267. data/spec/CArray/bug/store_spec.rb +0 -27
  268. data/spec/CArray/index/repeat_spec.rb +0 -10
  269. data/spec/CArray/method/eq_spec.rb +0 -80
  270. data/spec/CArray/method/is_nan_spec.rb +0 -12
  271. data/spec/CArray/method/ne_spec.rb +0 -18
  272. data/spec/CArray/method/round_spec.rb +0 -11
  273. data/spec/CArray/object/_attribute_spec.rb +0 -32
  274. data/spec/CArray/object/s_new_spec.rb +0 -31
  275. data/spec/CArray/serialize/Serialization_spec.rb +0 -89
  276. data/test/test_130.rb +0 -23
  277. data/test/test_ALL.rb +0 -51
  278. data/test/test_CABitfield.rb +0 -59
  279. data/test/test_CABlock.rb +0 -208
  280. data/test/test_CAField.rb +0 -40
  281. data/test/test_CAGrid.rb +0 -76
  282. data/test/test_CAMmap.rb +0 -11
  283. data/test/test_CARefer.rb +0 -94
  284. data/test/test_CARepeat.rb +0 -66
  285. data/test/test_CASelect.rb +0 -23
  286. data/test/test_CAShift.rb +0 -17
  287. data/test/test_CAWindow.rb +0 -55
  288. data/test/test_CAWrap.rb +0 -9
  289. data/test/test_CComplex.rb +0 -83
  290. data/test/test_CScalar.rb +0 -91
  291. data/test/test_attribute.rb +0 -281
  292. data/test/test_block_iterator.rb +0 -17
  293. data/test/test_boolean.rb +0 -99
  294. data/test/test_cast.rb +0 -33
  295. data/test/test_class.rb +0 -85
  296. data/test/test_complex.rb +0 -43
  297. data/test/test_convert.rb +0 -79
  298. data/test/test_copy.rb +0 -141
  299. data/test/test_creation.rb +0 -85
  300. data/test/test_element.rb +0 -146
  301. data/test/test_extream.rb +0 -55
  302. data/test/test_generate.rb +0 -75
  303. data/test/test_index.rb +0 -71
  304. data/test/test_mask.rb +0 -578
  305. data/test/test_math.rb +0 -98
  306. data/test/test_narray.rb +0 -64
  307. data/test/test_order.rb +0 -147
  308. data/test/test_random.rb +0 -15
  309. data/test/test_ref_store.rb +0 -211
  310. data/test/test_stat.rb +0 -414
  311. data/test/test_struct.rb +0 -72
  312. data/test/test_virtual.rb +0 -49
  313. data/utils/create_rdoc.sh +0 -9
  314. data/utils/make_tgz.sh +0 -3
  315. data/version.h +0 -18
@@ -1,13 +0,0 @@
1
- require 'mkmf'
2
-
3
- dir_config("carray", "../..", "../..")
4
-
5
- have_header("carray.h")
6
-
7
- if /cygwin|mingw/ =~ RUBY_PLATFORM
8
- have_library("carray")
9
- end
10
-
11
- create_makefile("carray/carray_imagemap")
12
-
13
-
@@ -1 +0,0 @@
1
- autoload :ImageMap, 'carray/graphics/imagemap'
@@ -1,273 +0,0 @@
1
- # ----------------------------------------------------------------------------
2
- #
3
- # carray/graphics/imagemap.rb
4
- #
5
- # This file is part of Ruby/CArray extension library.
6
- # You can redistribute it and/or modify it under the terms of
7
- # the Ruby Licence.
8
- #
9
- # Copyright (C) 2005 Hiroki Motoyoshi
10
- #
11
- # ----------------------------------------------------------------------------
12
-
13
- require 'carray'
14
-
15
- #
16
- # CArray#index_linear(scale)
17
- # CArray#index_log10(scale)
18
- #
19
- # ImageMap.create_image_map_linear(idx, n, palette, mask_color)
20
- # ImageMap.create_image_map_log10(idx, n, palette, mask_color)
21
- #
22
-
23
- class ImageMap < CArray
24
-
25
- def line (x0, y0, x1, y1, z)
26
- x0 = (x0 - @xoffset)*@xscale
27
- y0 = (y0 - @yoffset)*@yscale
28
- x1 = (x1 - @xoffset)*@xscale
29
- y1 = (y1 - @yoffset)*@yscale
30
- draw_line(y0, x0, y1, x1, z)
31
- end
32
-
33
- def polyline (x, y, z)
34
- x = (x - @xoffset).mul!(@xscale)
35
- y = (y - @xoffset).mul!(@xscale)
36
- draw_polyline(y, x, z)
37
- end
38
-
39
- def map_points (pw, pu, pv)
40
- lines, pixels = pw.dim0, pw.dim1
41
- gx = CArray.float(lines, pixels)
42
- gy = CArray.float(lines, pixels)
43
- yield(pu,pv,gx,gy)
44
- gx.sub!(@xoffset).mul!(@xscale)
45
- gy.sub!(@yoffset).mul!(@yscale)
46
- draw_points(gy, gx, pw)
47
- end
48
-
49
- def warp_image0 (pw, pu, pv)
50
- lines, pixels = pw.dim0, pw.dim1
51
- mu = CArray.float(lines, pixels+1)
52
- mv = CArray.float(lines, pixels+1)
53
- gu = CArray.float(lines+1, pixels+1)
54
- gv = CArray.float(lines+1, pixels+1)
55
- gx = CArray.float(lines+1, pixels+1)
56
- gy = CArray.float(lines+1, pixels+1)
57
- CArray.attach(pu,pv,pw) {
58
- [[pu,mu,gu], [pv,mv,gv]].each do |p,m,g|
59
- m[nil,1..-2] = 0.5*p[nil,0..-2] + 0.5*p[nil,1..-1]
60
- m[nil,0] = 1.5*p[nil,0] - 0.5*p[nil,1]
61
- m[nil,-1] = -0.5*p[nil,-2] + 1.5*p[nil,-1]
62
- g[1..-2,nil] = 0.5*m[0..-2,nil] + 0.5*m[1..-1,nil]
63
- g[ 0,nil] = 1.5*m[0,nil] - 0.5*m[1,nil]
64
- g[-1,nil] = -0.5*m[-2,nil] + 1.5*m[-1,nil]
65
- end
66
- yield(gu,gv,gx,gy)
67
- gx.sub!(@xoffset).mul!(@xscale)
68
- gy.sub!(@yoffset).mul!(@yscale)
69
- fill_rectangle_image(gy, gx, pw)
70
- }
71
- end
72
-
73
- def warp_image1 (pw, pu, pv)
74
- lines, pixels = pw.dim0, pw.dim1
75
- px = CArray.float(lines, pixels)
76
- py = CArray.float(lines, pixels)
77
- mx = CArray.float(lines, pixels+1)
78
- my = CArray.float(lines, pixels+1)
79
- gx = CArray.float(lines+1, pixels+1)
80
- gy = CArray.float(lines+1, pixels+1)
81
- CArray.attach(pu,pv,pw) {
82
- yield(pu,pv,px,py)
83
- px.sub!(@xoffset).mul!(@xscale)
84
- py.sub!(@yoffset).mul!(@yscale)
85
- [[px,mx,gx], [py,my,gy]].each do |p,m,g|
86
- m[nil,1..-2] = 0.5*p[nil,0..-2] + 0.5*p[nil,1..-1]
87
- m[nil,0] = 1.5*p[nil,0] - 0.5*p[nil,1]
88
- m[nil,-1] = -0.5*p[nil,-2] + 1.5*p[nil,-1]
89
- g[1..-2,nil] = 0.5*m[0..-2,nil] + 0.5*m[1..-1,nil]
90
- g[ 0,nil] = 1.5*m[0,nil] - 0.5*m[1,nil]
91
- g[-1,nil] = -0.5*m[-2,nil] + 1.5*m[-1,nil]
92
- end
93
- fill_rectangle_image(gy, gx, pw)
94
- }
95
- end
96
-
97
- def warp_grid (pw, pu, pv)
98
- lines, pixels = pw.dim0, pw.dim1
99
- gx = CArray.float(lines, pixels)
100
- gy = CArray.float(lines, pixels)
101
- CArray.attach(pu,pv,pw) {
102
- yield(pu,pv,gx,gy)
103
- gx.sub!(@xoffset).mul!(@xscale)
104
- gy.sub!(@yoffset).mul!(@yscale)
105
- fill_rectangle_grid(gy, gx, pw)
106
- }
107
- end
108
-
109
- def warp_image_gradation (pw, pu, pv)
110
- lines, pixels = pw.dim0, pw.dim1
111
- mu = CArray.float(lines, pixels+1)
112
- mv = CArray.float(lines, pixels+1)
113
- mw = CArray.float(lines, pixels+1)
114
- gu = CArray.float(lines+1, pixels+1)
115
- gv = CArray.float(lines+1, pixels+1)
116
- gw = CArray.float(lines+1, pixels+1)
117
- gx = CArray.float(lines+1, pixels+1)
118
- gy = CArray.float(lines+1, pixels+1)
119
- # bgx = gx[[0,2],[0,2]]
120
- # bgy = gy[[0,2],[0,2]]
121
- # bgw = gw[[0,2],[0,2]]
122
- # p("hello")
123
- CArray.attach(pu,pv,pw) {
124
- [[pu,mu,gu], [pv,mv,gv], [pw.float32,mw,gw]].each do |p,m,g|
125
- m[nil,1..-2] = 0.5*p[nil,0..-2] + 0.5*p[nil,1..-1]
126
- m[nil,0] = 1.5*p[nil,0] - 0.5*p[nil,1]
127
- m[nil,-1] = -0.5*p[nil,-2] + 1.5*p[nil,-1]
128
- g[1..-2,nil] = 0.5*m[0..-2,nil] + 0.5*m[1..-1,nil]
129
- g[ 0,nil] = 1.5*m[0,nil] - 0.5*m[1,nil]
130
- g[-1,nil] = -0.5*m[-2,nil] + 1.5*m[-1,nil]
131
- end
132
- yield(gu,gv,gx,gy)
133
- gx.sub!(@xoffset).mul!(@xscale)
134
- gy.sub!(@yoffset).mul!(@yscale)
135
- draw_rectangle_gradation_grid(gy, gx, gw)
136
- # pw.each_index do |i,j|
137
- # draw_rectangle_gradation(bgy.move(i,j), bgx.move(i,j), bgw.move(i,j))
138
- # end
139
- }
140
- end
141
-
142
- def warp_grid_gradation (pw, pu, pv)
143
- lines, pixels = pu.dim0, pv.dim1
144
- gx = CArray.float(lines, pixels)
145
- gy = CArray.float(lines, pixels)
146
- CArray.attach(pu,pv,pw) {
147
- yield(pu,pv,gx,gy)
148
- gx.sub!(@xoffset).mul!(@xscale)
149
- gy.sub!(@yoffset).mul!(@yscale)
150
- # bx = gx[[0,2],[0,2]]
151
- # by = gy[[0,2],[0,2]]
152
- # bw = pw[[0,2],[0,2]]
153
- draw_rectangle_gradation_grid(gy, gx, pw)
154
- # pw[0..-2,0..-2].each_index do |i,j|
155
- # draw_rectangle_gradation(by.move(i,j), bx.move(i,j), bw.move(i,j))
156
- # end
157
- }
158
- end
159
-
160
- def warp_grid_interpolation (pw, pu, pv, n, &block)
161
- lines, pixels = pu.dim0, pv.dim1
162
- m = 2+n
163
- hu = CArray.float(m, m)
164
- hv = CArray.float(m, m)
165
- hw = CArray.float(m, m)
166
- bu = pu[[0,2],[0,2]]
167
- bv = pv[[0,2],[0,2]]
168
- bw = pw[[0,2],[0,2]]
169
- CArray.attach(pu,pv,pw) {
170
- pw[0..-2,0..-2].each_index do |i,j|
171
- bu.move(i,j)
172
- bv.move(i,j)
173
- bw.move(i,j)
174
- hu.grid_new([[0,m-1],[0,m-1]])[] = bu
175
- hv.grid_new([[0,m-1],[0,m-1]])[] = bv
176
- hw.grid_new([[0,m-1],[0,m-1]])[] = bw
177
- [hu, hv, hw].each do |h|
178
- h[nil, 0].span!(h[0,0]..h[m-1,0])
179
- h[nil, m-1].span!(h[0,m-1]..h[m-1,m-1])
180
- h.dim0.times do |i|
181
- h[i, nil].span!(h[i,0]..h[i,m-1])
182
- end
183
- end
184
- warp_grid_gradation(hw, hu, hv, &block)
185
- end
186
- }
187
- end
188
-
189
- def initialize (*argv)
190
- super
191
- @yscale = 1
192
- @xscale = 1
193
- @yoffset = 0
194
- @xoffset = 0
195
- end
196
-
197
- attr_accessor :yscale, :xscale, :yoffset, :xoffset
198
-
199
- def set_xrange (rng)
200
- first = rng.begin
201
- last = rng.end
202
- if rng.exclude_end?
203
- @xscale = (dim1-1).to_f/(last-first).to_f
204
- else
205
- @xscale = dim1.to_f/(last-first).to_f
206
- end
207
- @xoffset = first
208
- end
209
-
210
- def set_yrange (rng)
211
- first = rng.begin
212
- last = rng.end
213
- if rng.exclude_end?
214
- @yscale = (dim0-1).to_f/(first-last).to_f
215
- @yoffset = last - @yscale
216
- else
217
- @yscale = dim0.to_f/(first-last).to_f
218
- @yoffset = last
219
- end
220
- end
221
-
222
- def warp (data, x, y, opt = {}, &block)
223
- option = {
224
- :grid => "point",
225
- :interpolation => nil,
226
- :gradation => false,
227
- }
228
- option.update(opt)
229
- case option[:grid]
230
- when /^area([01]|)$/, "pixel"
231
- if option[:interpolation]
232
- raise("can't interpolate in pixel mode")
233
- else
234
- if option[:gradation]
235
- warp_image_gradation(data, x, y, &block)
236
- else
237
- if option[:grid] == "area" or option[:grid] == "area0"
238
- warp_image0(data, x, y, &block) ### inter(extra)polation before
239
- else
240
- warp_image1(data, x, y, &block) ### inter(extra)polation after
241
- end
242
- end
243
- end
244
- when "point"
245
- if option[:interpolation]
246
- warp_grid_interpolation(data, x, y, option[:interpolation], &block)
247
- else
248
- if option[:gradation]
249
- warp_grid_gradation(data, x, y, &block)
250
- else
251
- warp_grid(data, x, y, &block)
252
- end
253
- end
254
- else
255
- raise("unknown grid type")
256
- end
257
- end
258
-
259
- def tfw (metric = 1.0)
260
- format([["%f"]*6].join("\n"),
261
- metric/@xscale,
262
- 0,
263
- 0,
264
- metric/@yscale,
265
- metric*(@xoffset+0.5/@xscale),
266
- metric*(@yoffset+0.5/@yscale))
267
- end
268
-
269
- end
270
-
271
- require 'carray/carray_imagemap.so'
272
-
273
-
@@ -1,4 +0,0 @@
1
- require 'carray/graphics/imagemap'
2
-
3
- warn "carray/image_map.rb will be obsolete"
4
-
@@ -1,83 +0,0 @@
1
- require 'carray'
2
- require 'carray/image_map'
3
- require 'sdi/hdf'
4
- require 'proj4'
5
-
6
- RESOLUTION = 5 # (km) on map projection
7
- UL_CORNER = [3000, 0] # (km) on map projection
8
- LR_CORNER = [0, 3000] # (km) on map projection
9
- SCANLINES = 2
10
-
11
- PROJECTION = "+proj=stere +lat_0=-90 +units=km"
12
-
13
- IMAGE_HEIGHT = (UL_CORNER[0] - LR_CORNER[0])/RESOLUTION.to_f + 1
14
- IMAGE_WIDTH = (LR_CORNER[1] - UL_CORNER[1])/RESOLUTION.to_f + 1
15
-
16
- SCALE_Y = -(UL_CORNER[0] - LR_CORNER[0])/(IMAGE_HEIGHT.to_f - 1)
17
- SCALE_X = (LR_CORNER[1] - UL_CORNER[1])/(IMAGE_WIDTH.to_f - 1)
18
- OFFSET_Y = UL_CORNER[0].to_f # LOWERRIGHT[0].to_f
19
- OFFSET_X = UL_CORNER[1].to_f
20
-
21
- image = ImageMap::Image.new(CA_INT32, [IMAGE_HEIGHT, IMAGE_WIDTH]) { -9999 }
22
- image.set_xrange(UL_CORNER[1]..LR_CORNER[1])
23
- image.set_yrange(LR_CORNER[0]..UL_CORNER[0])
24
-
25
- #
26
- #
27
- #
28
-
29
- sdi = SDI::HDF.open(ARGV[0])
30
- clat = sdi.open("Latitude")[]
31
- clon = sdi.open("Longitude")[]
32
- height = sdi.open("SolarZenith")[]
33
- sdi.close
34
-
35
- NSCANS = clat.dim0/SCANLINES
36
- NPIXELS = clat.dim1
37
-
38
- #cidx = clat.index
39
- #cidx = clat.template.random!(1000)
40
- cidx = height.copy
41
-
42
- proj = Proj4.new(PROJECTION)
43
-
44
- (NSCANS).times do |ns|
45
-
46
- printf("%i/%i\n", ns, NSCANS)
47
-
48
- line = ns * SCANLINES
49
-
50
- bidx = cidx[[line, SCANLINES], nil].copy
51
- blat = clat[[line, SCANLINES], nil].copy
52
- blon = clon[[line, SCANLINES], nil].copy
53
-
54
- image.warp(bidx, blon, blat, :grid => "area", :gradation => false) {
55
- # image.map_points(bidx, blon, blat) {
56
- |lo,la,gx,gy|
57
- proj.fwd_f(lo, la, gx, gy)
58
- }
59
-
60
- end
61
-
62
- puts "projection done"
63
-
64
- #image = height.map_object(image, -9999)
65
- image.mask = image < -9998
66
-
67
- puts "masking flipping done"
68
-
69
- sc = nil
70
- #sc = ImageMap::ColorScale.GMT("GMT_globe.cpt")
71
- #sc = ImageMap::ColorScale.GMT("DEM_poster.cpt")
72
-
73
- ImageMap.display(image, sc)
74
-
75
- open("projection.tfw", "w") { |io|
76
- io.printf([["%g\n"]*6].join(),
77
- RESOLUTION*1000,
78
- 0,
79
- 0,
80
- -RESOLUTION*1000,
81
- OFFSET_X*1000+RESOLUTION*500,
82
- OFFSET_Y*1000+RESOLUTION*500)
83
- }
@@ -1,99 +0,0 @@
1
- require 'carray'
2
- require 'carray/image_map'
3
- require 'sdi/hdf'
4
- require 'proj4'
5
- require 'benchmark'
6
-
7
- RESOLUTION = 10 # (km) on map projection
8
- UL_CORNER = [3000, -3000] # (km) on map projection
9
- LR_CORNER = [-3000, 3000] # (km) on map projection
10
- SCANLINES = 2
11
-
12
- PROJECTION = "+proj=stere +lat_0=-90 +units=km"
13
-
14
- IMAGE_HEIGHT = (UL_CORNER[0] - LR_CORNER[0])/RESOLUTION.to_f + 1
15
- IMAGE_WIDTH = (LR_CORNER[1] - UL_CORNER[1])/RESOLUTION.to_f + 1
16
-
17
- SCALE_Y = -(UL_CORNER[0] - LR_CORNER[0])/(IMAGE_HEIGHT.to_f - 1)
18
- SCALE_X = (LR_CORNER[1] - UL_CORNER[1])/(IMAGE_WIDTH.to_f - 1)
19
- OFFSET_Y = UL_CORNER[0].to_f # LOWERRIGHT[0].to_f
20
- OFFSET_X = UL_CORNER[1].to_f
21
-
22
- image = ImageMap.new(CA_INT32, [IMAGE_HEIGHT, IMAGE_WIDTH]) { -9999 }
23
- image.set_xrange(UL_CORNER[1]..LR_CORNER[1])
24
- image.set_yrange(LR_CORNER[0]..UL_CORNER[0])
25
-
26
- #
27
- #
28
- #
29
-
30
- sdi = SDI::HDF.open(ARGV[0])
31
- clat = sdi.open("Latitude")[]
32
- clon = sdi.open("Longitude")[]
33
- height = sdi.open("Height")[]
34
- sdi.close
35
-
36
- NSCANS = clat.dim0/SCANLINES
37
- NPIXELS = clat.dim1
38
-
39
- #cidx = clat.index
40
- #cidx = clat.template.random!(1000)
41
- cidx = height.copy
42
-
43
- proj = Proj4.new(PROJECTION)
44
-
45
- Benchmark.bm do |rep|
46
- rep.report do
47
- (NSCANS).times do |ns|
48
-
49
- # printf("%i/%i\n", ns, NSCANS)
50
-
51
- line = ns * SCANLINES
52
-
53
- bidx = cidx[[line, SCANLINES], nil].copy
54
- blat = clat[[line, SCANLINES], nil].copy
55
- blon = clon[[line, SCANLINES], nil].copy
56
-
57
- image.warp(bidx, blon, blat, :grid => "area", :gradation => false) {
58
- # image.map_points(bidx, blon, blat) {
59
- |lo,la,gx,gy|
60
- proj.fwd_f(lo, la, gx, gy)
61
- }
62
-
63
- end
64
- end
65
- end
66
-
67
- puts "projection done"
68
-
69
- #image = height.map_object(image, -9999)
70
- image.mask = image < -9998
71
-
72
- puts "masking flipping done"
73
-
74
- #sc = nil
75
- sc = ColorScale.GMT("GMT_topo.cpt")
76
- #sc = ImageMap::ColorScale.GMT("DEM_poster.cpt")
77
-
78
- #colors = CArray.float(5,5)
79
- #colors[0,nil] = [0, 0, 0, 0, 1]
80
- #colors[1,nil] = [1990, 0, 0, 0, 1]
81
- #colors[2,nil] = [2000, 1, 1, 1, 0]
82
- #colors[3,nil] = [2010, 0, 0, 0, 1]
83
- #colors[4,nil] = [5000, 0, 0, 0, 1]
84
- #sc = ImageMap::ColorScale.mapping(colors)
85
- #p sc.scale
86
- #p sc.palette
87
-
88
- ImageMap.display(image, sc)
89
-
90
- open("projection.tfw", "w") { |io|
91
- io.printf([["%g\n"]*6].join(),
92
- RESOLUTION*1000,
93
- 0,
94
- 0,
95
- -RESOLUTION*1000,
96
- OFFSET_X*1000+RESOLUTION*500,
97
- OFFSET_Y*1000+RESOLUTION*500)
98
- }
99
-