Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1001 1404 (N27W10 L=120)   453    156        0 ( 918/1377-1077/1443)  -9.5 W03-W15 N25-N28
      2  738 1411 (N35E07 L=103)   477    143        1 ( 669/1389- 808/1427) -16.7 E12-E02 N32-N36
      3 1469 1379 (N12W36 L=146)   914    306        0 (1324/1336-1603/1409) -34.5 W27-W44 N06-N17
      4 1815 1334 (S01W63 L=173)  1100    205        3 (1756/1300-1867/1364) -22.6 W56-W69 S04-N01
      5 1214 1276 (N14W18 L=128)   441     82        1 (1183/1250-1265/1310) -45.0 W15-W20 N11-N16
      6 1759 1097 (S12W48 L=158)  1630    384        4 (1733/ 914-1800/1280)  66.8 W43-W57 S23-S02
      7  226  924 (N18E56 L= 54)   759    190        2 ( 181/ 882- 286/ 978) -36.9 E59-E51 N14-N21
      8 1502  870 (S16W23 L=134)   631    184        2 (1414/ 821-1567/ 917) -55.0 W19-W26 S21-S11
      9  728  698 (S05E25 L= 85)   265     81        0 ( 691/ 674- 768/ 718) -36.9 E27-E23 S07-S03
     10 1164  438 (S33E09 L=101)   464     86        2 (1143/ 403-1186/ 474)  51.3 E11-E06 S35-S30
     11  633  363 (S23E46 L= 64)   892    229        1 ( 607/ 290- 663/ 437)  57.5 E49-E40 S29-S18
     12 1317  319 (S48E03 L=107)  3000    420        6 (1159/ 240-1415/ 410)   6.8 E18-W07 S53-S42
  

    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