Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1325 1241 (N17W14 L=282)  3528    766        4 (1032/ 944-1555/1528) -30.2 E10-W34 N03-N30
      2  313 1341 (N05E56 L=212)   386     70        1 ( 290/1309- 340/1364) -50.2 E57-E53 N02-N06
      3  137 1007 (S18E68 L=200)  1399    217        3 ( 111/ 948- 166/1076)  82.9 E71-E64 S21-S13
      4  208 1048 (S14E59 L=209)   545     50        1 ( 184/1023- 229/1072)  36.9 E61-E56 S15-S12
      5  249  981 (S18E54 L=214)   427    114        0 ( 225/ 932- 267/1030)  71.6 E55-E51 S20-S14
      6 1324  764 (S11W23 L=291)   707    181        1 (1229/ 722-1414/ 797)  35.0 W18-W28 S15-S08
      7  570  768 (S25E25 L=243)   466    131        1 ( 505/ 737- 629/ 789)   0.0 E29-E20 S26-S23
      8 1063  697 (S20W09 L=277)   761    177        1 ( 984/ 676-1161/ 733)   5.2 W03-W14 S21-S18
      9 1470  309 (S37W54 L=322)   892    129        2 (1445/ 279-1499/ 361) -29.1 W48-W57 S38-S33
     10 1287  306 (S41W40 L=308)  2347    405        4 (1106/ 277-1440/ 343)   9.8 W24-W53 S44-S36
  

    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