Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1240 1739 (N42W36 L= 79)   300     50        0 (1209/1721-1267/1755)   0.0 W32-W38 N40-N42
      2  990 1591 (N37W10 L= 53)   375    158        0 ( 915/1572-1057/1608) -21.8 W04-W14 N33-N38
      3  571 1551 (N43E23 L= 20)   278     53        0 ( 540/1523- 596/1571)  21.8 E26-E21 N41-N44
      4 1636  595 (S37W37 L= 80)   353     53        1 (1599/ 582-1664/ 615)  -9.5 W33-W39 S38-S35
      5  357  469 (S20E63 L=340)  2498    221        4 ( 302/ 386- 419/ 545) -90.0 E68-E59 S25-S14
      6  543  399 (S27E49 L=354)  1869    433        5 ( 411/ 280- 683/ 485) -62.2 E56-E43 S38-S19
      7 1198  376 (S43E02 L= 41)   564    101        2 (1160/ 343-1239/ 418) -56.3 E04-W00 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