Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  657 1724 (N26E47 L=285)   928    112        2 ( 620/1682- 699/1771) -15.3 E52-E41 N22-N29
      2  464 1704 (N21E63 L=269)  1765    151        6 ( 389/1661- 518/1746)  54.5 E69-E57 N17-N24
      3 1046 1437 (N17E09 L=323)   306     58        1 (1020/1402-1067/1465)  90.0 E10-E08 N13-N18
      4 1720 1233 (N24W39 L= 11)  1127    290        1 (1674/1103-1771/1353) -47.5 W33-W45 N18-N31
      5  923  957 (S13E04 L=328)   195     37        1 ( 903/ 933- 939/ 978)  71.6 E05-E03 S14-S11
      6 1714  642 (S08W56 L= 28)   379     64        1 (1691/ 621-1738/ 669) -26.6 W54-W59 S09-S06
      7  965  608 (S32W08 L=340)   291     66        0 ( 933/ 582- 995/ 632)  -9.5 W05-W11 S33-S30
      8 1661  516 (S16W60 L= 32)   434     36        1 (1638/ 492-1686/ 542) -71.6 W57-W61 S17-S14
      9  558  390 (S61E21 L=311)   460     61        1 ( 537/ 362- 576/ 428)  71.6 E23-E18 S63-S57
     10  720  322 (S61W01 L=333)   515     70        0 ( 696/ 293- 753/ 356) -20.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