Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1361 1784 (N56E04 L= 51)   382     41        0 (1343/1761-1379/1805) -53.1 E06-E01 N53-N58
      2  762 1650 (N23E36 L= 19)   739    170        0 ( 688/1618- 837/1680)  10.3 E41-E30 N21-N25
      3 1431 1572 (N40W09 L= 65)   728    129        2 (1383/1527-1488/1617)  -6.3 W05-W14 N37-N42
      4 1637 1314 (N29W31 L= 86)   769    259        1 (1567/1202-1693/1408) -36.9 W23-W37 N23-N33
      5 1369 1199 (N14W14 L= 70)   287    110        2 (1351/1143-1389/1237) -80.5 W13-W15 N10-N16
      6  482 1077 (S17E34 L= 21)   652    208        2 ( 392/1028- 555/1148)  45.0 E40-E30 S21-S10
      7 1798 1006 (N17W51 L=106)   548    112        2 (1774/ 949-1817/1051) -59.0 W48-W54 N14-N19
      8 1543  602 (S13W45 L=100)   786    260        2 (1494/ 496-1585/ 696) -65.2 W41-W48 S20-S07
      9 1420  340 (S31W54 L=109)   511     65        0 (1398/ 318-1450/ 371) -14.0 W50-W58 S32-S29
  

    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