Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1397 1740 (N36W56 L=143)  2412    301        2 (1332/1644-1480/1826)   1.9 W46-W76 N29-N39
      2 1734 1052 (S12W45 L=132)   287     37        1 (1714/1028-1750/1081) -90.0 W43-W46 S13-S09
      3  456  979 (N18E36 L= 50)   329     84        1 ( 422/ 949- 490/1011) -51.3 E38-E34 N15-N20
      4 1633  902 (S17W34 L=120)  1033    224        0 (1610/ 803-1652/1011)  57.5 W30-W37 S22-S11
      5  261  794 (N10E58 L= 28)  2109    353        3 ( 200/ 652- 329/ 933) -71.6 E62-E54 N01-N19
      6 1469  555 (S34W15 L=102)   295     50        0 (1440/ 529-1491/ 583)  26.6 W12-W17 S34-S31
      7  985  433 (S27E20 L= 67)   365     50        0 ( 954/ 406-1013/ 457)  11.3 E22-E17 S28-S25
      8 1371  304 (S51W01 L= 88)   968    250        0 (1262/ 233-1460/ 368)  12.8 E11-W11 S53-S47
      9  803  250 (S35E43 L= 43)   525     75        0 ( 771/ 229- 837/ 278)  12.5 E47-E38 S36-S33
  

    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