Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  326 1500 (N30E63 L=180)   467     58        1 ( 303/1481- 359/1522)  -5.2 E68-E57 N28-N31
      2  366 1265 (N15E48 L=195)   319     98        1 ( 351/1210- 393/1301) -71.6 E49-E46 N11-N17
      3  274  727 (S19E57 L=186)   443     94        0 ( 231/ 714- 318/ 740)   0.0 E62-E51 S20-S18
      4  294  642 (S25E58 L=186)   681    136        0 ( 278/ 578- 327/ 695)  42.0 E61-E51 S29-S20
      5  621  531 (S30E28 L=215)   377    105        0 ( 570/ 508- 669/ 552)  18.4 E33-E24 S32-S28
      6 1633  419 (S37W58 L=302)   498     87        0 (1617/ 397-1663/ 459) -24.4 W53-W65 S38-S33
      7 1436  394 (S39W37 L=280)   430    100        0 (1389/ 373-1485/ 419)  23.2 W33-W41 S40-S36
      8 1065  387 (S40W05 L=249)   429     56        0 (1038/ 364-1092/ 410)  53.1 W03-W07 S41-S37
      9  941  373 (S42E05 L=239)   657     93        2 ( 877/ 344- 984/ 396)   0.0 E09-E01 S44-S39
     10  664  362 (S44E30 L=214)   393     50        1 ( 631/ 343- 694/ 382)  -9.5 E32-E26 S45-S42
  

    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