Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1060 1705 (N38E10 L=248)   271     65        1 (1030/1684-1089/1730) -21.8 E12-E07 N36-N39
      2 1591 1470 (N32W32 L=291)   354     84        0 (1561/1429-1618/1505) -45.0 W29-W34 N29-N33
      3 1772 1249 (N21W48 L=307)   670    150        0 (1732/1174-1795/1306) -49.4 W44-W50 N16-N23
      4  212 1161 (S08E60 L=199)   520     56        2 ( 188/1127- 235/1191) -53.1 E62-E56 S09-S05
      5  543 1076 (S10E30 L=228)   565    147        0 ( 504/1020- 604/1134)  58.0 E32-E26 S13-S05
      6  225  921 (S22E57 L=201)   812    149        0 ( 197/ 856- 248/ 990)  76.0 E61-E54 S25-S17
      7  443  884 (S22E36 L=222)   287     89        0 ( 406/ 851- 482/ 912)  33.7 E39-E33 S24-S20
      8 1170  868 (S13W11 L=270)   221     96        0 (1130/ 842-1217/ 888)  45.0 W09-W13 S15-S11
      9  799  717 (S28E10 L=249)  1498    470        2 ( 620/ 620-1028/ 773)  30.0 E22-W04 S36-S20
     10 1621  632 (S18W48 L=306)  1179    332        0 (1458/ 594-1758/ 681)  24.4 W36-W59 S23-S12
     11 1683  431 (S28W68 L=327)  1922    252        3 (1601/ 374-1787/ 522)  57.3 W61-W74 S33-S19
     12  965  227 (S61W18 L=277)  1206    121        2 ( 924/ 187-1021/ 286) -16.4 W10-W27 S65-S56
  

    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