Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1060 1852 (N57E28 L= 35)  2098    192        4 ( 962/1823-1154/1879)  35.5 E35-E19 N50-N62
      2 1531 1203 (N19W28 L= 91)   768    275        2 (1388/1163-1621/1237)   3.2 W17-W35 N16-N22
      3 1656 1176 (N21W39 L=101)   285     71        1 (1631/1147-1678/1204) -33.7 W36-W40 N18-N21
      4 1676 1093 (N16W41 L=104)   915    193        2 (1618/1014-1731/1136)  33.7 W36-W45 N11-N19
      5 1829 1043 (N17W57 L=120)   418     73        1 (1816/ 998-1843/1076) -63.4 W55-W59 N14-N18
      6  166  764 (S35E72 L=351)   786     55        0 ( 138/ 744- 193/ 785)   8.5 E84-E64 S36-S33
      7 1629  552 (S16W54 L=117)   595    110        1 (1597/ 517-1659/ 610) -81.9 W51-W57 S19-S12
      8 1721  570 (S13W63 L=125)   815     85        1 (1699/ 519-1744/ 601) -59.0 W58-W66 S16-S11
      9 1385  490 (S26W38 L=100)   982    210        1 (1308/ 461-1483/ 534)  41.6 W32-W42 S29-S20
     10 1082  264 (S50W30 L= 93)   396     24        0 (1064/ 242-1103/ 288) -15.9 W27-W34 S51-S47
     11 1269  214 (S48W53 L=116)   744     57        2 (1235/ 185-1296/ 249) -13.0 W47-W60 S50-S44
  

    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