Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1439 1764 (N46W53 L=128)  2198    334        5 (1389/1690-1502/1838)  67.4 W49-W59 N39-N52
      2 1805 1330 (N11W63 L=138)  1544    271        4 (1740/1216-1878/1440) -60.6 W59-W69 N02-N18
      3 1726 1344 (N14W55 L=130)   365     86        2 (1703/1320-1757/1369)   9.5 W52-W58 N12-N15
      4  415 1337 (N32E41 L= 34)   291    104        2 ( 377/1324- 457/1358) -14.0 E44-E36 N30-N33
      5  240 1284 (N30E57 L= 18)  1990    455        6 ( 141/1253- 333/1321)  -2.7 E68-E47 N27-N31
      6 1455 1252 (N13W31 L=105)   239    102        0 (1408/1217-1489/1283)  18.4 W27-W33 N11-N14
      7  540 1057 (N13E31 L= 44)   218     45        1 ( 516/1033- 559/1073) -45.0 E32-E29 N11-N14
      8 1710  901 (S12W44 L=119)   386     54        1 (1691/ 877-1734/ 928) -56.3 W42-W46 S13-S10
      9 1775  669 (S29W51 L=126)  1392    194        3 (1715/ 594-1832/ 742)  38.7 W45-W57 S33-S24
     10 1588  553 (S33W35 L=109)   371     54        1 (1562/ 535-1614/ 576)  11.3 W32-W37 S34-S31
     11  913  497 (S26E15 L= 60)   865    241        1 ( 815/ 466-1005/ 547) -18.4 E20-E08 S27-S21
  

    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