Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  411 1604 (N15E62 L= 18)   515     73        1 ( 390/1570- 438/1638) -15.3 E67-E56 N12-N16
      2 1309 1531 (N34W03 L= 83)   261     58        0 (1286/1505-1334/1552) -26.6 W01-W05 N32-N35
      3  282 1523 (N07E72 L=  8)   655     66        2 ( 256/1494- 304/1551) -53.1 E74-E69 N05-N09
      4 1892 1032 (N22W64 L=144)  3333    283        3 (1823/ 911-1952/1131)  42.0 W57-W84 N12-N30
      5  629 1023 (S14E23 L= 57)   240     90        1 ( 588/1009- 668/1048)   0.0 E26-E20 S15-S12
      6  539  702 (S36E24 L= 56)  2211    405        6 ( 441/ 590- 668/ 842)  57.5 E32-E16 S45-S23
      7 1574  645 (S10W45 L=125)   792    265        0 (1472/ 591-1654/ 700)  50.7 W40-W49 S15-S04
      8 1261  490 (S28W32 L=112)   899    215        2 (1161/ 434-1353/ 543)  59.7 W28-W36 S34-S22
      9 1548  492 (S19W52 L=132)   787    142        2 (1505/ 439-1613/ 537)  63.4 W48-W56 S23-S15
     10 1109  423 (S37W26 L=106)   511     83        2 (1078/ 394-1141/ 462) -23.2 W22-W29 S39-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