Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  655 1731 (N26E47 L=285)   936    110        1 ( 618/1691- 697/1777) -24.4 E52-E41 N23-N29
      2  464 1709 (N21E63 L=269)  2162    164        7 ( 390/1664- 526/1775)  77.5 E70-E58 N17-N27
      3  641 1605 (N18E40 L=292)   262     47        0 ( 629/1579- 655/1637) -36.9 E42-E38 N15-N19
      4 1043 1440 (N17E10 L=323)   339     62        1 (1024/1411-1064/1475) -45.0 E11-E08 N14-N18
      5 1719 1229 (N24W39 L= 12)  1170    280        1 (1673/1099-1770/1349) -47.5 W33-W45 N18-N30
      6  921  959 (S13E04 L=328)   220     39        0 ( 900/ 934- 939/ 978)  63.4 E05-E03 S14-S11
      7 1716  643 (S08W57 L= 29)   330     27        0 (1699/ 623-1736/ 668) -26.6 W54-W58 S09-S06
      8  976  602 (S32W09 L=341)   408    119        0 ( 929/ 574-1033/ 634)  -7.1 W05-W13 S33-S30
      9 1661  515 (S16W60 L= 32)   373     51        0 (1640/ 493-1683/ 546) -71.6 W57-W61 S17-S14
     10  558  391 (S61E21 L=311)   479     39        1 ( 537/ 363- 577/ 426)  71.6 E23-E18 S63-S57
     11  719  323 (S61W01 L=333)   527     75        0 ( 695/ 293- 752/ 358) -26.6 E02-W06 S62-S57
  

    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