Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1117 1525 (N27E08 L= 33)   333     53        2 (1092/1503-1145/1548) -14.0 E10-E06 N26-N29
      2 1161 1410 (N21E02 L= 39)   206     63        1 (1118/1388-1187/1434)   0.0 E04-W00 N20-N22
      3  700 1069 (S10E20 L= 22)   204     86        1 ( 670/1033- 721/1116)  90.0 E20-E18 S12-S06
      4 1854  972 (N16W58 L=100)  2040    242        1 (1812/ 865-1880/1063) -42.3 W51-W62 N09-N21
      5 1651  981 (N10W39 L= 81)   304    111        1 (1617/ 928-1673/1036) -35.5 W35-W42 N07-N12
      6  333  938 (S26E43 L=358)   662    156        1 ( 282/ 916- 388/ 983)  36.9 E48-E39 S27-S21
      7  735  800 (S24E11 L= 30)   334     94        1 ( 676/ 779- 781/ 821)   8.1 E15-E08 S24-S23
      8  418  701 (S39E34 L=  8)  1338    216        4 ( 331/ 614- 492/ 787) -23.6 E42-E26 S43-S34
      9 1319  712 (S14W26 L= 67)   410    147        1 (1290/ 636-1369/ 772) -83.7 W24-W27 S18-S09
     10  632  628 (S38E15 L= 27)   536    115        1 ( 572/ 611- 701/ 651)  10.3 E20-E09 S39-S35
     11 1109  626 (S24W16 L= 58)   343    132        0 (1043/ 611-1186/ 637)  24.0 W12-W21 S26-S22
     12 1541  396 (S25W56 L= 98)  1205    269        1 (1419/ 358-1643/ 451)  40.2 W48-W62 S30-S18
  

    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