Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1249 1813 (N51E13 L=276)   534     74        3 (1222/1789-1277/1840)  45.0 E15-E10 N48-N54
      2  465 1711 (N21E61 L=228)   558     89        0 ( 441/1690- 497/1726)  21.8 E64-E57 N19-N22
      3  659 1535 (N13E36 L=253)   855    337        1 ( 612/1417- 725/1656) -39.1 E43-E27 N07-N20
      4  971 1411 (N13E13 L=275)   280     91        1 ( 947/1355- 989/1451) -71.6 E14-E11 N09-N15
      5  251 1259 (S11E58 L=230)   341     25        1 ( 234/1236- 268/1283) -71.6 E60-E56 S12-S09
      6  289 1205 (S13E53 L=236)   480     83        1 ( 271/1160- 308/1246) -68.2 E55-E50 S15-S10
      7  175 1199 (S15E65 L=224)  1113    202        4 ( 141/1161- 214/1236) -24.0 E70-E59 S17-S12
      8  209 1113 (S20E58 L=231)   481     71        1 ( 185/1086- 229/1143)  76.0 E60-E56 S21-S17
      9  736 1088 (S11E18 L=271)   356     80        2 ( 700/1061- 768/1118)  63.4 E20-E16 S12-S08
     10  485 1074 (S17E34 L=255)   611    184        2 ( 420/1050- 585/1099)   9.5 E38-E26 S19-S14
     11 1160  976 (S06W09 L=297)   316    103        1 (1129/ 931-1203/1024) -23.2 W05-W12 S08-S04
     12  117  972 (S29E66 L=223)   752    223        0 (  83/ 941- 152/1004)  26.6 E69-E61 S30-S27
     13  236  985 (S27E53 L=236)   458    125        1 ( 197/ 969- 274/1002)   8.1 E56-E49 S28-S25
     14  775  898 (S20E11 L=277)  1320    424        1 ( 595/ 806- 970/ 979)   9.5 E22-W02 S24-S16
     15  409  879 (S30E36 L=252)   864    181        0 ( 340/ 848- 482/ 929)  -4.8 E42-E30 S31-S28
     16 1182  694 (S21W18 L=307)   339     92        0 (1141/ 663-1226/ 725)   0.0 W15-W21 S21-S19
     17  868  586 (S36W03 L=292)   739    193        2 ( 770/ 545- 967/ 615)   7.6 E04-W11 S38-S34
     18 1747  549 (S10W65 L=354)   858    163        4 (1707/ 462-1780/ 591)  90.0 W62-W68 S15-S07
     19 1684  526 (S14W60 L=348)   993    169        4 (1648/ 474-1716/ 577) -80.5 W57-W64 S16-S10
  

    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