Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1614 1663 (N32W61 L= 30)   553    192        0 (1587/1627-1633/1699)  14.0 W56-W64 N31-N34
      2 1396 1534 (N30W37 L=  6)  1140    206        2 (1322/1486-1485/1580)   0.0 W30-W45 N27-N32
      3  471 1481 (N44E37 L=292)   375     77        1 ( 436/1462- 507/1496)   9.5 E39-E33 N42-N45
      4  799 1198 (N20E11 L=318)   272    125        0 ( 737/1170- 858/1227) -26.6 E15-E07 N17-N22
      5  688 1077 (N15E20 L=309)   215     88        2 ( 651/1061- 739/1096)   0.0 E22-E16 N13-N16
      6 1762 1018 (S08W50 L= 19)   647    162        4 (1715/ 961-1807/1064) -48.8 W46-W53 S12-S04
      7 1670  936 (S11W41 L= 10)   608    173        1 (1645/ 855-1699/1024)  63.4 W37-W44 S15-S05
      8  213  920 (N10E61 L=268)  1220    229        2 ( 162/ 876- 267/ 986)  23.2 E68-E53 N07-N13
      9 1529  885 (S11W29 L=358)   992    193        3 (1478/ 801-1572/ 970)  52.1 W24-W33 S15-S06
     10 1837  894 (S18W56 L= 25)   557     86        3 (1825/ 845-1851/ 935)  59.0 W54-W57 S20-S15
     11 1825  725 (S29W55 L= 24)   636    104        2 (1803/ 697-1845/ 770) -71.6 W52-W57 S31-S25
     12 1275  419 (S37W07 L=336)  1525    379        4 (1151/ 314-1362/ 556)  37.9 E05-W14 S43-S28
  

    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