Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1033 1598 (N35W10 L=151)   582    143        1 ( 967/1572-1108/1617) -18.4 W05-W14 N32-N37
      2  846 1544 (N35E04 L=137)   604    220        2 ( 749/1485- 923/1607)  24.8 E11-W02 N31-N38
      3 1814  939 (S16W54 L=195)   704    186        3 (1778/ 876-1861/1034)  52.1 W50-W60 S19-S10
      4 1568  860 (S17W32 L=174)  1386    442        3 (1448/ 719-1654/1015)  41.4 W23-W40 S24-S08
      5  557  687 (S14E35 L=106)   524     92        2 ( 514/ 644- 590/ 731)  45.0 E38-E31 S16-S11
      6 1017  624 (S24E06 L=135)   614    242        1 ( 969/ 516-1080/ 715)  47.7 E11-W00 S30-S19
      7  270  598 (S15E64 L= 77)  2109    128        3 ( 225/ 537- 318/ 657)  76.0 E70-E59 S19-S11
      8  267  520 (S20E71 L= 71)   843     54        0 ( 239/ 497- 290/ 546) -45.0 E74-E66 S21-S17
      9  739  441 (S32E30 L=111)   452    121        1 ( 708/ 400- 786/ 492)  26.6 E34-E24 S34-S29
     10  618  289 (S41E50 L= 91)   507     49        2 ( 590/ 257- 640/ 316) -90.0 E52-E47 S43-S38
  

    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