Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  521 1707 (N46E51 L=140)  1786    373        2 ( 477/1634- 550/1765) -29.4 E60-E44 N40-N50
      2  761 1492 (N31E21 L=170)   233    113        0 ( 730/1462- 813/1522) -26.6 E24-E16 N29-N33
      3  247 1290 (N15E60 L=131)  4928    647       12 ( 172/1161- 329/1486)  81.9 E70-E53 N06-N27
      4  369 1262 (N14E47 L=144)   680    144        3 ( 335/1190- 411/1319)  69.4 E49-E44 N09-N17
      5  755 1047 (N03E17 L=174)   219     71        0 ( 726/1011- 783/1080) -45.0 E18-E14 N00-N04
      6  239  823 (S14E58 L=134)  2510    385        6 ( 190/ 643- 285/ 924)  81.0 E63-E52 S26-S07
      7 1747  831 (S07W52 L=243)   424     80        0 (1722/ 799-1762/ 863) -45.0 W49-W53 S09-S05
      8 1085  453 (S35W07 L=198)   535     86        0 (1040/ 435-1134/ 473)  -8.1 W03-W10 S36-S33
      9  720  422 (S39E21 L=170)   829    175        0 ( 641/ 396- 808/ 448)  14.9 E28-E13 S41-S36
  

    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