Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  673 1628 (N38E35 L=339)  1522    283        4 ( 540/1586- 792/1668)  -4.2 E50-E23 N36-N40
      2  308 1520 (N27E68 L=307)   706     72        2 ( 277/1494- 341/1537)  -7.6 E76-E61 N25-N28
      3 1280 1384 (N27W15 L= 29)  1568    393        4 (1193/1294-1342/1505) -56.3 W08-W19 N20-N34
      4  910 1468 (N29E12 L=  3)   268     43        1 ( 882/1452- 936/1491)   0.0 E13-E09 N27-N30
      5  457 1363 (N18E43 L=332)   554    167        0 ( 401/1295- 497/1403) -28.6 E49-E38 N14-N20
      6 1737 1259 (N21W51 L= 66)   327     76        0 (1723/1218-1751/1306) -90.0 W50-W52 N18-N24
      7 1579 1272 (N21W37 L= 51)   314     61        1 (1557/1243-1604/1302) -63.4 W35-W38 N19-N23
      8  918 1120 (N07E07 L=  7)   782    245        1 ( 835/1048-1017/1182)  42.0 E11-E01 N02-N11
      9 1283  565 (S25W22 L= 36)   393     89        1 (1240/ 543-1329/ 595)  -8.1 W18-W25 S26-S23
  

    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