Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1823 1349 (N17W65 L=351)   565    104        1 (1795/1321-1849/1368)  -8.1 W61-W69 N15-N18
      2  339 1009 (N06E47 L=238)   284     42        1 ( 318/ 983- 362/1030)  56.3 E49-E45 N04-N07
      3 1826  804 (S16W59 L=344)  1870    413        3 (1708/ 660-1907/ 964)  40.4 W47-W70 S24-S07
      4 1016  932 (S01E01 L=285)   203     43        1 ( 992/ 911-1047/ 951)  26.6 E02-W00 S02-N00
      5  437  832 (S05E40 L=246)   534     93        1 ( 407/ 793- 473/ 880) -56.3 E42-E37 S07-S01
      6 1561  584 (S28W36 L=322)   316     58        0 (1534/ 552-1588/ 612)  45.0 W34-W38 S29-S25
      7  796  531 (S26E19 L=267)   240     71        0 ( 780/ 499- 824/ 562) -68.2 E20-E17 S28-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