Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1455 1693 (N31W51 L=295)   353     21        0 (1434/1674-1478/1712) -56.3 W49-W53 N29-N32
      2  550 1565 (N47E22 L=222)   347     46        1 ( 519/1545- 579/1582)  -9.5 E25-E19 N46-N49
      3 1446 1544 (N23W41 L=284)   467    109        1 (1405/1502-1487/1578) -71.6 W38-W42 N19-N25
      4  413 1534 (N49E37 L=206)   572     97        2 ( 387/1491- 444/1570)  59.0 E40-E33 N45-N51
      5 1333 1521 (N24W32 L=276)   538    137        2 (1294/1472-1362/1568)  26.6 W27-W35 N21-N26
      6  931 1476 (N31W05 L=249)  1109    333        2 ( 774/1447-1090/1496) -16.7 E05-W15 N28-N34
      7 1440 1435 (N17W35 L=279)   459    189        0 (1365/1390-1506/1486)   8.7 W28-W41 N15-N18
      8 1235 1206 (N09W16 L=260)  1122    308        2 (1112/1174-1381/1237) -14.0 W08-W24 N04-N11
      9  358 1014 (N16E42 L=201)  1009    224        3 ( 269/ 931- 442/1065) -36.9 E50-E36 N10-N19
     10 1109  795 (S11E00 L=243)   334    110        2 (1076/ 747-1144/ 850) -90.0 E01-W01 S14-S07
     11 1683  665 (S34W37 L=281)   617    129        0 (1627/ 638-1736/ 698)   0.0 W31-W42 S35-S32
     12  665  613 (S12E31 L=213)   299     95        0 ( 613/ 598- 708/ 637)   0.0 E34-E27 S13-S11
     13  549  584 (S12E40 L=204)   274     87        0 ( 531/ 550- 576/ 621)  68.2 E40-E37 S14-S09
     14  661  502 (S19E35 L=208)   339     69        1 ( 642/ 463- 679/ 542)  59.0 E36-E33 S21-S16
     15  447  508 (S14E51 L=193)   420     66        2 ( 407/ 487- 481/ 527)  -9.5 E54-E47 S15-S13
     16  537  331 (S27E55 L=188)  1578    179        1 ( 463/ 287- 624/ 365)   2.5 E68-E45 S29-S24
     17  995  245 (S45E27 L=217)  1502    211        1 ( 908/ 196-1079/ 305)   3.4 E35-E18 S47-S41
  

    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