Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  484 1656 (N37E57 L=342)  1095    117        3 ( 431/1631- 544/1684)  -3.2 E66-E48 N35-N39
      2  598 1539 (N30E38 L=  2)   577    104        1 ( 571/1486- 642/1587) -45.0 E42-E34 N26-N33
      3  935 1414 (N25E10 L= 30)  1052    319        1 ( 845/1341- 984/1532)  85.6 E14-E05 N20-N33
      4 1753 1374 (N29W56 L= 96)  4265    602        4 (1545/1200-1923/1516) -23.2 W37-W81 N18-N38
      5  348 1369 (N17E53 L=346)   437    101        0 ( 309/1348- 382/1391)  18.4 E56-E49 N15-N18
      6 1245 1302 (N21W12 L= 52)   588    233        2 (1173/1259-1331/1359)  35.0 W07-W17 N17-N25
      7 1784 1212 (N19W56 L= 95)  1261    245        3 (1733/1095-1834/1309) -70.3 W50-W60 N11-N25
      8 1705 1151 (N15W47 L= 86)   461    123        0 (1674/1116-1748/1204) -39.8 W44-W50 N12-N17
      9 1410 1121 (N11W24 L= 63)   213     86        0 (1371/1094-1439/1147)  45.0 W21-W25 N08-N12
     10  826  734 (S18E10 L= 29)   570    206        0 ( 794/ 645- 891/ 787)  68.2 E12-E06 S23-S13
     11 1749  640 (S17W60 L= 99)  1003    139        2 (1705/ 595-1808/ 698)  41.2 W55-W64 S19-S12
  

    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