Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1353 1831 (N42W65 L=171)   514     44        0 (1339/1802-1365/1858)  24.0 W60-W71 N39-N43
      2 1550 1535 (N19W51 L=156)   742    188        1 (1516/1456-1586/1612)  59.0 W48-W55 N13-N23
      3  651 1562 (N50E10 L= 95)   307     40        1 ( 629/1537- 670/1587) -76.0 E12-E08 N47-N51
      4  337 1557 (N59E55 L= 51)   496     59        1 ( 302/1543- 372/1567)  -5.7 E66-E46 N57-N60
      5 1504 1532 (N20W47 L=153)   385     60        2 (1485/1490-1523/1561)  68.2 W45-W49 N17-N22
      6  271 1295 (N41E52 L= 53)   640    128        0 ( 259/1228- 281/1355) -83.7 E55-E50 N36-N45
      7  509 1092 (N23E30 L= 75)   896    259        1 ( 415/1003- 603/1158) -42.5 E36-E24 N16-N28
      8  239 1050 (N26E55 L= 50)  1336    146        4 ( 183/1011- 296/1111)  24.4 E62-E49 N23-N29
      9 1626  969 (S14W34 L=140)   963    259        0 (1503/ 903-1720/1024)  -6.3 W24-W42 S14-S12
     10  276  981 (N21E52 L= 53)   512     94        0 ( 257/ 945- 296/1023)  53.1 E54-E49 N19-N23
     11 1859  979 (S21W57 L=162)   853     99        0 (1828/ 935-1885/1015)  14.0 W52-W60 S23-S18
     12  263  920 (N18E55 L= 51)  1049    186        2 ( 213/ 842- 322/ 982) -68.2 E60-E49 N12-N22
     13 1385  628 (S26W11 L=116)   899    282        0 (1315/ 524-1435/ 750)  40.2 W03-W16 S31-S19
     14 1140  543 (S24E07 L= 99)  1456    562        0 ( 972/ 358-1318/ 714)  22.1 E24-W08 S32-S18
     15 1600  682 (S30W27 L=133)   305     85        0 (1569/ 650-1631/ 711)  18.4 W24-W30 S31-S28
     16 1780  637 (S41W46 L=152)  1606    233        1 (1694/ 574-1860/ 710)   5.0 W35-W58 S43-S38
     17  499  342 (S21E61 L= 45)   756     49        0 ( 472/ 317- 525/ 373)  10.3 E66-E55 S23-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