Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1593 1634 (N26W58 L=200)   599     42        1 (1566/1609-1624/1660)   7.1 W54-W62 N25-N27
      2  909 1337 (N22E01 L=141)   193     14        0 ( 888/1317- 928/1353) -63.4 E02-W00 N20-N23
      3  246 1233 (N28E54 L= 88)   524     73        2 ( 225/1205- 266/1275) -45.0 E57-E51 N26-N31
      4 1794  900 (S21W51 L=193)   937    161        0 (1723/ 844-1852/ 962)  17.1 W43-W58 S22-S17
      5 1604  924 (S15W34 L=176)   355     72        0 (1578/ 889-1628/ 955)  36.9 W32-W36 S16-S13
      6  480  620 (S12E43 L= 98)  1004    175        1 ( 455/ 529- 510/ 697)  48.0 E48-E38 S17-S07
      7 1707  649 (S35W43 L=185)   364     66        1 (1677/ 628-1739/ 670) -20.6 W39-W47 S37-S32
      8 1448  565 (S35W20 L=161)   852    230        1 (1358/ 528-1561/ 611) -32.0 W13-W29 S39-S29
      9 1593  552 (S39W33 L=174)   309     68        1 (1570/ 525-1618/ 580)  36.9 W30-W34 S40-S37
     10 1217  470 (S36W01 L=143)   828    177        0 (1139/ 408-1284/ 537)  24.4 E05-W06 S39-S32
     11 1058  490 (S31E09 L=132)   586    103        1 (1022/ 448-1097/ 533)  26.6 E11-E05 S34-S28
  

    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