Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  540  924 (N00E32 L=194)   223     63        0 ( 519/ 894- 563/ 950) -45.0 E33-E30 S01-N01
      2  998  879 (S04E02 L=223)   494    202        1 ( 909/ 848-1085/ 928) -33.7 E07-W02 S06-N00
      3  997  634 (S20E04 L=222)   205    103        1 ( 980/ 601-1019/ 688) -80.5 E04-E02 S22-S16
      4  248  649 (S17E65 L=161)   545     60        0 ( 223/ 620- 268/ 674) -53.1 E68-E61 S19-S15
      5  489  613 (S20E40 L=186)   788    143        2 ( 451/ 550- 537/ 660)  33.7 E44-E35 S23-S16
      6 1732  627 (S27W54 L=280)   488     50        1 (1709/ 602-1754/ 658) -63.4 W51-W56 S28-S24
      7  430  512 (S27E50 L=176)   676     79        2 ( 404/ 485- 470/ 539) -36.9 E53-E46 S28-S24
      8 1487  458 (S36W33 L=259)   346     69        0 (1451/ 434-1515/ 484)  31.0 W30-W35 S38-S34
      9  436  425 (S33E55 L=171)   603     75        1 ( 419/ 392- 466/ 466)  36.9 E59-E51 S36-S30
  

    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