Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1198 1751 (N53W01 L=161)   582     92        1 (1159/1730-1251/1767)  35.5 E03-W05 N49-N55
      2  589 1728 (N39E52 L=108)   515     84        1 ( 556/1711- 640/1745)  23.2 E55-E47 N36-N40
      3  900 1523 (N29E16 L=144)  2011    436        3 ( 725/1426-1083/1601)  -6.1 E30-E02 N25-N32
      4 1706 1549 (N45W55 L=215)   939     91        1 (1662/1530-1754/1583)  10.0 W47-W64 N43-N47
      5 1211 1322 (N21W07 L=167)   194     60        0 (1179/1299-1239/1338) -26.6 W05-W09 N19-N21
      6 1103 1260 (N15W02 L=162)   468    221        0 (1029/1165-1158/1297) -29.1 E03-W06 N10-N17
      7  186 1227 (N00E66 L= 94)  1683    144        3 ( 146/1156- 231/1289)  74.1 E72-E61 S03-N04
      8 1725 1145 (N17W47 L=207)   378     61        0 (1697/1127-1758/1161)  21.8 W44-W49 N15-N18
      9 1087  785 (S13W07 L=167)  1003    370        1 ( 940/ 729-1235/ 855)  34.5 E01-W15 S18-S07
     10  483  788 (S22E33 L=127)   344     74        1 ( 452/ 761- 518/ 823)  45.0 E35-E30 S23-S19
     11 1199  627 (S21W17 L=177)  2145    602        1 ( 949/ 569-1398/ 779)  35.8 W01-W26 S28-S09
     12 1068  268 (S50W21 L=181)   383     59        0 (1042/ 239-1089/ 291) -33.7 W17-W23 S53-S48
  

    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