Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  955 1814 (N59E17 L=167)   366     48        3 ( 934/1789- 976/1840) -71.6 E20-E14 N56-N62
      2  460 1606 (N36E55 L=129)   614    108        1 ( 433/1563- 495/1650) -19.7 E62-E48 N33-N39
      3 1370 1467 (N33W23 L=207)  1518    423        4 (1246/1324-1530/1557) -31.8 W13-W34 N23-N38
      4  684 1404 (N24E26 L=158)   544    176        0 ( 644/1330- 731/1478) -45.0 E30-E21 N19-N28
      5 1677 1259 (N20W46 L=230)   539    131        6 (1624/1239-1718/1307)  -7.1 W41-W50 N18-N23
      6 1733  823 (S07W52 L=236)  1288    170        3 (1668/ 762-1805/ 875)  35.0 W46-W58 S10-S03
      7 1485  625 (S20W35 L=219)   570    149        1 (1452/ 568-1527/ 688)  76.0 W32-W36 S24-S16
      8 1630  635 (S19W47 L=231)   307     62        0 (1593/ 612-1653/ 656)  21.8 W43-W48 S20-S17
      9  493  546 (S33E38 L=146)   676    133        1 ( 435/ 528- 557/ 566)   0.0 E44-E32 S33-S30
  

    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