Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1209 1594 (N39W01 L=195)   251     42        0 (1189/1577-1234/1618) -45.0 W00-W03 N37-N40
      2 1562 1033 (N09W34 L=228)   555    108        1 (1512/1012-1626/1062)  20.6 W30-W38 N07-N10
      3  123 1020 (S17E75 L=119)   901     68        4 ( 107/ 987- 143/1057)  45.0 E79-E70 S18-S14
      4  545  845 (S20E28 L=166)  1933    521        7 ( 374/ 717- 732/1023)  45.0 E42-E17 S31-S06
      5 1710  650 (S10W57 L=251)   402     49        1 (1692/ 625-1733/ 673) -71.6 W54-W59 S11-S08
      6 1439  623 (S17W35 L=229)   294     30        0 (1416/ 606-1466/ 644)   0.0 W33-W37 S18-S15
      7 1789  581 (S13W71 L=266)   908     77        2 (1755/ 564-1836/ 597)  10.0 W66-W83 S14-S11
      8 1519  463 (S26W49 L=243)   376    107        0 (1478/ 430-1551/ 511)  71.6 W47-W50 S28-S22
      9 1045  233 (S54W27 L=221)  1694    318        3 ( 906/ 210-1189/ 258)  24.6 W13-W40 S60-S48
  

    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