Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  983 1767 (N38E29 L= 98)  2613    460        7 ( 857/1647-1062/1864) -21.8 E44-E19 N29-N43
      2  460 1678 (N18E61 L= 67)  3301    363        6 ( 381/1611- 538/1756)  45.0 E67-E53 N13-N24
      3  549 1667 (N19E52 L= 76)   744     98        2 ( 506/1631- 591/1709)  71.6 E54-E49 N16-N22
      4  548 1557 (N12E46 L= 82)  1425    318        6 ( 460/1469- 628/1615)  28.6 E53-E40 N07-N17
      5  609 1528 (N12E40 L= 87)   471     85        1 ( 586/1481- 640/1564) -71.6 E41-E37 N08-N14
      6  249 1501 (N03E74 L= 53)   885    129        2 ( 223/1464- 273/1531)  68.2 E79-E71 N00-N05
      7 1367 1306 (N20W12 L=139)   629    134        2 (1309/1266-1424/1344)  23.2 W08-W16 N16-N22
      8 1536  499 (S19W50 L=177)   294     77        0 (1508/ 485-1578/ 515)  45.0 W47-W52 S19-S16
      9 1152  290 (S44W39 L=166)  1187    236        1 (1049/ 262-1260/ 331)  17.4 W29-W45 S46-S39
  

    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