Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1171 1488 (N26W15 L= 50)   997    231        2 (1052/1435-1264/1519) -19.7 W06-W20 N21-N28
      2  483 1348 (N24E33 L=  2)   904    181        1 ( 406/1298- 569/1392) -26.6 E40-E26 N20-N27
      3 1013 1157 (N07W01 L= 36)   279    100        1 ( 975/1116-1050/1207) -63.4 E01-W02 N04-N10
      4  779 1074 (N04E14 L= 21)   244     95        0 ( 761/1021- 795/1115)  68.2 E15-E13 N01-N06
      5 1862 1058 (S06W61 L= 96)   468     75        2 (1847/1017-1876/1097)  76.0 W58-W62 S08-S04
      6 1908  911 (S16W68 L=103)  1395    113        3 (1867/ 868-1943/ 955)  14.0 W61-W74 S18-S12
      7 1687  931 (S13W43 L= 78)   252     64        1 (1656/ 903-1713/ 947)   0.0 W40-W45 S14-S11
      8 1491  671 (S27W28 L= 63)   260     94        0 (1480/ 622-1499/ 719) -81.9 W27-W29 S30-S23
      9  620  515 (S28E34 L=  1)  1287    273        0 ( 552/ 408- 685/ 627)  37.4 E44-E26 S34-S21
  

    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