Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1719 1469 (N20W60 L= 85)   548     81        2 (1695/1431-1743/1500)  78.7 W58-W62 N17-N22
      2  669 1334 (N30E20 L=  5)  1311    325        2 ( 521/1291- 818/1388)   0.0 E31-E08 N28-N32
      3 1285 1216 (N14W19 L= 44)   515    150        3 (1228/1181-1356/1242) -32.0 W15-W23 N10-N15
      4  172 1044 (N16E65 L=320)  1028    184        2 ( 137/1006- 212/1093) -35.5 E69-E59 N13-N18
      5 1652  974 (S07W40 L= 65)   834    245        4 (1609/ 857-1697/1082)  50.2 W35-W45 S13-S01
      6  658  714 (S08E28 L=357)  1403    467        7 ( 528/ 566- 709/ 876)  50.2 E40-E23 S16-N02
      7  918  794 (S06E10 L= 15)   545    198        2 ( 835/ 763-1002/ 822) -18.4 E15-E05 S09-S04
      8 1739  741 (S24W46 L= 71)   961    139        3 (1688/ 706-1784/ 776)  -5.7 W41-W51 S26-S21
      9 1541  711 (S22W30 L= 54)   370    119        2 (1513/ 660-1566/ 761) -76.0 W27-W31 S25-S17
     10  832  710 (S10E16 L=  8)   201     45        1 ( 818/ 674- 844/ 748)  68.2 E17-E15 S12-S07
     11 1342  641 (S23W15 L= 39)   228     79        1 (1308/ 624-1387/ 662)   9.5 W12-W18 S23-S21
     12  644  310 (S36E44 L=341)  1231    203        1 ( 603/ 257- 679/ 384)  67.4 E48-E40 S41-S29
  

    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