Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  399 1566 (N42E55 L=244)   458     78        0 ( 356/1550- 434/1590)  12.5 E61-E51 N41-N44
      2 1041 1558 (N39W05 L=305)   238     83        0 (1003/1541-1080/1576)   9.5 W02-W08 N37-N40
      3 1786 1442 (N23W68 L=  7)  2633    264        4 (1696/1373-1858/1532) -45.0 W59-W79 N18-N30
      4  345 1180 (N17E47 L=253)   431     70        2 ( 324/1137- 374/1216)  53.1 E48-E44 N14-N19
      5 1172 1130 (N10W10 L=309)   677    260        0 (1079/1028-1241/1180) -41.6 W04-W13 N03-N13
      6  176 1108 (N12E65 L=234)   571     60        1 ( 156/1080- 203/1143) -26.6 E69-E61 N10-N14
      7  822  923 (S00E13 L=287)  1514    450        0 ( 743/ 779- 946/1098) -59.9 E17-E06 S09-N10
      8  495 1049 (N09E34 L=265)   226     51        0 ( 483/1020- 510/1079)  76.0 E35-E33 N06-N10
      9 1609 1028 (N01W38 L=338)   321    118        0 (1586/ 972-1634/1075)  51.3 W36-W40 S02-N03
     10 1649  965 (S04W41 L=341)   805    228        1 (1567/ 920-1725/1045)  15.9 W34-W48 S05-N00
     11 1681  846 (S11W44 L=343)   343     97        0 (1649/ 809-1725/ 873)  26.6 W40-W47 S13-S09
     12 1776  775 (S17W54 L=353)   369     92        0 (1750/ 739-1790/ 819)  59.0 W51-W55 S19-S14
     13  823  740 (S12E14 L=285)   246     89        0 ( 810/ 690- 839/ 786)  80.5 E15-E12 S14-S08
     14 1070  702 (S15W01 L=301)   247    128        0 (1016/ 668-1120/ 731) -31.0 E02-W03 S17-S13
     15  380  663 (S15E48 L=252)   627    102        2 ( 343/ 619- 416/ 712) -80.5 E50-E45 S18-S12
     16 1417  569 (S27W24 L=324)   299     64        0 (1389/ 551-1451/ 587)  26.6 W22-W26 S28-S25
     17  716  552 (S23E24 L=276)   228     60        0 ( 684/ 540- 750/ 563)  11.3 E26-E21 S24-S22
     18 1198  334 (S44W09 L=309)  2292    214        7 (1126/ 263-1273/ 402)  13.0 W02-W15 S50-S37
     19  691  256 (S48E40 L=260)  3021    286        6 ( 567/ 216- 822/ 291)   8.3 E66-E24 S52-S44
  

    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