Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1039 1687 (N50W13 L=140)   592     93        1 (1011/1643-1065/1723)  45.0 W09-W16 N46-N52
      2 1268 1668 (N44W31 L=158)   357     74        2 (1245/1638-1295/1689)  21.8 W28-W33 N41-N45
      3  965 1606 (N44W04 L=131)   491     81        2 ( 933/1576- 996/1635) -68.2 W01-W06 N41-N46
      4  582 1457 (N38E28 L= 99)   286     44        1 ( 567/1434- 603/1483) -45.0 E29-E26 N36-N39
      5 1246 1050 (N05W14 L=141)   486    158        1 (1171/1021-1325/1096) -29.1 W09-W18 N02-N08
      6 1071  483 (S30E04 L=124)  1084    341        1 ( 908/ 408-1190/ 589)  20.9 E16-W05 S33-S23
      7  360  486 (S23E58 L= 69)  1095    327        0 ( 317/ 439- 412/ 525)  21.8 E66-E51 S27-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