Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  594 1676 (N26E49 L=247)  1585    445        4 ( 533/1554- 632/1786) -35.0 E60-E39 N18-N32
      2  814 1735 (N35E37 L=259)   549    105        1 ( 773/1708- 864/1761)  -6.3 E41-E32 N33-N36
      3  536 1410 (N08E40 L=256)   508    136        2 ( 495/1382- 574/1443)  -8.1 E44-E37 N05-N09
      4  235 1194 (S11E59 L=237)  1082    143        3 ( 210/1123- 268/1271)  90.0 E61-E56 S15-S05
      5  178 1130 (S16E64 L=232)   654    106        2 ( 163/1096- 198/1161) -53.1 E66-E60 S17-S13
      6  643  999 (S14E22 L=274)   439     85        1 ( 605/ 964- 683/1042)  80.5 E23-E19 S16-S10
      7  188  835 (S34E58 L=238)   557     88        2 ( 168/ 808- 214/ 864) -56.3 E60-E54 S35-S32
      8  613  756 (S29E19 L=276)   316     66        0 ( 584/ 724- 642/ 780)  51.3 E21-E17 S31-S26
      9  461  621 (S42E31 L=265)   708    143        1 ( 412/ 586- 508/ 662)  37.9 E35-E26 S45-S38
     10 1554  369 (S27W59 L=354)   468    109        0 (1524/ 353-1586/ 389)   8.1 W54-W61 S27-S25
  

    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