Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  649 1600 (N36E35 L= 66)   354    105        0 ( 614/1557- 678/1638)  74.1 E36-E34 N32-N39
      2  299 1530 (N27E70 L= 31)  2129    136        5 ( 250/1483- 336/1594) -38.7 E80-E61 N23-N31
      3  964 1443 (N28E07 L= 94)  1008    199        6 ( 901/1362-1027/1509) -48.4 E12-E02 N22-N31
      4 1251 1389 (N26W13 L=114)   405    143        0 (1210/1319-1278/1436) -63.4 W09-W14 N21-N29
      5 1542 1310 (N23W34 L=135)   247     50        1 (1520/1284-1563/1337) -56.3 W32-W35 N21-N24
      6 1819 1139 (N14W59 L=160)   517     86        2 (1797/1096-1842/1178) -90.0 W56-W61 N11-N16
      7  392 1115 (N02E43 L= 58)   259     43        1 ( 374/1091- 418/1139) -45.0 E44-E41 N01-N03
      8  142  914 (S12E71 L= 30)  1288    141        3 ( 112/ 863- 178/ 977) -41.2 E77-E65 S15-S08
      9  180  937 (S11E64 L= 37)   544     80        2 ( 154/ 918- 214/ 973) -15.9 E68-E59 S11-S08
     10  369  923 (S10E44 L= 57)   345    132        0 ( 322/ 877- 422/ 958)  33.7 E48-E39 S12-S06
     11 1179  673 (S19W13 L=114)  2810    704        7 ( 911/ 520-1397/ 899)  46.0 E04-W24 S31-S03
     12  709  806 (S14E18 L= 83)   241     79        0 ( 663/ 789- 746/ 819)   0.0 E21-E15 S14-S13
     13 1761  730 (S12W57 L=158)   406     75        2 (1736/ 702-1791/ 751)  36.9 W53-W59 S13-S10
     14  294  627 (S30E57 L= 44)  1764    208        2 ( 231/ 542- 354/ 686)  22.1 E79-E47 S37-S24
     15 1716  574 (S22W58 L=160)  3033    368        9 (1635/ 508-1789/ 678) -31.4 W48-W74 S26-S15
     16  593  439 (S41E31 L= 71)  1091    126        2 ( 543/ 391- 653/ 484)  26.6 E36-E24 S45-S36
     17  758  283 (S53E19 L= 83)   429     42        1 ( 738/ 259- 782/ 307) -63.4 E20-E16 S55-S50
  

    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