Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  882 1574 (N21E25 L=281)   501    108        1 ( 842/1524- 906/1618) -33.7 E29-E21 N17-N22
      2 1679 1278 (N28W34 L=340)  1598    379        0 (1645/1137-1746/1430) -53.1 W29-W42 N19-N37
      3 1566 1315 (N26W25 L=330)   258     89        0 (1550/1268-1589/1351)  90.0 W23-W25 N23-N29
      4 1870 1246 (N34W57 L=  3)   874    148        1 (1825/1175-1901/1300) -20.0 W50-W62 N30-N36
      5  178  995 (S29E63 L=242)  1645    196        1 (  86/ 970- 280/1019)   4.8 E86-E50 S30-S25
      6  839  832 (S22E05 L=301)   286    117        0 ( 819/ 786- 867/ 869) -45.0 E06-E02 S24-S19
      7 1716  843 (N05W48 L=354)   595     89        0 (1669/ 822-1763/ 861)  29.7 W44-W52 N02-N06
      8 1130  647 (S24W18 L=324)   611    224        0 (1030/ 622-1240/ 677)  36.0 W12-W23 S28-S19
      9 1470  551 (S19W43 L=349)   411     66        2 (1446/ 521-1491/ 585) -76.0 W40-W45 S20-S16
     10 1556  499 (S18W52 L=358)   299     54        1 (1529/ 480-1581/ 517)  90.0 W49-W54 S19-S16
     11 1545  423 (S22W57 L=  3)  1520    176        2 (1445/ 394-1607/ 467)  38.7 W50-W63 S27-S18
  

    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