Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  932 1871 (N62W10 L=161)   431     68        0 ( 909/1843- 967/1892)  14.0 W07-W15 N58-N64
      2 1027 1542 (N32W07 L=158)   510    156        0 ( 967/1486-1087/1594) -53.1 W04-W10 N27-N36
      3 1868 1310 (N06W68 L=219)  1027     73        2 (1845/1279-1903/1354) -33.7 W64-W74 N04-N09
      4  389 1245 (N21E40 L=110)  1649    217        6 ( 294/1187- 491/1293)   3.4 E48-E31 N17-N23
      5 1855  986 (S12W59 L=210)  1106    136        2 (1797/ 951-1908/1043)  -7.6 W52-W67 S14-S09
      6 1590  869 (S16W35 L=185)  1614    390        1 (1543/ 698-1626/1031)  78.7 W31-W37 S26-S06
      7 1881  852 (S21W64 L=215)  3610    338       12 (1822/ 729-1933/ 949)  85.6 W56-W73 S28-S15
      8 1481  772 (S21W26 L=177)   272    106        0 (1448/ 731-1504/ 822)  51.3 W23-W28 S23-S18
      9 1709  595 (S36W49 L=199)   902    163        0 (1649/ 534-1766/ 637)  31.0 W43-W54 S39-S33
     10  593  579 (S22E35 L=115)   230     79        0 ( 575/ 542- 622/ 612)  33.7 E37-E31 S23-S19
     11  312  569 (S19E60 L= 90)   481     44        2 ( 292/ 546- 333/ 595) -76.0 E62-E57 S20-S16
  

    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