Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  483 1521 (N18E44 L=327)   547    125        2 ( 425/1497- 547/1538)  11.3 E48-E38 N15-N18
      2  829 1475 (N18E18 L=352)   510    131        1 ( 773/1441- 873/1501)   0.0 E21-E15 N16-N19
      3  860 1385 (N13E15 L=356)  1225    285        7 ( 750/1313- 959/1451)  45.0 E20-E09 N07-N18
      4  529 1420 (N12E37 L=333)   247     67        1 ( 509/1386- 550/1444)  71.6 E38-E35 N09-N13
      5  626 1385 (N10E29 L=341)   380    112        1 ( 587/1335- 653/1417) -38.7 E32-E26 N07-N12
      6  794 1036 (S08E14 L=357)   288    112        1 ( 755/ 999- 823/1082) -36.9 E16-E11 S10-S05
      7  625 1037 (S10E24 L=346)   326    118        0 ( 593/ 984- 647/1081)  90.0 E25-E22 S13-S07
      8 1499  997 (S01W29 L= 39)   199     77        0 (1488/ 955-1512/1040) -78.7 W27-W29 S03-N02
      9  150  900 (S22E66 L=304)  1978    203        5 (  87/ 840- 210/ 970)  16.9 E82-E57 S25-S17
     10  559  788 (S26E27 L=343)  1229    330        3 ( 413/ 746- 711/ 841)  10.3 E38-E16 S28-S21
     11 1304  691 (S21W23 L= 33)  1065    414        1 (1186/ 569-1469/ 814)  57.7 W18-W30 S30-S11
     12  993  771 (S21W02 L= 12)   676    298        1 ( 859/ 747-1130/ 794)  21.8 E06-W09 S24-S17
     13 1537  635 (S21W40 L= 50)   680    144        1 (1509/ 571-1558/ 715) -66.0 W37-W42 S24-S15
     14 1734  591 (S19W59 L= 70)   379     68        2 (1704/ 568-1761/ 619)  53.1 W56-W61 S20-S16
     15  974  559 (S35W04 L= 14)   848    228        0 ( 858/ 539-1087/ 581)   7.1 E04-W12 S36-S33
     16 1657  528 (S24W55 L= 66)   308     42        1 (1631/ 509-1683/ 544)   0.0 W52-W58 S25-S23
     17 1399  294 (S45W50 L= 60)  1915    227        5 (1298/ 249-1481/ 344)  12.5 W37-W61 S47-S40
     18  956  283 (S56W11 L= 21)   305     73        1 ( 926/ 263- 995/ 307) -12.5 W06-W15 S56-S53
     19 1545  261 (S43W72 L= 82)   915     70        1 (1532/ 222-1562/ 302) -14.7 W64-W85 S44-S39
  

    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