Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  316 1527 (N09E67 L=113)  1171    172        2 ( 272/1450- 361/1597) -78.7 E71-E64 N03-N13
      2  701 1509 (N16E33 L=147)   304     96        0 ( 678/1466- 720/1544) -38.7 E36-E30 N13-N17
      3 1481 1376 (N29W19 L=200)  1846    482        2 (1329/1245-1646/1514) -13.6 W05-W34 N25-N34
      4  296 1375 (S01E60 L=121)   518     50        1 ( 272/1351- 318/1408)  45.0 E61-E57 S02-N01
      5  852 1126 (S03E12 L=168)   530    159        3 ( 829/1046- 881/1200) -49.4 E15-E08 S06-N01
      6  197  965 (S27E61 L=120)   513     82        0 ( 185/ 921- 212/1008)  90.0 E62-E58 S29-S24
      7  773  674 (S31E06 L=175)   599    154        1 ( 716/ 646- 846/ 699)  18.4 E09-W00 S33-S28
      8  852  580 (S34W02 L=183)   305    126        0 ( 814/ 531- 877/ 633)  90.0 W00-W02 S38-S30
      9 1635  489 (S17W60 L=241)   412     47        1 (1608/ 475-1665/ 504)  18.4 W57-W63 S17-S15
  

    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