Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1130 1877 (N57E31 L=212)  1541    145        1 (1068/1834-1188/1918) -15.3 E43-E21 N53-N61
      2  741 1707 (N29E42 L=202)  1758    392        2 ( 651/1584- 840/1835) -20.4 E63-E28 N22-N36
      3 1334 1792 (N57E05 L=238)   386     49        0 (1313/1770-1362/1809)  59.0 E07-E02 N54-N59
      4 1159 1551 (N30E07 L=236)  1056    208        2 (1106/1449-1206/1625) -42.3 E12-E01 N24-N34
      5  820 1494 (N17E25 L=218)   445    170        1 ( 767/1434- 874/1531) -16.7 E29-E19 N14-N19
      6 1106 1318 (N14E03 L=240)   451    136        0 (1080/1261-1150/1372)  90.0 E04-E01 N10-N18
      7 1693 1170 (N23W39 L=282)   332     92        0 (1664/1129-1719/1211) -26.6 W35-W41 N21-N24
      8  419 1049 (S18E39 L=205)   533    114        0 ( 362/1020- 471/1068)   0.0 E43-E34 S18-S16
      9  499  941 (S22E30 L=213)   543    111        1 ( 481/ 886- 521/ 990) -80.5 E31-E28 S25-S19
     10  582  867 (S25E23 L=220)   236     75        1 ( 546/ 850- 623/ 887)  31.0 E25-E20 S26-S22
     11  823  844 (S20E07 L=237)   661    215        1 ( 760/ 766- 897/ 878)  56.3 E09-E02 S25-S16
     12 1306  730 (S13W25 L=268)   411    136        0 (1270/ 671-1332/ 784) -45.0 W21-W28 S15-S10
  

    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