Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  849 1624 (N49W03 L= 98)   583     70        2 ( 827/1584- 881/1658)  51.3 W00-W05 N46-N51
      2  876 1349 (N29E02 L= 93)   458    167        1 ( 832/1284- 930/1402)  32.0 E06-W02 N26-N31
      3  428 1161 (N26E38 L= 58)  1226    365        1 ( 300/1088- 533/1216) -14.0 E50-E30 N19-N30
      4  240 1065 (N22E57 L= 38)   823    170        2 ( 198/ 991- 284/1127)  35.0 E62-E52 N18-N25
      5 1023  868 (S02E03 L= 92)   210     65        0 (1003/ 835-1046/ 900) -76.0 E04-E02 S04-N00
      6  732  665 (S09E25 L= 71)   230     42        1 ( 706/ 645- 755/ 683) -56.3 E26-E23 S10-S07
      7  526  433 (S21E49 L= 46)   485     72        1 ( 504/ 407- 559/ 456)   9.5 E52-E46 S22-S18
      8  670  360 (S28E42 L= 54)  1318    235        2 ( 608/ 271- 741/ 450)  76.0 E44-E39 S36-S20
  

    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