Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1279 1502 (N24W29 L=213)  3245    775        9 ( 959/1446-1567/1530) -24.8 W09-W48 N15-N33
      2  774 1335 (N29E07 L=177)   408    114        1 ( 742/1294- 814/1376) -11.5 E08-E05 N25-N31
      3 1881 1177 (S11W64 L=247)  1849    255        5 (1847/1114-1918/1247)  40.6 W57-W69 S14-S06
      4  278 1214 (N34E47 L=137)   396    101        1 ( 240/1188- 307/1230) -33.7 E50-E44 N31-N35
      5 1814 1135 (S11W54 L=238)   679    111        2 (1777/1098-1843/1171)  63.4 W50-W58 S13-S09
      6 1924 1067 (S19W65 L=249)  2512    332        0 (1863/ 997-1956/1147)  66.0 W58-W70 S23-S14
      7 1590  740 (S27W28 L=211)  2548    604        1 (1413/ 531-1694/1028)  37.2 W11-W40 S36-S13
      8 1860  992 (S21W55 L=239)   532     59        0 (1846/ 956-1873/1025) -90.0 W53-W57 S23-S18
      9  473  814 (N06E38 L=145)   328     63        1 ( 448/ 789- 504/ 842) -14.0 E40-E36 N04-N07
     10  653  707 (S04E30 L=154)   587    229        2 ( 572/ 623- 718/ 804) -77.0 E31-E27 S10-N03
     11 1421  697 (S24W15 L=199)   584    161        0 (1357/ 610-1486/ 784)  32.5 W09-W20 S28-S20
     12 1094  499 (S27E11 L=173)   906    352        2 (1034/ 352-1144/ 640)  46.8 E20-E03 S35-S19
     13  820  550 (S17E25 L=159)   288     99        0 ( 800/ 508- 839/ 602)  63.4 E26-E22 S19-S13
     14  704  353 (S27E41 L=143)   965    276        2 ( 661/ 284- 733/ 442)  32.7 E48-E34 S31-S20
     15  887  265 (S38E35 L=149)   556     86        0 ( 851/ 239- 927/ 288)   7.1 E38-E30 S39-S36
     16  753  238 (S36E46 L=138)   566    103        1 ( 721/ 211- 780/ 268)  20.6 E50-E42 S37-S34
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y