Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1096 1731 (N39E20 L= 64)   532    115        0 (1051/1695-1140/1770) -11.3 E25-E15 N36-N41
      2  506 1744 (N23E65 L= 19)   700     52        1 ( 474/1721- 534/1764)  -5.2 E72-E60 N21-N24
      3 1304 1493 (N29W03 L= 87)  1009    249        1 (1219/1394-1348/1576) -32.5 E04-W07 N23-N33
      4  967 1419 (N15E14 L= 70)   228    123        0 ( 950/1363- 992/1473) -74.1 E15-E12 N11-N18
      5 1081 1327 (N12E05 L= 79)   390    133        3 (1035/1271-1111/1363) -38.7 E07-E02 N09-N14
      6 1624 1237 (N24W31 L=115)   272     41        2 (1597/1212-1644/1257)  36.9 W28-W32 N22-N25
      7 1330 1174 (N11W13 L= 97)   679    285        1 (1223/1063-1372/1244) -21.0 W04-W18 N05-N13
      8 1608 1079 (N14W33 L=117)   238     45        0 (1593/1049-1623/1108)  76.0 W32-W33 N12-N16
      9 1768  984 (N14W48 L=132)  1663    385        1 (1687/ 905-1861/1051)   0.0 W40-W60 N08-N17
     10 1695  469 (S14W68 L=152)  1590    240        5 (1639/ 387-1733/ 555) -57.5 W61-W76 S20-S09
     11 1220  510 (S28W29 L=113)   352     86        0 (1169/ 498-1263/ 523)  38.7 W26-W31 S30-S26
     12  636  476 (S50E10 L= 74)   435     87        1 ( 612/ 452- 669/ 515) -18.4 E12-E06 S52-S47
     13  769  488 (S45W01 L= 85)   247     28        0 ( 750/ 468- 795/ 506)  63.4 W00-W02 S46-S42
     14 1204  359 (S38W38 L=122)   469    131        0 (1151/ 312-1259/ 401)  90.0 W36-W39 S42-S33
     15 1329  331 (S35W49 L=133)  1083    263        0 (1216/ 270-1451/ 383)  57.3 W44-W53 S41-S27
     16 1457  293 (S32W63 L=147)  1308    169        4 (1419/ 236-1495/ 360) -16.7 W56-W76 S34-S27
  

    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