Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  812 1858 (N72E21 L=151)   605     43        1 ( 791/1826- 830/1885) -45.0 E28-E15 N67-N76
      2 1339 1499 (N32W28 L=200)   413    136        0 (1265/1484-1398/1512) -11.3 W22-W32 N30-N33
      3  837 1480 (N36E10 L=162)   629    156        1 ( 768/1432- 887/1505) -18.4 E15-E06 N32-N37
      4 1437 1436 (N27W34 L=206)   430    125        0 (1407/1387-1456/1482)  56.3 W30-W36 N24-N30
      5 1237 1438 (N29W19 L=191)   340     82        2 (1219/1406-1269/1478) -90.0 W17-W20 N27-N32
      6 1593 1419 (N24W46 L=218)   381     96        3 (1557/1398-1648/1446) -33.7 W43-W50 N21-N25
      7  230 1312 (N27E65 L=107)  1215    142        2 ( 168/1274- 291/1343)  -7.1 E79-E55 N24-N28
      8 1860 1258 (N10W69 L=241)  1662    166        1 (1832/1184-1894/1330) -90.0 W65-W73 N04-N13
      9  922 1185 (N16E05 L=167)   655    245        3 ( 858/1109- 994/1244) -32.0 E09-W00 N11-N19
     10 1751  858 (S13W50 L=222)   506      5        0 (1725/ 837-1777/ 877)  -9.5 W47-W53 S14-S11
     11 1577  780 (S15W35 L=207)   305    119        0 (1532/ 746-1625/ 820)  20.6 W31-W39 S17-S13
     12 1701  662 (S25W48 L=220)  2034    288        2 (1613/ 539-1757/ 769)  57.3 W41-W52 S32-S18
     13  322  646 (S15E56 L=116)  1497    242        3 ( 266/ 577- 369/ 717) -66.0 E60-E51 S19-S10
     14 1138  639 (S20W04 L=176)   247     71        1 (1111/ 613-1172/ 658)  26.6 W02-W06 S21-S18
     15  394  599 (S18E50 L=122)   390     62        0 ( 368/ 579- 415/ 625)  26.6 E53-E47 S19-S16
     16  674  410 (S33E33 L=139)   934    261        2 ( 566/ 371- 780/ 454)  10.8 E44-E23 S35-S30
     17  995  217 (S54E14 L=158)   479     68        1 ( 962/ 202-1036/ 239)  -7.1 E17-E09 S55-S51
  

    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