Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  588 1610 (N52E29 L=252)   614    110        0 ( 538/1580- 647/1631) -19.7 E35-E21 N48-N53
      2  504 1601 (N52E39 L=242)   328     24        1 ( 488/1579- 521/1622) -76.0 E42-E36 N50-N54
      3 1688 1406 (N17W55 L=336)   784    175        1 (1615/1359-1749/1461) -45.0 W50-W60 N13-N21
      4  998 1294 (N23W03 L=284)   304    122        0 ( 942/1265-1052/1323) -45.0 W00-W05 N20-N25
      5 1491 1239 (N11W33 L=314)   235     82        0 (1464/1218-1526/1269) -53.1 W31-W35 N09-N13
      6  837  805 (S04E15 L=266)   239     81        0 ( 791/ 786- 876/ 819) -11.3 E17-E12 S05-S02
      7  999  440 (S30E12 L=269)  2470    531        4 ( 790/ 316-1153/ 617)  25.9 E32-W01 S37-S20
      8 1518  463 (S40W28 L=309)   633    111        1 (1474/ 442-1554/ 496)  15.9 W24-W31 S41-S37
  

    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