Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1346 1688 (N35W47 L=149)  3342    531        6 (1203/1527-1548/1763)  -7.1 W30-W71 N27-N39
      2  601 1134 (N23E23 L= 79)   324     39        0 ( 577/1107- 626/1164) -53.1 E24-E21 N21-N25
      3  196 1012 (N24E62 L= 40)   705    124        0 ( 143/ 968- 247/1040) -18.4 E70-E55 N20-N26
      4 1729  891 (S21W42 L=144)   891    200        0 (1646/ 862-1795/ 931) -28.3 W35-W48 S23-S15
      5  238  878 (N16E59 L= 44)   338     59        0 ( 222/ 847- 253/ 906)  63.4 E60-E56 N13-N17
      6 1859  851 (S28W56 L=158)  1276    159        0 (1799/ 815-1924/ 900)  -5.7 W48-W68 S30-S26
      7  170  800 (N12E70 L= 32)  2814    183        2 ( 113/ 742- 223/ 869) -29.7 E80-E63 N07-N15
      8 1427  754 (S19W17 L=119)   399    100        1 (1409/ 698-1450/ 799)  56.3 W14-W18 S22-S16
      9 1768  697 (S35W44 L=147)  1246    147        4 (1713/ 626-1808/ 763)  31.0 W38-W49 S39-S31
     10  379  572 (S05E56 L= 47)   543    110        0 ( 323/ 554- 420/ 592)  -5.7 E60-E50 S06-S03
     11 1276  475 (S33W00 L=102)   314     90        0 (1248/ 439-1297/ 519)  45.0 E02-W02 S34-S29
     12  760  352 (S27E38 L= 64)   346     46        0 ( 731/ 326- 789/ 376)  11.3 E40-E35 S28-S25
  

    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