Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  588 1717 (N31E53 L=344)  1202    126        2 ( 542/1656- 621/1761) -18.4 E61-E46 N27-N34
      2  392 1607 (N20E63 L=334)  1275    172        3 ( 326/1566- 475/1657)  36.0 E71-E56 N16-N24
      3  169 1219 (S07E68 L=328)  2385    154        5 ( 119/1156- 224/1286)  38.7 E77-E62 S11-S02
      4 1372 1185 (N14W17 L= 54)   369    133        2 (1343/1136-1393/1246) -81.9 W15-W18 N11-N18
      5 1046 1075 (N01W00 L= 37)   966    244        2 ( 972/ 978-1118/1168)  72.9 E02-W02 S05-N08
      6  179  976 (S21E63 L=334)   494     68        0 ( 158/ 940- 203/1006)  38.7 E65-E59 S23-S19
      7 1302  863 (S06W20 L= 56)   606    137        0 (1247/ 825-1363/ 919)  60.3 W16-W22 S08-S01
      8 1587  880 (N02W38 L= 74)   280     48        1 (1570/ 854-1609/ 913)  90.0 W36-W39 N00-N03
      9 1780  756 (S01W58 L= 95)  1265    142        2 (1728/ 707-1818/ 814)  74.1 W54-W61 S05-N02
     10  223  781 (S33E58 L=339)  1259    104        1 ( 170/ 754- 280/ 809)   3.6 E67-E50 S34-S30
     11  610  575 (S38E19 L= 18)   244     52        0 ( 578/ 554- 640/ 597)  38.7 E21-E16 S40-S36
  

    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