Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  363 1579 (N09E64 L=269)   622     80        2 ( 349/1539- 379/1627) -40.6 E68-E60 N06-N12
      2  444 1587 (N11E56 L=277)   359     51        1 ( 425/1559- 463/1613) -63.4 E58-E53 N09-N13
      3  572 1588 (N14E45 L=287)   297     58        0 ( 553/1566- 595/1609)  71.6 E47-E44 N12-N15
      4 1743 1479 (N45W41 L= 14)  1860    233        0 (1667/1386-1826/1557) -18.4 W31-W53 N41-N48
      5 1655 1357 (N32W31 L=  4)   871    157        0 (1607/1292-1691/1414) -24.0 W26-W35 N28-N35
      6 1587 1236 (N22W28 L=  1)   415     99        2 (1558/1187-1606/1287) -38.7 W24-W29 N19-N24
      7  328 1221 (S13E50 L=283)   282     36        0 ( 308/1203- 352/1243)   0.0 E52-E47 S13-S11
      8 1900 1061 (N24W62 L= 35)   653     68        0 (1886/1027-1916/1095) -36.9 W59-W65 N21-N25
      9  951  934 (S14E02 L=331)   214     55        0 ( 921/ 912- 980/ 956)  63.4 E03-W00 S15-S11
     10  786  825 (S24E08 L=325)   310    157        0 ( 712/ 801- 857/ 845)  36.9 E12-E04 S27-S21
     11 1229  685 (S19W22 L=355)   287    121        1 (1189/ 664-1278/ 722)  63.4 W19-W24 S21-S15
     12 1261  574 (S24W29 L=  1)   822    257        0 (1156/ 525-1378/ 607)  42.3 W23-W35 S30-S19
  

    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