Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  465 1666 (N25E56 L=148)   493     48        0 ( 445/1641- 488/1692)  90.0 E57-E53 N23-N27
      2 1588 1588 (N41W33 L=237)   934    111        2 (1540/1549-1636/1622) -12.5 W27-W37 N38-N43
      3 1090 1430 (N19E03 L=201)   295     83        1 (1046/1410-1128/1450)   0.0 E06-W00 N17-N19
      4 1457 1413 (N25W21 L=225)   281     82        0 (1416/1390-1494/1440)  -9.5 W17-W23 N23-N26
      5 1657 1324 (N24W37 L=241)   431    126        0 (1617/1269-1694/1369) -33.7 W33-W40 N21-N26
      6  145 1099 (S13E68 L=136)  3123    181        4 (  91/1036- 185/1170)  76.0 E79-E62 S16-S08
      7 1314  761 (S16W22 L=227)   530    216        0 (1238/ 691-1371/ 841)  65.6 W19-W24 S21-S10
      8  339  803 (S30E46 L=158)  1349    217        4 ( 242/ 764- 444/ 834)   8.1 E57-E36 S32-S27
      9 1628  667 (S15W47 L=251)   803    222        1 (1576/ 622-1690/ 722)  49.4 W43-W50 S18-S10
     10 1552  528 (S25W48 L=252)  3856    702        1 (1284/ 445-1882/ 686)  31.0 W30-W76 S35-S07
     11  337  459 (S54E63 L=142)   633     77        1 ( 307/ 427- 368/ 494)  12.3 E76-E53 S55-S50
  

    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