Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  680 1559 (N16E36 L=322)   673     96        1 ( 636/1531- 736/1584)  33.7 E39-E32 N13-N18
      2 1415 1368 (N23W15 L= 12)   995    350        0 (1351/1209-1477/1514) -47.0 W07-W21 N14-N31
      3 1672 1192 (N20W36 L= 34)  1068    299        1 (1597/1073-1740/1284) -31.6 W29-W44 N15-N24
      4  453  889 (S28E34 L=324)   252     49        0 ( 437/ 853- 466/ 919) -53.1 E35-E32 S30-S26
      5 1033  591 (S32W13 L= 11)   219     36        1 (1008/ 571-1058/ 608)   0.0 W11-W14 S32-S30
      6  925  541 (S38W08 L=  6)   414     86        0 ( 879/ 522- 976/ 561)  29.7 W04-W11 S40-S35
      7  822  455 (S47W04 L=  2)   294     56        0 ( 795/ 435- 848/ 484)  68.2 W01-W05 S48-S43
      8  397  446 (S61E48 L=310)   729     69        1 ( 360/ 419- 424/ 479)  17.5 E59-E40 S63-S57
      9  547  384 (S62E24 L=334)   541     61        0 ( 513/ 362- 584/ 412)  -5.7 E28-E18 S63-S59
  

    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